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The Role of Transition from Laminar to 
Turbulent Flow in Fluid Mechanics 

J ^ 

By 

HUGH L. DRYDEN, Ph.D.* 
INTRODUCTION 

Fluid mechanics is one of the basic sciences which finds practical 
application in such diverse fields as aeronautics, hydraulic en- 
gineering, ventilation, chemical engineering, sanitary engineering, 
marine engineering, automotive engineering, railroad engineering, 
turbine design, and ballistics. In these many fields, flow relation- 
ships of such bewildering diversity are exhibited that at first sight 
the observed phenomena seem to have little in common. In reality 
the relationships are no more complicated than those involved in 
avoiding the failure of structures or structural parts under service 
conditions. The engineer soon learns to recognize certain char- 
acteristic types of behavior occurring in failures of a great variety 
of structures, and to study these types of failure under the simplest 
possible conditions. Thus, buckling is recognized as a type of 
failure characteristic of thin-walled structures or parts of struc- 
tures, and a study of the buckling of a simple strut under axial 
loading aids in the understanding of the buckling of more compli- 
cated structures. 

In a similar way it is possible to pick out certain characteristic 
flow phenomena which can be recognized in very complex flow 
patterns and in diverse types of apparatus and structures. The 
study of these characteristic phenomena under the simplest pos- 
sible conditions aids greatly in the understanding of the flow of 
fluids in general. Some of these phenomena will be discussed in 
this paper. 

* Chief, Mechanics and Sound Division, National Bureau of Standards. 
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TYPES OF FLOW AROUND A CYLINDER 

The principal phenomena may be illustrated by the flow about 
a very long circular cylinder immersed in a fluid and subjected to 
flow normal to its axis at successively higher speeds. It is well 
known that the fields of flow about cylinders of various diameters 
in various fluids and at various speeds are similar if the ratio of 
the product of diameter D by speed U to the kinematic viscosity v 
is the same. This ratio TJD/v is called the Reynolds Number, and 
is the only significant parameter provided that D is not comparable 
with the mean free path of the molecules of the fluid, that U is not 
comparable with the speed of sound in the fluid, and that the 
cylinder is not near the free surface of a liquid. For a cylinder one 
inch in diameter in air at normal pressure and at a temperature of 
78°F, a Reynolds Number of i corresponds to a speed of 0.002 
ft/sec, 100 to 0.2 ft/sec, 10,000 to 20 ft/sec, etc. For the same 
cylinder in water at 90°F a Reynolds Number of i corresponds to 
10*“^ ft/sec, 100 to 10-2 ft/sec, 10,000 to i ft/sec, etc. It is hardly 
possible to cover experimentally a range in Reynolds Number 
from I to 300,000 using one cylinder because of various practical 
difficulties, yet it may be helpful to consider an increasing Reyn- 
olds Number as corresponding to an increasing speed. 

Fig. I is a concise summary of the known facts. Six general 
types of flow are distinguished and illustrated by sketches and 
graphs of pertinent measurements. The sketches for A, E, and F 
are schematic only, but the others are taken from published ex- 
perimental results. 

The six types of flow correspond to definite characteristics in 
the behavior of the total force exerted on the cylinder by the fluid 
stream. The force is usually expressed in terms of the drag coeffi- 
cient, C^, defined as the ratio of the force to the product of pro- 
jected area (diameter times length) and velocity pressure (one- 
half the product of density by square of the speed) . The variation 
of the drag coefficient with Reynolds Number is shown in Fig. 2 
with logarithmic scales, the solid curve corresponding to a fluid 
stream of reasonably small initial turbulence and the dotted curve 
to a stream of large initial turbulence (i).* The turbulence of the 
stream referred to here and elsewhere in the paper is a small-scale 


• A number in parentheses refers to the list of references at the end of the paper. 



FROM LAMINAR TO TURBULENT FLOW 


3 



Fig. I. Types of flow around an infinitely long circular cylinder. 

A. At Reynolds Numbers less than i. 

B. At Reynolds Number of 20. S is the separation point of the laminar boundary 
layer. 

C. At Reynolds Number of 1 74. One stage of the flow which varies with the time 
as eddies are shed alternately. Traverse with total head tube along AA' shown at 
the right, with width of free laminar layer indicated. 

D. Boundaries of free layer at Reynolds Numbers of 5000 and 14480, and graphs 
of variation of width of free layer with distance X downstream. D is the diameter 
of the cylinder. The solid curve in the graph is for a Reynolds Number of 5000; 
the dotted curve is for a larger Reynolds Number or for a stream of greater turbu- 
lence at a Reynolds Number of 5000. 

E. Location of separation of laminar boundary layer and general character of 
flow at Reynolds Number of 80000. Transition occurs simultaneously with separation. 

F. Location of separation of turbulent boundary layer and general character of 
flow at Reynolds Number of 1000000. Transition occurs in the laminar boundary 
layer before separation. 
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irregular random motion superposed on the mean flow of the 
fluid by honeycombs, grids, or other upstream obstacles. The 
fluctuations of speed are of fairly high frequency and amount to 
only a few per cent of the mean speed at the most. A turbulence 
of only a few tenths of a per cent of the mean speed often has a 



Fig. 2. Variation of drag coefficient of infinitely long circular cylinder with the 

Reynolds Number. 


pronounced effect on the flow pattern. The curves shown corre- 
spond to intensities of turbulence of the order of i and 3 per cent. 

Fig I A shows the general character of the flow at extremely 
low Reynolds Numbers, i.e., less than i. The inertia forces are 
small in comparison to the viscous forces in the region near the 
cylinder, and the flow closes in smoothly behind the cylinder. The 
flow resembles superficially the potential flow of a perfect fluid. 
The stream lines, however, for the same flux are at a greater dis- 
tance from the cylinder than in the case of potential flow. The 
frictional effects extend to a distance of many times the diameter 
of the cylinder (2). 

The drag coefficient is very large but decreases rapidly with 
increasing Reynolds Number. The force itself varies more nearly 
as the first power of the speed than as the second. 

As the Reynolds Number is increased, the inertia forces become 
comparable with and finally larger than the viscous forces, and the 
character of flow gradually changes. Fig. iB shows the flow at a 
Reynolds Number of 20 as computed by Thom (3) from the 
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Navier-Stokes equations of motion of a viscous fluid. Photographs 
by Thom, by Nisi, and by Page (4) of the flow of air around a 
cylinder at Reynolds Numbers from 8 to 170 confirm the presence 
of the two stationary symmetrical eddies which become unstable 
as the Reynolds Number is increased, first vibrating irregularly, 
then breaking away alternately from the two sides. In the neigh- 
borhood of the point S occurs tfie first of thp characteristic phe- 
nomena — that called separation. The main flow separates from 
the surface of the cylinder. Close to the cylinder the direction of 
the flow reverses on passing the neighborhood of S. At these 
Reynolds Numbers the rate of decrease of the drag coefficient be- 
comes less rapid. 

At a Reynolds Number only a little above 20, the exact value 
apparently depending on various factors such as the initial turbu- 
lence and the distance from the cylinder to other objects, the 
eddies begin to break away alternately on either side in a periodic 
manner to form the well-loiown K^rmdn vortex trail (5). Fig. iC 
shows one stage of the flow obtained from one frame of a motion 
picture taken at a Reynolds Number of 174 by Nayler and Frazer 
(6). It would require many sketches to show the complete flow 
since the flow varies with the time. At this Reynolds Number the 
effect of viscosity is largely confined to a region near the cylinder 
about one-tenth of the diameter of the cylinder in thickness. This 
region is called the boundary layer. The speed decreases to zero 
at the surface of the cylinder. The thickness of the boundary layer 
decreases with increasing Reynolds Number. 

If a traverse is made with a pitot tube behind the cylinder as 
indicated by AA' in Fig. iC, the reference pressure being the pitot 
pressure in the approaching fluid stream, results are obtained (7) 
as shown in the auxiliary figure accompanying Fig. iC. The 
falling total pressure in the wake of the cylinder is associated with 
a decrease in the speed of flow. The rapid changes are confined to 
a limited region whose width may be defined as indicated. The 
width of this mixing region, which is a continuation of the bound- 
ary layer beyond the separation point and accordingly designated 
a “free” layer, increases slowly with distance from its point of 
origin at the separation point S. While the flow is not steady with 
time, the flow in this free layer as well as in the attached layer is 
essentially of the type described by the Navier-Stokes equations 
in that the lateral mixing is due solely to the action of viscosity, 
the unsteadiness being in the nature of a wavering of the layer as 
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a whole. A type of motion in which lateral mixing is due to viscous 
action is usually termed laminar, and hence both the free layer 
and the attached layer are said to be laminar layers. 

At Reynolds Numbers from 5,000 to 1 5,000 a new phenomenon 
was observed by Schiller and Linke (7). The thickness of the free 
layer after increasing slowly with distance from the separation 
point for some distance begins to increase at a more rapid rate, 
and as the Reynolds Number is increased the place at which the 
rapid increase in thickness begins moves closer and closer to the 
cylinder. Fig. iD shows the boundaries of the free layer for Reyn- 
olds Numbers of 5,000 and 14,480 as determined by Schiller and 
Linke. We have learned to associate this rapid increase in thick- 
ness with transition from laminar to turbulent flow, i.e., from a 
type of motion in which the lateral transfer of momentum is by 
viscous shear to a type in which the lateral transfer occurs by the 
irregular motion of small portions of fluid. The change in char- 
acter of flow in the wake modifies the pressure distribution in such 
a way that the drag coefficient increases with increasing Reynolds 
Number. 

Schiller and Linke made the observation that the transition 
occurs closer to the cylinder at a given Reynolds Number when 
the turbulence of the stream was increased by inserting a wire 
screen of small mesh. Under the same condition the rise in the 
drag coefficient occurred at a lower Reynolds Number. This sen- 
sitiveness to the effect of turbulence is also a well-known char- 
acteristic of the transition from laminar to turbulent flow. 

At some value of the Reynolds Number between 5,000 and 
20,000, the exact value being dependent on the initial turbulence, 
the point of transition in the free layer reaches the laminar sepa- 
ration point and for a considerable range of Reynolds Number 
transition begins at the separation point. The drag coefficient is 
constant throughout this range. The location of the separation 
point and the general nature of the flow are indicated schemati- 
cally in Fig. lE. Periodic fluctuations of speed of the type usually 
associated with a Kdrmdn vortex trail are found in the wake. 
However, no K^rm^n trail has been seen or photographed at 
Reynolds Numbers above 2,500, i.e., above those at which the 
trail results from the rolling up of a laminar layer. When turbu- 
lence sets in, the vorticity is rapidly diffused as is any substance 
used as an indicator in photographing the flow. Vortices, however, 
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break away periodically for Reynolds Numbers up to 100,000 or 
more. 

When the Reynolds Number reaches 100,000 or more, the value 
again being dependent on the initial turbulence, transition begins 
in the boundary layer ahead of the laminar separation point. 
Separation is then delayed beegwse of the more intense mixing 
which sweeps away a large part of the wake. Separation of the 
turbulent boundary layer still occurs, though delayed. The peri- 
odicity in the wake disappears and the whole wake becomes turbu- 
lent, the motions being random and highly irregular. The drag 
coefficient decreases to a much lower value. The flow is indicated 
schematically in Fig. iF. 

CHARACTERISTIC PHENOMENA 

The important characteristic phenomena which are noted in 
the flow about a cylinder and which can be found in many other 
flow patterns may be listed as follows: 

1. Formation of free and boundary layers^ i.e., regions in which in- 
tense changes in velocity are concentrated, the velocity gradients 
being much larger than elsewhere. At high Reynolds Numbers 
boundary layers are found on the upstream portions of the sur- 
faces of immersed solid bodies and often at other boundary sur- 
faces. A free layer may occur elsewhere in the flow, usually as a 
continuation of a boundary layer attached to a surface. It can be 
recognized by the presence of velocity gradients much greater 
than elsewhere in the flow. 

2. Laminar flow, i.e., a flow in which the lateral transfer of mo- 
mentum occurs solely through the action of viscosity. It can be 
recognized by the small rate of diffusion of material particles, 
streams of smoke or dye maintaining their identity for relatively 
long distances, by the relatively slow obliteration of velocity gradi- 
ents, and the relative thinness of the free layers and boundary 
layers. The average flow conforms to that computed from the 
Navier-Stokes equations. 

3. Turbulent flow, i.e., a flow in which the lateral transfer of mo- 
mentum occurs by the irregular motion of small portions of fluid. 
It can be recognized by the presence of irregular and random 
velocity fluctuations of relatively high frequency, by the rapid 
rate of diffusion of material particles or heat, and by the relatively 
thick free layers and boundary layers. 
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4. Transition^ i.e., the change from laminar to turbulent flow 
occurring in some limited part of the field of flow. It can be recog- 
nized by the change in the rate of thickening of free or attached 
layers, by a marked increase in the rate of diffusion of material 
particles and heat, by a marked increase in frequency and loss of 
periodic characteristics in the velocity fluctuations, and by an 
acceleration of the flow close to the boundaries of immersed 
objects. 

5. Separation^ i.e., a separation of the main flow from a boundary 
accompanied by a reversal in direction of the flow very close to 
the boundary behind the separation point and by the formation 
of a wake in which the velocity is much reduced. It is recognized 
most easily by the reversal which can be detected by emitting 
suitable tracer substances from openings in the boundary. 

It should be noted that transition may occur either in a free 
layer or in an attached boundary layer, but that it occurs more 
readily in free layers. So far as is known, turbulence in streams of 
constant destiny arises only in free or attached layers. 

Likewise, separation may occur either with a laminar or a 
turbulent layer, but occurs more readily in laminar layers. 

THE STUDY OF TRANSITION 

The most common occurrence of transition in the experience of 
the average person is that in the rising column of smoke from a 
cigarette lying on an ash-tray in a quiet room. The smoke rises for 
some distance in smooth filaments which may wave around but do 
not lose their identity. This flow is a laminar flow. At some distance 
above the cigarette, the filaments suddenly break up into a con- 
fused eddying motion, a turbulent motion. The transition moves 
closer to the cigarette when the air of the room is disturbed. 

The simplest flow in which transition has been observed and 
studied quantitatively is that in the boundary layer of a thin 
smooth flat plate set parallel to the fluid flow, the so-called skin- 
friction plate (8). It cannot yet be said even in this simple case 
that transition is fully understood, although much progress has 
been made. Techniques have been developed not only to detect 
the occurrence of transition by the rapid thickening of the layer 
and by the increase of the velocity close to the plate, but also to 
measure the distribution of the mean velocity and the intensity 
of the velocity fluctuations. The experimental data at hand (9) 
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show that in this case the principal parameters which determine 
the location of transition are the speed and viscosity of the fluid 
and the intensity and scale of the initial turbulence. 

If the plate is placed in a stream in which the static pressure in- 
creases or decreases, the location of transition depends on the sign 
and magnitude of the pressure gradient. A pressure decreasing 
downstream is favorable to the maintenance of laminar flow, while 
a pressure increasing downstream accelerates transition. 

Experiments under other conditions indicate that curvature of 
the surface and surface roughness also influence transition, con- 
cave or rough surfaces causing transition further upstream, while 
convex surfaces show transition further downstream than plane 
surfaces. 

The attempt to generalize the quantitative information now 
available has not been fully successful. Transition does not seem 
to depend solely on the local condition of the boundary layer as 
described by its thickness, the speed of the fluid stream just out- 
side it, the intensity and scale of the turbulence just outside, the 
local pressure gradient, and the curvature and roughness of the 
surface. The subject is under active study and it is hoped that at 
least the simplest case of transition will soon be understood. 

As an illustration of a more detailed study of transition and the 
recognition of the several characteristic phenomena, Fig. 3 sum- 
marizes the results of a survey of the speed distribution in the 
boundary layer of a cylinder of elliptic cross section with major 
and minor axes 11.78 and 3.98 inches, respectively, subjected to 
an air flow parallel to the major axis at a mean speed of 70 ft/sec 
with an initial turbulence amounting to 0.85 per cent of the mean 
speed. The curves in the figure are contour lines of equal mean 
speed derived from measurements with a hot-wire anemometer 
which is insensitive to the direction of air flow. 

The interpretation and a full discussion of these measurements 
are given in Technical Report No. 652 of the National Advisory 
Committee for Aeronautics. The boundary layer is laminar or 
nearly so to a distance of about 10.8 inches from the nose. Separa- 
tion of the laminar boundary layer occurs at about 10 inches from 
the nose as compared to a value of 9.5 inches computed from the 
theory of laminar boundary layers. Transition of the free layer is 
evidenced by the rapid approach of the contours toward the sur- 
face beginning at a distance of 10.8 inches from the nose. In this 
region the speed near the surface increases on proceeding down- 
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stream. The layer in fact reattaches itself to the surface and sepa- 
rates again as a turbulent layer at 12.1 inches from the nose. 



Fig. 3. Distribution of mean speed in the boundary layer of a 3.98 by 1 1 .78 inch 
elliptic cylinder in an airstream speed 70 ft/sec, intensity of turbulence 0.85 per- 
cent, scale of turbulence 0.26 inch. 

Uo is the speed at a large distance from the cylinder, u is the local speed. The 
curves are contours of constant speed. 

Laminar separation occurs at a; 10, transition of the free laminar layer at 10.8 
followed by reattachment of the layer and separation of the turbulent layer at 12.1. 

ROLE OF TRANSITION 

In general terms transition may be stated to play both beneficial 
and detrimental rdles in fluid mechanics. It is beneficial in that it 
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delays separation of the flow and promotes more rapid diffusion 
of momentum, matter, and heat. It is detrimental in that it in- 
creases skin-friction and often modifies the pressure distribution 
in such a way that drag is increased. These general effects may 
be observed in practically every field of application of fluid 
mechanics. j ^ 

In aeronautics, for example, the lift of the wings of an airplane 
is very much reduced and their drag greatly increased if separa- 
tion of the flow occurs on the upper surface. The occurrence of 
transition is beneficial in that it causes the air flow to close in 
around the wings delaying separation until close to the trailing 
edge. The lift of the wings is maintained to relatively high angles 
of attack and the drag is much less than would be the case with 
separation of a laminar boundary layer. The desired performance 
of wings, struts, and fuselages could not be obtained if separation 
were not delayed by the transition in the boundary layer from 
laminar to turbulent flow. On the other hand, transition is detri- 
mental in that it increases the skin-friction of wings, fuselages, and 
propeller blades. High-speed performance is greatly affected by 
the location of transition on these members. In aeronautical en- 
gineering it is desirable therefore to be able to control the position 
of transition so that it is delayed as long as possible to avoid in- 
crease in the frictional resistance but that it occurs early enough 
to avoid flow separation. 

Another important r61e of transition in aeronautics is in its re- 
lation to the interpretation of model experiments. Small-scale 
turbulence is absent in the atmosphere but is present in varying 
degree in wind tunnels. Because of the effect of turbulence on 
transition great care must be exercised in the conduct and inter- 
pretation of wind tunnel tests. 

In marine, automotive, and railroad engineering transition 
plays the same r61e in determining the resistance to forward mo- 
tion as in aeronautics. Towing basins as well as wind tunnels have 
varying degrees of initial turbulence, and it is the practice in many 
towing basins to control the position of transition on ship models 
by adding local surface roughness in the form of sand grains or’ 
wire bands at the desired transition point. 

In chemical engineering transition plays an important part in 
problems of the flow of fluids in pipes or other enclosed channels, 
in problems of diffusion and mixing, and in problems of heat 
transfer. The relationships are very different in laminar and turbu- 
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lent flow. The occurrence of transition increases pressure losses in 
flow in pipes, increases the rate of diffusion of matter and the rate 
of transfer of heat. 

In hydraulic engineering, transition plays a part in the transport 
of sand and sediment, the scouring of river beds, the separation of 
the flow around obstacles, around the blades of water turbines, 
and in draft tubes and other diffusers, friction in pipes and chan- 
nels, and various mixing processes such as those involved in the 
disappearance of density currents. 

Aiiy of these applications could be described in some detail, al- 
though in many cases the concepts outlined herein have not been 
applied to guide the experimental work. Special techniques are 
required to identify the occurrence of the characteristic phe- 
nomena, but it is by no means necessary to make tedious velocity 
traverses very close to the surface as illustrated by Fig. 3. Such 
traverses supply a great deal of additional information but simple 
methods are sufficient for mere identification of separation and 
transition. 
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Problems of Flow in Compressible Fluids 

By 

J * 

THEODORE VON KARMAN, Ph.D., Dr* Ing., Sg.D.* 

In many applications of fluid mechanics the assumption of incom- 
pressibility or constant density of the fluid yields results of sufficient 
accuracy for practical purposes. In general it is assumed that this 
approximation is justified if the velocity, density, and pressure 
changes are relatively small, as, for example, in most problems of 
practical hydraulics and aeronautics. Often it is stated that the 
dynamics of incompressible fluids will give satisfactory results if 
the velocities involved in the problem are small in comparison 
with the velocity of sound. However, if we look over the entire 
field of fluid mechanics somewhat more carefully, we find that 
besides the high-speed phenomena there are other cases in which 
density variation or elasticity of the fluid cannot be neglected. In 
all, the problems in which compressibility enters as a governing 
factor can be classified in the following groups: 

a) Pressure propagation. In an incompressible fluid the pressure 
changes propagate with infinite velocity. Hence, if we are con- 
cerned with pressure oscillations in fluids, we have to consider the 
elasticity of the fluid as the factor governing the phenomena, irre- 
spectively of the magnitude of the velocities or the pressure changes 
involved. Problems of acoustics, water hammer in pipes, oscilla- 
tions in gas conduits, engine manifolds, etc., belong in this class 
of problems. 

b) Density currents. In the problems belonging in this group we 
are concerned with a stratified fluid medium, as the atmosphere 
or water carrying silt. In such problems the velocities and pressure 
changes caused by the motion may be small, but the phenomena 
are governed by the gravity field produced by density differences. 

c) Slow motion of fluids with large density changes. In this class of 
problems the pressure and density changes may be produced by 
some mechanism, for example, by the pistons of a reciprocating 

* Director of the Daniel Guggenheim Aeronautical Laboratory and Professor of 
Aeronautics, California Institute of Technology. 
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engine. In other cases the change of density is caused by large 
pressure gradients, as, for example, in the flow of gases through 
granular media. 

d) High-speed flow and high-speed motion of solids in fluid media. In 
flow problems belonging in this class the density and pressure 
differences are produced or accompanied by large variations of 
the velocity. This class includes the flow of elastic fluids through 
ducts, vane systems, etc., with velocities of the order of or higher 
than the velocity of sound, and the motion of solid bodies with 
such velocities. Such problems occur in the design of gas and 
steam turbines, compressors, propellers, jet propulsors, high-speed 
airplanes, and in particular in exterior ballistics. 

This paper is restricted to a short review of a few problems be- 
longing to item d) . They deserve a great deal of attention because 
a better understanding of the high-speed phenomena is essential 
for further progress in machinery dealing with elastic fluids, in 
high-speed aeronautics, and in shell and bomb design. 

We shall use the term “gas dynamics” for this group of fluid 
mechanics problems. It is understood that gas dynamics is based 
on combination of the principles of fluid mechanics and thermo- 
dynamics. 

I. GAS DYNAMICS IN “HYDRAULICS” FASHION 

It has become fashionable in the last few years to make a dis- 
tinction between hydraulics and hydrodynamics (or fluid me- 
chanics). It is hard to say exactly where hydraulics ends and fluid 
mechanics begins. Some authors consider hydraulics as a purely 
empirical science and hydrodynamics as a purely mathematical 
discipline dealing with ideal fluids. I believe this is a rather errone- 
ous view, since hydraulics is based on the same principles of 
mechanics which are the essential content of the equations of 
hydrodynamics, and hydrodynamics has never restricted its do- 
main to ideal fluids. It seems to me that what we might call treat- 
ment of fluid mechanics problems in the “hydraulic manner” is 
something similar to the approximate treatment of problems of 
elastic bodies by the ordinary beam theory. In both cases field 
distributions, i.e., two- and three-dimensional velocity and stress 
distributions, are approximated by essentially one-dimensional 
patterns, using certain standard distributions and average values 
across the cross sections of a pipe, channel, or beam. 
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Let us consider, for example, the flow of a fluid through a nozzle 
of variable cross section. To fix the ideas we assume that the fluid 
is at rest in a container under certain pressure and flows into a 
space of zero absolute pressure. The mass flow through the nozzle 
is given by the product pAv where A is the cross section area, v the 
mean velocity of the fluid through the cross section, and p its den- 
sity. Hence, the condition of continuity requires that pAv = const. 
In the case of an incompressible fluid the density is constant and 
therefore we have Av = const, or 


dA ^ 

1 ~ “7 

In the case of an elastic fluid 

dA dv ^ 

A V p 


(I) 


( 2 ) 


If we neglect friction losses, the dynamic equation for stationary 
flow can be written in the form {p = pressure) 

dp dp 

vdv + -j- — — o 
dp p 


or 


P 


dp 

dp 


vdv 


( 3 ) 


It is known that the quantity ~ is equal to the square of the 

dp 

velocity of propagation of infinitesimally small pressure difference, 
in other words it is equal to the square of the velocity of sound a 
in a gas of the given pressure p and density p. Hence 


dp V 
— ^dv 


p 

Substituting Eq. (4) in Eq. (2) we obtain 

U 

A 




or introducing the ratio - as the “local Mach’s number” Af, 


It 


) 


( 4 ) 


(5) 


( 6 ) 


is seen 
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1 ) that the approximation of incompressibility holds if i . 

2) that the cross section corresponding to a stationary flow 
through the nozzle decreases with increasing velocity if M^di 
and increases with increasing velocity if > i . 

3) that at the cross section of minimum area, i.e., “at the 
throat,” M = 1; the velocity is equal to the “local velocity of 
sound.” 



Fig. I 


Figure i shows the distribution of cross section area as function 
of the velocity v first for an incompressible fluid (line a) then for 
an elastic fluid (line b) under the assumption of adiabatic change. 
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In this case where 7 = ^ is the ratio between 

2 Ct, 

the specific heats for constant pressure p and constant volume F; 

is the velocity of sound corresponding to the initial state at rest. 

Let us now assume that a certain nozzle has been designed for a 

pressure ratio ~ where is the ratio between initial pressure 

p9 

and the pressure at the throat. The question arises what happens 
if the pressure is maintained at the entrance of the nozzle, but 
the pressure at the exit is varied between the limits p^ and p^. 

First, it is evident that the fluid, after reaching the velocity of 
sound at the throat, can either continue the isentropic expansion 
or can be isentropically compressed in the diverging portion of 
the nozzle. Evidently the flow capacity of the nozzle is the same 
in both cases; however, in the second case the gas leaves the exit 
with a velocity vj < a. This velocity and the pressure at the exit are 
the same as in the cross section of the same area in the converging 
portion of the nozzle. Denoting this pressure by pj we see that we 
obtain smooth flow with maximum capacity for two values of the 
pressure, namely p^ and pj. If the exit pressure is higher than />/, 
i.e., it is between the limits p^ and pj^ we obtain subsonic isen- 
tropic flow whose capacity increases with the pressure difference 
between entrance and exit. If the exit pressure is between pj and 
p^, no isentropic flow is possible that is consistent with both the 
pressure difference and the exit cross section area. In such cases 
a more or less discontinuous change of pressure takes place some- 
where between the throat and the exit or in the free jet behind 
the exit. Such discontinuous change is called a “shock wave” for 
reasons to be explained in the next section. 

If we assume that the discontinuous pressure change takes place 
at a certain cross section and behind this section the process is 
isentropic again, we can compute isentropic flows with new values 
of the entropy so that they are consistent with the given end pres- 
sure and the exit cross section area. Figure 2 shows a family of 
pressure curves corresponding to such isentropic flows. Figure 3 
shows the corresponding states of the fluid in the pressure volume 
diagram. If the discontinuous transition ends on the upper branch 
of one of these pressure curves, the gas is compressed and its flow 
is subsonic downstream of the shock; a jump ending on a lower 
branch would mean that the gas expands and flows with super- 
sonic velocity. It will be seen later that the second case is physi- 
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Fig 3 


V 
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cally not possible. In fact, there is a minimum value pj' of the end 
pressure, which can be obtained by such a shock; this is the value 
obtained when the shock takes place at the exit. If the end pressure 
is smaller than the shock occurs in the free jet in the form of 
so-called “oblique” shock waves. It can also occur that the stream 
separates from the wall. Which process takes place depends prob- 
ably on the shape of the nozzle and frictional effects. If th 6 end 
pressure is between pj and /?/', a plane shock wave is observed in 
the diverging part of the nozzle. The magnitude and location of 
this discontinuity depends on the dynamical relations to be dis- 
cussed below. 


2. SHOCK WAVES 

The velocity of sound in a fluid is the velocity of propagation 
of infinitesimally small pressure changes. It has been mentioned 
that the square of the velocity of sound is equal to the derivative 
of the pressure p with respect to the density p. In an ideal case 

under the assumption of isentropic change “ = 7-, where 7 is 

dp p 

the ratio between the specific heats Cp and c^, Eamshaw and Rie- 
mann investigated the propagation of finite disturbances. In par- 
ticular Riemann emphasized that in a compressible fluid finite 
discontinuities of pressure can be developed by a process which 
itself is continuous. 

Let us consider one dimensional wave motion in an elastic fluid 
originally at rest. The dynamic equation can be written in the 
form 

du du I dp 
or 

du du dp I dp 
The equation of continuity is given by 

^+£(p«) = o (2.3) 

If we consider small perturbations, i.e., we assume that p' = p — po 


(2.1) 

(2.2) 
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(po = density at rest) and u are small quantities, and put = a*, 

dp 

we obtain by neglecting terms of higher order: 

^ _ d^ ^p 
dt Po dx 

dp' du 

dt ‘"’‘dx 

It is easily seen that for example, 

u == Uo<t>{x — at) 

p' = Po~<l>{x - at) (2.5) 

where is an arbitrary function, is a solution of this system of 
differential equations. The equations (2.5) represent a plane wave 
of arbitrary shape proceeding in the positive x direction with the 
velocity of propagation a. 

It can be easily verified that the wave motion represented by 
equations (2.5) satisfies the equation: 


and that — = -• Now if we assume that also in the case of a finite 

Po a 

disturbance the variation of the density is a function of the velocity 
only, it can be shown that equation (2.6) will be replaced by 

du du 

(2.7) 

Consequently, if the distribution of u{x) is given by the first equa- 
tion (2.5) at ^ = o, for example, the shape of the velocity distribu- 
tion curve will be deformed for ^ > o. Whereas in the case of 
infinitesimally small waves the velocity distribution curve is dis- 
placed unchanged by a length aAt in the time A/, in the case of 
the finite wave the displacements of the positive ordinates are 
larger and those of the negative ordinates are smaller than aAL 

Hence, the velocity gradient ^ and the density gradient ^ in- 
crease during the propagation of the wave. If we apply this con- 
sideration, for example to a wave of sinusoidal shape [i.e., if 
(a: — at) is a sine function], we find that the tendency to develop 
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an abrupt change in density occurs at places where the fluid 
moves in the direction of increasing density and pressure. After 

a certain time — becomes infinite and thereafter the equations 

fail to give a continuous solution for the velocity and density 
distributions. 

The process is somewhat analogous to the collapse of ocean 
waves breaking on the shore due to the fact that the velocity of 
propagation of the wave crest is greater than that of the trough. 

Assuming now that pressure, density, and velocity undergo a 
discontinuous change at a plane x = constant, the velocity of 
propagation of the discontinuity can be computed by satisfying 
the following conditions: (a) conservation of matter through the 
discontinuity surface, (b) conservation of momentum, and (c) 
conservation of energy. These relations are established in the 
easiest way if the discontinuity surface is considered stationary 
and the fluid passing through it. In other words we consider the 
discontinuity surface as the front of a standing wave known as a 
shock wave. In this case pressures and densities on the two sides 
of the shock wave satisfy the following relation: 


^ = A Z A = yh±ji (2.8) 

Ap P2 — Pi P2 + Pi 

(The subscript i refers to the state before, 2 behind the shock 
wave.) 

It is seen that if pi-->p2— >p and this relation becomes 


dp ^ p 


(2.9) 


This expression is the differential equation for isentropic change. 
However, if and Ap are finite, the entropy of the gas is larger 
on the side of the shock wave where the pressure is larger. There- 
fore according to the second law of thermodynamics the gas passes 
through the shock wave from the low-pressure to the high-pressure 
side. In other words the progressing wave with finite discontinuity 
is always a compression wave: it propagates from regions of high 
pressure to regions of low pressure. It can also be shown that the 
velocity of propagation of a finite discontinuity relative to the 
fluid on the low-pressure side is always greater than the velocity 
of sound in that medium. Therefore, a standing shock wave can 
exist only if the fluid upstream of the shock wave moves with 
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supersonic velocity. The velocity of propagation relative to the 
high pressure side is smaller than the velocity of sound in that 
medium. In other words the fluid passing through a standing 
shock wave changes its condition from supersonic to subsonic. 
Thus it is clear why no plane shock wave can occur between the 
initial pressure curve and one of the lower branches in the nozzle 
considered in the preceding section. 

For a real fluid no mathematical discontinuity can exist. The 
thickness of the transition layer has been estimated by several 
authors, for example, by L. Prandtl, Lord Rayleigh, R. Becker, 
and G. I. Taylor. It appears that it is determined chiefly by the 
viscosity and the heat conduction of the gas. The calculations of 
the authors mentioned give for large compression ratios ex- 
tremely small values for the thickness, values of the order of 
10-6 — io “7 cm, so that it is doubtful whether in this case the 
hydrodynamic equations and the equations of heat conduction 
can be applied inside of the transition layer. For pressure ratios 
near to unity, i.e., for the case that the change of velocity. Aw, of 
the gas passing through the shock wave is small in comparison to 
the velocity with which the gas enters the shock wave, the thick- 
ness of the transition layer can be expressed, according to G. I. 

Taylor, by the simple relation 5 = const — where the constant 

Aw 

is of the order 6.5 and v is the kinematic viscosity of the gas and 

it is assumed that the ratio ^ (jit = coefficient of viscosity, 

Cp = heat capacity, X = coefficient of heat conduction) is con- 
stant and is slightly lower than unity. For example, for air of 
atmospheric pressure and for a velocity change of Aw = 30 m/sec, 
8 = 0.003 

These calculations are based on the assumption that the statis- 
tical equilibrium of the gas is reached instantaneously, i.e., we 
can assume a statistical equilibrium no matter how fast the change 
of state takes place. Now the time in which the gas passes through 
the transition layer is, for example, in the case considered above, 
of the order of io~* sec if the velocity of the gas is of the order of 
the velocity of sound. For larger compression ratios and velocities 
the time will be of the order of 10“^® sec. This time interval appears 
to be too small to secure statistical equilibrium. In particular the 
energy distribution of the intermolecular degrees of freedom will 
have considerable time lag, which is, for example, of the order of 
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io“® sec for nitrogen. It appears that an adequate analysis of the 
structure of the shock wave requires a treatment from the point 
of modern statistical mechanics. 


3. THE ANALOGY BETWEEN OPEN CHANNEL FLOW 
AND FLOW OF EJ.^STIG FLUIDS 


The phenomenon of the shock wave is to some extent analogous 
to the hydraulic jump observed in open channel flows. The hy- 
draulic jump is an abrupt transition between so-called super- 
critical and subcritical flow. The flow of water in an open channel 
is supercritical if the mean velocity v is larger than the velocity 
of long waves in a shallow channel, which is given by \/ gh 
{g = acceleration of gravity, h = depth of the flow). The flow is 
called subcritical if the mean velocity is under this critical value. 
The continuity equation for stationary flow has the form 


= constant (3.1) 

This corresponds to the continuity equation of an elastic gas 
through a pipe of constant cross section: 


vp = constant (3.2) 

if we let the density p of the gas correspond to the depth h of the 
flow. 

The dynamic equation for channel flow is 

gdh + = o (3.3) 

The dynamic equation for a gas: 


dp dp 
dp p 


+ vdv — o 


(34) 


We remember that ^ is the square of the velocity of sound; hence 
dp 

^ corresponds to gh^ and replacing p by A we see that ^ ~ will 
dp dp p 

be replaced by gdh. 

It is seen that both the continuity equation and the dynamic 
equation are analogous, if we consider h and p as corresponding 

, o}^ 

quantities. The quantity F = which is the resultant of the 

static pressure acting on a cross section of the fluid, corresponds 
to the pressure p of the gas. 
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The analogy between the one-dimensional flow of an elastic 
fluid and open-channel flow of an incompressible fluid can be ex- 
tended to the two-dimensional case. We find that the equations 
of the strictly two-dimensional flow of an elastic fluid and the 
approximate equations for the motion of a layer of an incompres- 
sible heavy fluid of variable thickness are identical if we replace 
the density of the gas p by the thickness h of the incompressible 

pfj^ 

fluid layer and the pressure of the gas by the quantity P = ^ — 

2 

However, assuming adiabatic change, the relation between the 
pressure and density of the gas has the form P = const. pT, whereas 
for the fluid layer P = const.A^. Hence, the flow of the fluid layer 
corresponds to the flow of a gas with 7=2. Nevertheless, the 
analogy is quite useful for quantitative study of some complicated 
flow problems. This analogy was first given by E. Jouguet; it was 
developed by D. Riabouchinsky and E. Preiswerk. The present 
writer applied the same analogy to illustrate wave phenomena in 
curved channels. 

4. PROBLEMS OF SUBSONIC FLOW IN AERONAUTICS 

One of the most surprising developments in modern fluid me- 
chanics is the successful application of the theory of potential 
motion of ideal incompressible fluids to the actual flow of air 
around airfoils and streamlined bodies. A few decades ago the 
potential theory was considered as a purely mathematical disci- 
pline. Later it was understood that it might be a valuable guide 
for understanding of the general laws which govern the lift, mo- 
ment, induced drag, etc., of airfoils. However, as aerodynamical 
design became more refined, it was found that the quantitative 
predictions of the theory are of a remarkable accuracy in spite of 
its idealized conception. As a matter of fact, the aeronautical 
designer has today a great confidence in the results of the theory 
of potential motion. It is possible to design airfoil sections and 
wings for certain given pressure and lift distributions. One notable 
limitation is the phenomenon of stalling, which essentially depends 
on the velocity distribution in the immediate neighborhood of 
the solid surface, in other words, on the behavior of the boundary 
layer. 

However, when the velocity of an aircraft reaches about one- 
half of the velocity of sound, discrepancies appear due to the 
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neglection of the compressibility of the air, and it becomes im- 
perative to extend the theory to elastic fluids. 

Let us consider two-dimensional flow. In this case the motion 
of the fluid is completely defined by the stream function and the 
velocity potential function or, in other words, by the streamlines 
and the lines of equal potential. These two families of curves are 
orthogonal, and in the case of an incompressible fluid they con- 
stitute a pattern of geometrically similar square-surface elements. 
The mathematical expression of this condition can be written in 
the form: 


d(t> d\l/ d(t> d\// 

dx dy dx ( 4 *^) 

where (i>{x^y) is the potential function and ^p{x^y) is the stream 
function. The corresponding equations for an elastic fluid have 
the form 


p d(j> d\p p d<l) d\p 

Po dx dy po dy dx ( 4 - 2 ) 

where p is the variable density of the fluid po is its density at rest. 
The density is a function of the absolute value of the velocity; it is 
given by the relation 


/ 


dp dp 
dp p 



= const. 


( 4 - 3 ) 


is equal to the square of the velocity of sound and is a given 
dp 

function of the density p. 

The solution of the system of equations (4.2) and (4.3) is in- 
finitely more complicated than the solution of the system (4.1). 
The latter leads immediately to Laplace’s equation both for </> and 
and we can use the powerful and elegant methods of conformal 
transformation. In the case of (4.2) and (4.3) the differential equa- 
tion for 0 has the form 


/ U^\d‘^<j) 2UV d^(t> / _ / \ 

\ a^/dx^ dxdy \ a^/dy^ ^ ( 4 * 4 ) 


where u = v 


— and (u^ + v^). The quan- 

dx dy 2 


tity is a constant, equal to the square of the velocity of sound 
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in the gas at rest. The complicated nonlinear equation (4.4) can- 
not be expected to allow elegant and easy methods of solution. 

The equation for becomes linear if we consider u and i.e., 
the velocity components, as independent variables. In other words 
we draw equipotential lines = const, in the “hodograph plane.” 
The equation for as function of the velocity components result- 
ing from this transformation is a linear one, and it is possible to 
find particular solutions for it. For example, if we introduce polar 
co-ordinates in the hodograph plane, putting u ^ q co% 0 ^ v ^ q 
sin 0, where q is the absolute magnitude of the velocity and e is 
the inclination between the velocity vector and the x=axis, the 
differential equation for has particular solutions of the form 
0 = f{q)e^^^ where the function f{q) can be expressed by hyper- 
geometrical functions and can be calculated with sufficient 
accuracy by means of convergent power series. However, to obtain 
the corresponding flow pattern and velocity distribution we have 
to transfer the streamlines into the physical y plane. Then we 
sometimes shall find that solutions which look beautiful in the 
hodograph plane give, for example, overlapping flows in the 
physical plane. In other words our method is an inverse one, and 
an assumed relationship between the stream function and the 
velocity components does not necessarily give a flow with physi- 
cally possible boundaries. However, in a few cases the method has 
yielded useful results. A further promising application of the hodo- 
graph method to approximate solutions of flow problems is dis- 
cussed later in this section. 

Turning our attention to approximate methods of solution, the 
first idea that occurs is to apply the process of iteration. We might 
start from equation (4.1) putting p= constant, say in (4.2). 
Then we calculate p{x^y) from equation (4.3) and substitute this 
function in equation (4.2). This general idea was applied in differ- 
ent forms by Lord Rayleigh, O. Janzen, S. G. Hooker, L. Poggi, 
C. Kaplan, and I. Imai analytically, and by G. 1 . Taylor experi- 
mentally. The last-mentioned author used an electric computing 
device for the solution of (4.2). This device consists essentially of 
a layer of an electrolyte of variable thickness, the variable thickness 
being determined by the function p(x, j). The analytical iteration 
process employed by the other authors requires relatively much 
labor and it is difficult to have an estimate of the degree of approxi- 
mation achieved by the consecutive steps. 

A rather simple and ingenious method of approximation is 
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based on the assumption that the flow of a fluid around an airfoil 
or a streamlined body can be considered as a small perturbation 
of a parallel flow. If we denote the velocity and density of the un- 
disturbed flow by U and p®, and assume that u and p are only 
slightly different from U and p® respectively and v is small in com- 
parison to we can write i/jc + \l/ — poUy + and 

neglect terms containing squares and products of <^'5 p' and 
their derivatives. 

In this case the equations (4.2) become 

^ ^ 4. ^ 

Po Bx po By , po By 
and equation (4.3) takes the form: 


B\p^ 

Bx 


(4.6) 


a 



(4-7) 


Substituting (4.7) in the first equation (4.6) we obtain 

V a°^)dx ~ dj> 

d<l>' ^ 

By Bx 

Eliminating from these two equations we have 

/ t/2\dV' . 

V ° 


(4.8) 


(4-9) 


Equation (4.9) can be applied in the following way. Suppose 
we have solved the problem of the flow of an incompressible fluid 
around a cylindrical body of a certain cross section. Then the same 
velocity potential gives the flow of a compressible fluid around a 
cylindrical body whose cross section is affine to the cross section 
of the first body and can be derived from it by expansion in the 

f/2 

flow direction in the ratio i : \/i — where == — and U is 

the velocity of the undisturbed fluid. This consideration was sug- 
gested by H. Glauert and L. Prandtl. 

In the case of a streamlined body or an airfoil section the de- 
signer is mainly interested in finding the magnitude and the loca- 
tion of the highest velocity. It is assumed in general that if the 
velocity at any place reaches the velocity of sound (more exactly 
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the local velocity of sound), shock waves appear, the lift of the 
airfoil decreases, and its drag increases very considerably. If the 
thickness of the airfoil section and the camber of its center line are 
small in comparison to the chord length, and the section is con- 
tracted in the ratio i in the length direction, the maxi- 

mum value of the suction which determines the maximum velocity 

increases in the ratio — Hence, if the maximum suction 

Vi- 

f/2 

produced in an incompressible fluid is equal to Ppo — , we obtain 

2 

3 t/2 

in the case of a compressible fluid a suction equal to — ~ — • 

V I — 2 

Denoting the velocity increase by w, we have approximately the 

1 • rr ^ 

relation p^Uu = -r- — - ™po — or 

^ Vi - 2 


g U 
2V1 - 


(4.10) 


The velocity will reach the velocity of sound, U u = a% or 

— = I — Af . Hence, the critical Mach’s number, i.e., the ratio — 

for which the local velocity reaches the velocity of sound is given 
by the relation 


I - Af= - 


^ M 


1 




(4.1 i) 


or 


(i - + Af) = 


m 

2 


(4.12) 


For very small value of i — Af is small and we have 

(4-13) 

The comparison with the experimental evidence shows that equa- 
tion (4.12) underestimates the effect of compressibility; the equa- 
tion gives too large values for the critical Mach’s number. C. 
Kaplan suggested using a more exact relation between the velocity 
increase and suction. Then, he believes, the discrepancy is elimi- 
nated; however, it appears from the experiments that this correc- 
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tion is not sufficient. Moreover, from the theoretical point of view 
it is rather unsafe to mix up in the same consideration a linear 
approximation and relations which take higher terms into account. 

The present writer believes that closer and more consistent ap- 
proximation can be obtained by an adequate adaptation of a 
method used by A. Chaplygin, B. Demtchenko, A. Busemann and 
others for the flow of compressible fluids at low velocities. This 
adaptation has been suggested by the writer and worked out in 
an elegant way by H. S. Tsien. The method makes use of the trans- 
formation discussed above as the hodograph method. We consider 
the potential function and the stream function as functions of the 
velocity components. For sake of convenience we use the absolute 
magnitude, q and the inclination 6 of the velocity vector (i.e., the 
angle between the direction of the undisturbed flow and the 
velocity vector) as polar co-ordinates in the hodograph plane. Then 
the equations for and ^ in the case of an incompressible fluid 
have the form 


d<i> dy}/ d<t> d\l/ 


In other words, the complex quantity </> + ixp is an analytical func- 
tion of the complex variable log q — id. Let us now approximate 

the isentropic curve p = const in the neighborhood of the 

point po by a linear approximation 





(4-15) 


Then it can be shown that the equations for a compressible fluid 
whose physical behavior corresponds to equation (4.15) appear 
without further approximation in the form: 

p® d<t) dp d(j) p® dp 

If we now introduce a new variable 



CO = 


(4-17) 
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which is a generalization of log the variable used above, we have 

/o\ 

d(a dB^ 66 do) ' 


Hence, + iyf/ will be an analytical function of w — id. In this case 
we are able to find solutions of the system of equation (4.18) using 
the method of conformal transformation as it is used in the case 
of incompressible fluids. 

This method can be employed for the determination of the pres- 
sure distribution in a compressible fluid, provided the pressure 
distribution in an incompressible fluid is known. Let us assume 
that <p{o)y o) and d) are the potential function and the stream 
function respectively for the flow of an incompressible fluid around 
a given cross section, where w is the logarithm of the magnitude 
and 6 the angle of inclination of the velocity. Now if the values of 
q corresponding to w according to equation (4.17) are taken as 
the values of the velocity, we obtain a flow pattern satisfying the 
equations (4.16) of the compressible fluid. The relation (4.17) 
enables us to establish a simple relation between the magnitudes 
of the suction produced in an incompressible and a compressible 

fluid; this relation depends on the Mach’s number M De- 
noting the maximum suction produced in an incompressible fluid 
by — we have for the suction produced in a compressible fluid 






Vi — + 




2(1 + Vi - M ^) 


Ap 


(4-19) 


It is seen that for very small values of Ap the result is identical 
with the result obtained by use of the Glauert-Prandtl method. 
However, for larger values of the suction the effect of Mach’s 
number appears intensified by the influence of the second term 
in the denominator. As a matter of fact A^-~>oo for a value of M 
which is smaller than unity. 

If we return to the physical plane and draw the flow pattern 
mentioned above, we notice that the section undergoes a slight 
deformation. This transformation has been carried out by Tsien 
for several cases, and it seems that the influence of this slight dis- 
tortion of the cross section does not change the results to a signifi- 
cant extent. Moreover, the correction is likely to be balanced by 
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the deviation between the exact isentropic relation and the ap- 
proximate relation (4.15) used in this theory. It would be of great 
interest to check the predictions of the theory sketched above with 
a great number of measurements. Insofar as measurements are 
available on elliptic cylinders and airfoils, good agreement has 
been found. ^ ^ 

Flow problems in steam and gas turbines offer an as yet unex- 
plored field for applications of the theory of flow of compressible 
fluids. Many design problems related to such machinery are con- 
cerned with the flow of a compressible fluid through a lattice. 
Hence they require an extension of the methods discussed above 
and applied heretofore to single airfoils to an infinite system of 
identical airfoils. 

5. APPLICATIONS OF THE THEORY OF COMPRES- 
SIBLE FLUIDS TO EXTERIOR BALLISTICS 

The problem of the drag of projectiles is one of the oldest prob- 
lems of fluid mechanics and has attracted the attention of a great 
number of mathematicians and physicists, beginning with Isaac 
Newton. Newton’s theory of air resistance was later abandoned. 
As a matter of fact, the basic discrepancy between experience and 
d’Alembert’s paradox, according to which a body moving with 
uniform velocity through a non-viscous fluid does not meet any 
resistance, to some extent blocked the way to a rational theory of 
drag. The various fantastic propositions for best shapes of pro- 
jectiles reproduced in Figs. 4 and 5 show the utmost lack of any 
theoretical guide to empirical speculation.^ To be sure the funda- 
mental distinction between drag at supersonic and subsonic speeds 
has been pointed out by quite a number of physicists and ballis- 
ticians dealing with the problem; namely, that the energy loss 
corresponding to the drag at supersonic speeds is dissipated in the 
waves accompanying the projectile, in particular in the head wave 
emanating from the nose of the bullet, whereas the energy loss at 
subsonic velocities is mainly dissipated in eddies produced in the 
rear of the body. However, even an approximate quantitative 
theory of the wave resistance of shells was first attempted only a 
few years ago, and an exact theory does not exist as yet. 

1 From Th. H. Green, Army Ordnanci, May-Jime 1932. 
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The calculation of the retardation of projectiles by air resistance 
requires the solution of the problem of axially symmetric flow of 



Fio. 4. Small arms projectile designs prior to the American Civil War. 

I. Prussian; 2. Swiss; 3. Russian; 4. Prussian; 5. Arrow Ball; 6. Spanish; 7. Hunt- 
ing; 8. Russian; 9. Colts; 10. Russian; ii. Mangeot; 12. Richards; 13. Peters; 
14. Beckwith; 15. Russian; 16. American; 17. English; 18. Sardinia; 19. Enfield; 
20. Pomme k Pin; 21. Mexican; 22. Whitworth; 23. Swiss; 24. Prussian; 25. Beck- 
with; 26. French; 27. Mangeot; 28. Tamisien. 


an clastic fluid around a body of revolution for various velocities. 
It is not difficult to write up the differential equation for such a 
flow, as we did in the case of two-dimensional flow (cf. Eq. (4.4)). 
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However, the equation for the case of axial symmetry is even more 
complicated than the equation (4.4). Hence, we will consider first 



Fig. 5. The form of some of the projectiles patented m Great Britain. 

I. Brandon; 2 Farini; 3. Grenfell; 4. Trebor; 5. Unge; 6 Reichwald; 7. Hin|on; 
8. Parry; 9. Shipright; 10. Cole; ii. Simpson, 12. Hadfield; 13. Edmonds; 14. Lin- 
dcnufer; 15. Johnson; 16. Fisher; 17. Parsons; 18. Palliser; 19. Mennens; 20. Lake; 
21. Hazeltine; 22. Palliser; 23. Longsdon; 24. Hope; 25. Farini; 26. Day; 27. Page. 

the approximate linearized equation which corresponds to equa- 
tion (4.9). Then we easily obtain the following differential equa- 





FLUID MECHANICS 


36 

tion for the potential as function of the cylindrical co-ordinates 
r and x: 


( lP\d^<t> , av . I H 


( 5-0 


If U the solution of this equation can be reduced to the solu- 
tion of Laplace’s equation, as was explained in the previous sec- 
tion. However, ii U>a^ as in most problems in exterior ballistics, 
the equation is of the type of the wave equation: 


I aV ^ a^<^ 4. ^ 
df dr^ r dr 


which controls the propagation of small disturbances in a two- 
dimensional medium. 

In order to visualize the physical meaning of equation (5.2), let 
us assume that in a two-dimensional medium initially at rest a 
small “radial burst” of the character of a source at r = o occurs 
at the time / == o. The effect of this source will propagate in the 
radial direction so that at the time t the portion of the medium for 
which r>€t will be undisturbed and the effect of the disturbance 
will be restricted to the portion rdcU Therefore, if we choose t 
and r as cylindrical co-ordinates, the effect of the disturbance 
applied at r = o at the time ^ = o will be restricted to the inner 
portion of a cone whose vertex is at r = / = o and whose half 
vertex angle is equal to tarr^ c. The surface of the cone contains 
the circles r = which form the wave front at various times. If 
we apply this conclusion to equation (5.1), we have to keep in 

mind that c is replaced by — . Therefore, if we assume a 

stationary disturbance at an arbitrary point on the x-axis, say 
f = o, X = Xo the effect of such a disturbance will be restricted to 
the inner portion of a cone whose vertex is at r = o, x = Xo and 
whose half vertex angle is equal to 






(5-3) 


The angle /3 is called Mach’s angle. The relation (5.3) can be 
easily visualized in the following way. Assume that the fluid is at 
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rest and the disturbance moves in the negative x direction with 
the velocity U. Then at any instant a spherical wave will be 
created by the disturbance. Evidently the common front of these 
spherical waves is a cone with the half vertex angle equal to 

sin~i since in the time t the disturbance, i.e., the vertex of the 

cone, is displaced by the length Ut and the radius of the spherical 
wave is equal to aH. 

The solution of the wave equation of the type (5.2) has been 
worked out by T. Levi-Civita and H. Lamb. Their solution has 
been used by the present writer in a joint paper with N. Moore 
for the computation of the flow around a projectile in the following 
manner: We assume continuously distributed sources along the 
AT-axis between the nose and the base of the projectile and deter- 
mine the distribution of sources in such a way that the velocity 
normal to the surface of the shell vanishes, at least with the ap- 
proximation used in this linear theory. It is quite remarkable that 
the values of the drag calculated from the pressure distribution 
resulting from this simplified theory is in fair accordance with the 
drag vs. velocity curve derived from firing experiments. To be 
sure the accordance is only good in the case of slender and sharp- 
nosed projectiles. In the case of shells with blunt noses no such 
accordance can be expected since in that case no justification 
exists for the assumption of small disturbances. 

This linearized theory has been used by the writer for the com- 
putation of the best shape of an ogive; recently C. Ferrari made 
interesting contributions to the same problem, which first was 
treated by Newton. H. S. Tsien has extended the theory to the 
flow around projectiles whose axis is inclined at a small angle to 
the flow direction. He used a distribution of doublets along the 
axis of the shell, and computed the coefficients of the side force 
and pitching moment. These quantities have important bearing 
on problems in exterior ballistics. 

If we give up the assumption that the flow can be approximated 
by small deviations from a uniform parallel flow, we obtain in 
general at least one discontinuity surface which has the character 
of a standing oblique shock wave. A standing oblique shock wave 
is a discontinuity surface which the fluid enters in a direction dif- 
ferent from the normal. In this case only the normal component 
of the velocity suffers an abrupt change on passing through the 
shock wave; the tangential component remains unaltered. There- 
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fore, the streamlines show a break at the discontinuity surface. 
We already mentioned such oblique shock waves in connection 
with the nozzle problem. In fact, if the pressure at the exit of a 
nozzle is slightly higher than the calculated end pressure, the gas 
may continue its isentropic expansion to the end section; then at 
the edges of the end section the higher value of the external pres- 
sure creates oblique shock waves which are reflected at the free 
surface of the jet, until the surplus energy is dissipated by impact 
losses and friction. 

The exact shape of the oblique shock wave produced at the 
sharp edge of an obstacle can be computed exactly in the two- 
dimensional case by the methods of characteristics. Such compu- 
tation was, for example, carried out by L. Crocco for a sharp- 
nosed half-body. As a matter of fact in the case of two-dimensional 
flow the solution of the exact equation of an elastic fluid is much 
simpler for supersonic than for subsonic flow due to the fact that 
in the supersonic case the equation has real characteristics. 

However, in the case of three-dimensional flow, even if axial 
symmetry can be assumed, the solution is very laborious and the 
only case in which it has been worked out in detail is that of a 
conical obstacle of infinite extension. This problem was first treated 
by A. Busemann, who noticed that the partial differential equa- 
tion can be reduced to an ordinary differential equation. G. I. 
Taylor and J. W. Maccoll calculated numerically the solutions 
of this equation and discussed the various flow patterns in great 
detail. It appears that one obtains conical standing shock waves 
above a certain Mach’s number, whose magnitude depends on 
the vertex angle of the cone. It also appears that if the half vertex 
angle exceeds 57.6° no conical shock wave can be found at any 
velocity. In the range in which conical shock waves exist, there 
are two solutions corresponding to two different values of the in- 
clination of the shock waves. However, it seems that only one of 
the two shock waves is stable. 

The observations of the flow around sharp-nosed shells confirm 
the conclusion of the theory that shock waves can be found at the 
tip only above a certain minimum value of Mach’s number. For 
smaller Mach’s number it appears from the observations that the 
shock wave is detached from the nose of the shell and lies ahead of 
the nose of the projectile. It seems that the flow between the shock 
wave and the shell is subsonic at least near the axis of symmetry. 
It has not been possible as yet to compute the shape of such a 
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detached wave surface and the velocity and pressure distribution 
at the shell surface resulting from such flow conditions. Ferrari 
applied the method of characteristics to integration of the exact 
diflFerential equation of the flow with axial symmetry around a 
shell for the case in which the flow is supersonic everywhere. This 
method requires cumbersome calculations, but leads in principle 
to the solution of the problem. For the mixed supersonic-subsonic 
case, which occurs at lower Mach’s number, no proper mathe- 
matical method has even been indicated. 

Optical observations show behind a projectile a system of shock 
and expansion waves and some eddy formations of periodic char- 
acter (Figs. 6 and 7) ; there are only qualitative descriptive analyses 
available for this complex flow pattern.^ 

It appears that although the simplified mathematical analysis 
is a valuable guide for the solution of problems of high-speed 
aerodynamics related to exterior ballistics, the difficulties of an 
exact mathematical analysis and the complex character of the 
physical phenomena make further experimental research impera- 
tive. It is the firm conviction of the present writer that as in aero- 
nautics proper the wind tunnel — ^in this case the supersonic wind 
tunnel — will establish itself as a most valuable tool of research. It 
is bound to supplement the results of the firing tests and contribute 
new suggestions for further improvements in design. 

* Figs. 6 and 7 are reproduced from Crane’s textbook on Ballistics. Fig. 6 shows 
oblique shockwaves starting from the nose Fig. 7 (projectile fired backwards) 
shows a shockwave detached from the shell. 
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The movement by a fluid of solid particles heavier than the fluid 
itself is an important phenomena in many engineering problems. 
The transportation of dust and sand by wind and air currents 
occurs in problems with which the meteorologist and the me- 
chanical engineer are concerned. The movement of sediment by 
flowing water has been and will be more so in the future, one of 
the major problems of the hydraulic engineer. No doubt there are 
also many chemical engineering problems where the transporta- 
tion of various suspensions in liquids is of significance. 

The mechanics of flowing fluids must of course be first under- 
stood and studied if we are to know more about the action of the 
fluid on foreign particles. The interaction between the moving 
fluid and the heavier foreign particles is quite complex and has led 
to the formulation of a great variety of empirical rules and for- 
mulas. Considering the transportation of sand and silt by flowing 
water, experimental observations indicate that it is possible to 
describe three different methods by which this movement takes 
place. In order to provide a general background for the specific 
problem with which this paper is concerned, a brief description 
of these various methods or phases of sediment transport will be 
given. 

As water flows over a layer of unconsolidated sand grains a suf- 
ficient force may be exerted on the individual grains so as to dis- 
place them and cause them to roll along the channel bed. Such 
movement is referred to as “bed-load” transportation. Depending 
on the shape of the particles and the water velocity, some of them 
may become raised off the bed due to various hydrodynamic forces 
set up as the water flows over a particle. Once lifted off the bed 
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the particle will travel with the water for a short distance and then 
drop down. On striking the bed it may cause another particle to 
be dislodged off the bed and raised a short distance into the main 
stream. This particular process is referred to as “saltation,” and 
though little is known about the process quantitatively, it appears 
to be always present when sediment transport takes place. It was 
mentioned by Gilbert (i) in the report of his investigations, and 
has been vividly described by Bagnold (2) in regard to the trans- 
port of sand in deserts by high velocity winds. 

Of course, as soon as a particle of sand is raised off the bed it 
immediately tends to fall back due to the action of gravity. How- 
ever, this immediate dropping back to the bed may not occur if 
the fall velocity is counterbalanced by an upward component 
equal or greater in magnitude. Such upward velocity components 
normal to the mean direction of flow occur in turbulent flow. If 
the particles are kept in the main body of the liquid due to the 
turbulent velocity components, then the material transported in 
this fashion is called “suspended-load.” In most rivers and canals 
in the midwestern and western states the major portion of the 
sediment is transported in this manner. 

Since material is kept in suspension due to the turbulence of the 
fluid, it is logical to suppose that any analyses or studies of sus- 
pended material transpoitation will depend on our understanding 
of the turbulence mechanism. Within recent years our knowledge 
of this mechanism has increased tremendously, and because of 
this fact we have been able to carry forward our analytical and 
experimental studies of suspended material transportation at a 
rapid pace. 

A. GENERAL THEORY OF DIFFUSION IN A TURBU- 
LENT FLUID 

The simplest concept of how settleable particles are kept in sus- 
pension is to visualize the mass of fluid as filled with numerous 
eddies of various sizes whirling about in an apparently unpredict- 
able and random fashion. These eddies are capable of transferring 
heat, mass, momentum, foreign particles, etc., from one region of 
the fluid to another. The assumption has been made that the dif- 
fusion process due to the turbulence eddies is similar in general 
to molecular diffusion; the individual eddies performing the func- 
tion of the molecules. Experience seems to indicate that there is a 
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close similarity between the two phenomena if the scale of the 
turbulent diffusion process is large compared to the size of the 
eddies. 

The simplest case of diffusion which it is possible to study ex- 
perimentally is that of the spread from a line or point source of 
heat, or of matter of the same specific gravity as the turbulent 
water. The general equation of diffusion in a single direction for 
a steady state is: 


n = —D-j-f 
ay 


(0 


where n is the amount of matter transferred across a unit area in 
unit time, N is the concentration of the matter per unit volume 
and Z) is a “diffusion coefficient.’’ For the case of molecular dif- 
fusion it has been shown that: 


D = 1/2* (2) 

where ^ is the mean square displacement of the individual par- 
ticles (such as Brownian particles) after a time t. 

If particles having the same specific gravity as the fluid are in- 
jected into a turbulent stream of the fluid moving at a uniform 
velocity, their spread can be definitely calculated if the laws of 
molecular diffusion apply. If the concentration at the source is 
M then the concentration at a point a distance x downstream and 
a distance jv normal to the horizontal line passing through the 
point of injection is given by: 




2T 


(3) 


where ay is or the root mean square of the transverse travel 
of the particles. The value of cry depends on the downstream dis- 
tance, X. However, for molecular diffusion, as indicated in Eq. 2, 
f is equal to 2 i)/, where t can in this case be taken as x/U^ where 
U is the mean velocity of the fluid stream. In other words the 
value ofJj2 varies directly as the downstream distance. 

Observations of the spread of small immiscible droplets injected 
into a turbulent water stream indicate that the value of y^ does 
not vary directly as the downstream distance. In Fig. i is shown 
the actual variation. Note that the variation ofj^ for small values 
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of X is parabolic and gradually becomes linear for larger values 

of X, In other words ^ tends to attain a constant value. 
ax 



Fio. I . Experimental data on diffusion in a turbulent water stream.'^I IBM 

In order to explain analytically this phenomena Taylor (3) in- 
troduced the concept of a correlation coefficient, i?,, between the 
velocity of a fluid particle at any instant and its velocity at a time 
interval, t, later. With this concept Taylor develops the following 
relation for the variation of y with t: 



where is the mean square value of the velocity component 
normal to the direction of flow. Thus the diffusion process is in 
essence reduced to a study of the variation of Rt with time. It ap- 
pears plausible to suppose that in true turbulence /?, will approach 

zero as t increases, and when i?, is zero, the value of j' Rtdt be- 
comes constant. This means that the variation of yi with t (or 
with X for a uniform mean velocity, U) becomes linear. 
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A diffusion coefficient is defined for turbulence thus: 

= ( 5 ) 

By obtaimng motion pictures of the spread of small immiscible 
droplets injected at a point into turbulent water it is possible to 
calculate values ofj^ for various distances x, and thus determine 
the turbulent diffusion coefficient directly. Fig. 2 shows typical 
data of the variation of D from top to bottom in the center vertical 
section of an open water channel. The channel was 2.5 ft. wide 
and i.o ft. deep, and the mean velocity in the vertical section 
was .865 ft. per sec. 



Fio. 2. Variation of mass and momentum diffusion coefficients and of mean 
velocity in the center of a smooth rectangular open channel. 


The theory of momentum exchange due to the eddies in turbu- 
lent flow indicates that the shear in fully developed turbulence is: 


r 


dU 


(6) 


where e is the momentum exchange coefficient and p the unit 
density. Whether the coefficient D is equal to e has not been defi- 
nitely proven as yet, however, there is every reason to suppose 
that the two coefficients are at least proportional. In Fig. 2 are 
shown calculated values of e for the conditions for which D was 
determined experimentally. The computation of € was not very 
accurate since some assumptions had to be made regarding the 
value of the shear. The proximity of the numerical values of D 
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and € indicates that these two coefficients are probably identical. 
More data is being obtained to check this further under various 
conditions. 


B. THEORY OF SUSPENDED MATERIAL TRANSPORTA- 
TION IN AN OPEN WATER CHANNEL 


The analysis of the diffusion of sediment in turbulent water is 
complicated by the fact that it is constantly settling. If settleable 
material is placed in a turbulent liquid at a point, its diffusion 
will be quite similar to that of non-settling material, (Eq. 3), ex- 
cept that the point of maximum concentration instead of remain- 

ing at the same height is shifted down a distance equal to ^ > 


where c is the velocity of fall of the sediment in still water. 

The general differential equation controlling the phenomena 
of settleable material transportation in the main body of a turbu- 
lent fluid, considering just the two-dimensional case, is: 


aw aw bN _ UdN 

dt ~ df ^ dx’ 


( 6 ) 


where N is the mean concentration of sediment at any point {x^y) 
and at any time, t. The terms Dy and are the diffusion coeffi- 
cients in the designated directions. For the case when the sediment 
concentration in a channel does not vary on an average with time, 
and does not change downstream along the same elevation, Eq. 6 
reduces to: 


(PN 

dy ^ ^^dy^ ’ 


(7) 


or, 


cN^D,-^ ( 8 ) 

Eq. 8 states that the net diffusion of sediment across a unit area 
in the j-direction by the turbulence is equal to the amount settling 
due to gravity. This is the condition of equilibrium; that is, the 
suspended sediment concentration does not vary in the direction 
of flow. In other words there is no filling or scouring of the channel 
bottom. 
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Integration of Eq. 8 gives: 



(9) 


where is the concentration a distance a above the bottom. In 

^ * 

order to evaluate Dy it is necessary to relate it Jto known hydraulic 
factors of open channel flow. For turbulent uniform flow with a 
mean velocity gradient normal to the direction of flow the unit 

shear is represented as: r = where e is the coefRcient of mo- 

mentum transfer. The assumption is made that e is equal to Dy, 
Then knowing the shear and the velocity gradient the value of 
Dy can be calculated. Mathematically this can be done by ex- 
pressing the shear and the mean velocity as functions of the dis- 
tance j from the bottom. For the case of a wide channel, r = 

(i ~ where To is the bottom shear and is equal to ydS^ where 
7 is the specific weight of the fluid, S the channel slope and d its 
depth. The Prandtl-von K^rm^n law for velocity distribution in 
rough conduits, as transformed for open channels, is: 

U/U^ = I + ^^^(1 + log.j/t/) (10) 


= mean velocity in vertical section, 
k = von K^rm^n’s unusual constant, having a value between 
.38 and .40. 


The value of using Eq. 10, is \/ 7 j]>fky. On substitution in 
dy 

the equation for shear the value of the exchange coefficient is 
obtained: 


e = kdy/ To/p ( I - z)z, (11) 

where, z—y/d. 

Note that e is zero at the bottom and the surface, and is a maxi- 
mum at the mid-depth. Obviously, € does not become zero at the 
surface of most rivers, particularly those that show a definite 
“boiling.” Also, e is not zero at the zero level of suspended ma- 
terial, since at this level the suspended material merges with the 
bed-load and no definite solid surface exists. 



48 


FLUID MECHANICS 


Substituting Eq. 1 1 for Dy in Eq. 9 and performing the integra- 
tion we obtain: 



t = c/y/jo/p 

Though there is practically no experimental data for checking this 
relationship, Rouse (4) has presented some which appear to fol- 
low this relation quite closely. From a practical engineering stand- 
point Eq. 12 is not very convenient to use since the concentration 
is calculated as zero at the surface and infinity at the bottom. 
River measurements do not indicate that the suspended sediment 
concentration approaches zero at the surface. A more practically 
useful relation, and one that seems to approximate actual condi- 
tions quite well, is obtained if for the value of the diffusion co- 
efficient, Dy, in Eq. 9, the average value of the coefficient in the 
vertical section is calculated. This average value can be calculated 
from Eq. 1 1 : 

eav.. = rfVi^/15 (13) 

(The factor 15 is obtained if k is assumed to be 0.40). 

The equation for distribution of suspended material then be- 
comes: 

(14) 

This relation plots as a straight line on semi-logarithmic paper, 
and the concentration at any point in the vertical section can be 
determined by simply plotting the known concentration, JV„, and 
drawing a straight line through this point at a slope as determined 
from — 1 5^. In Fig. 3 are some actual measurements of suspended 
material concentrations in the Mississippi River at Muscatine, 
Iowa, as obtained with the assistance of the Rock Island U. S. 
Army Engineer Office. The straight lines drawn through the 
points indicating the concentration of sediment in a certain size 
range have a slope as calculated from Eq. 14, where t is equal to 
cl-y/T^, To is taken as ydS, and c the average velocity of fall of 
the sediment in any particular size range. Other field data for 
the Missouri River and smaller rivers also seem to check the vari- 
ation of suspended material as given by Eq. 14 very well. 

In order to obtain the total amount of suspended sediment 
transported in unit time per unit width of a section of channel it 
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is necessary to multiply the concentration per unit volume at any 
point by the mean velocity at that point. This involves combining 
Eq. 14 with the velocity distribution equation (Eq. 10). The total 
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Fig. 3. Distribution of suspended material in a vertical section of the 
Mississippi River. 


amount transported per unit time per unit width of any size range, 
characterized by the fall velocity, r, is then equal to: 


SjV. = (15) 

The relative roughness ratio comes into the function, </>, 

since it can be shown that in Eq. lo the quantity 
= l•^o^y'g when using Manning’s formula for the mean 

velocity, and 0.40 for the value of k. The function, n/d^^^) 
which we will call P, is plotted in Fig. 4. 

For the data shown in Fig. 3 the total amount of sediment 
carried of each size range was computed by multiplying the concen- 
tration at each point by the measured velocity, and was also com- 
puted using the bottom concentration as shown in Fig. 3 and then 
applying Eq. 15. The actual sum total of all sizes transported was 


50 


FLUID MECHANICS 


found to be i .85 lbs. per second per unit width. The work involved 
in obtaining this figure was quite laborious since it involved a 
graphical integration for each size range. The value calculated by 



Fig. 4. Field data on relation between suspended material, bed composition, and 
hydraulic characteristics of the sediment and the river or canal. 


use of Eq. 15 was 1.89 lbs. per second, and this result was obtained 
with the use of the diagram of Fig. 4 and a few simple calculations. 

C. RELATION BETWEEN SUSPENDED MATERIAL AND 
BOTTOM COMPOSITION 

It should be recalled that in the above analyses suspended ma- 
terial transportation under equilibrium conditions was discussed. 
Of course, in actual practice this simplified condition may not 
always be approached, especially for the fine sizes of silt. The non- 
equilibrium condition is very complex and was briefly analyzed 
in a previous paper by the author (5). For equilibrium conditions 
the concentration of any sediment increases toward the bottom, 
and therefore it is necessarily logical to suppose that the volume 
concentration on the bottom is greater than in any other part of 
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the stream. If this is not so, then our theory of suspended material 
distribution is invalid. 

The diffusion process due to the turbulence near a boundary, 
such as the bottom, is probably not similar in all respects to that 
in the main body of the stream. There is, of course, a transfer by 
the turbulence of material from the bottom into the main stream, 
otherwise all the suspended material would eventually settle out. 
The rate at which material settles out depends on the concentra- 
tion of the suspended material near the bottom and its velocity 
of fall. Under equilibrium conditions this rate of settling must be 
equaled by the rate of transfer of material upwards from the 
bottom by the eddies. Thus for any size of material there must 
exist a relation between the concentration of the material in sus- 
pension near the bottom, its velocity of fall, the amount of this 
material in the bottom, and some parameter or parameters char- 
acterizing the turbulence. 

It is assumed that the intensity of the turbulence is the important 
item which determines the ability of the turbulence to place ma- 
terial in suspension. Since experiments indicate that the velocity 
fluctuations tend to be statistically distributed according to the 
normal error law, it is possible to characterize the intensity of the 
velocity fluctuations by the root-mean square value of the fluctua- 
tions; that is, by the term \/p> where v is the transverse (upward) 
component of the velocity. The shear in turbulent flow is 

r = puv -f However, when sediment is being placed in sus- 

pension the liquid is quite turbulent and therefore p-j- is small 

compared to the first term. If it can be assumed that the correla- 
tion in the bottom region between u and v is in general the same 
for various flow conditions, and that and maintain a definite 
relationship, then To a Thus the bottom intensity of the turbu- 
lence is assumed to be proportional to \/r^, the so-called 
“friction velocity.” The various assumptions in the analysis, of 
course, need experimental confirmation. Such data is being ob- 
tained, but no definite general results are as yet on hand. 

If the value of y/r^jp characterizes the ability of the turbulence 
to pick sediment off the bottom, then a functional relationship 
can be set up between the suspended sediment concentration at 
the bottom, jV'o? the bottom concentration of any size of material, 
JVJ,, the velocity of fall of the material, and y/ro/p- Such a rela- 
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tionship was first developed by Lane and the author (6) and it is 
as follows: 


Mo/M= <t>{c/V7J~p) ( 1 6 ) 

Note that c/\/toIp is the parameter which we have been designat- 
ing by t in previous analyses. If a functional relationship such as 
Eq. 1 6 exists, then a plotting of No/N^ against t should give a 
smooth curve. No laboratory data exists which permits such a 
plotting; however, there are some meager field data. Using the 
best field data available, a plotting was made as shown in Fig. 5, 
and though there is considerable scatter of the points, there does 
seem to be a strong tendency toward a grouping about a definite 
curve. Considering the inadequateness of the data and the widely 
varying conditions under which it was obtained by various in- 
dividuals, the plotting in Fig. 5 seems quite significant. in 
Fig. 5 is in p.p.m. and in per cent.) 



Fig. 5. Curves for obtaining values of function in the expression for total rate of 
suspended material transport. 

The practical usefulness of a functional relationship such as 
shown in Fig. 5 is quite important. Such a relationship, after it 
has been definitely confirmed in the laboratory, permits the pre- 
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diction of whether under certain flow conditions the suspended 
load will be increased by scour or whether it will partly settle. It 
is noted that no material appears to be in suspension if the param- 
eter t is greater than about unity. In other words, when tht fall 
velocity of sediment in still water approaches the ‘Triction ve- 
locity,” transportation of such material in suspension becomes 
negligible. 


D. CONCLUSIONS 

A study of the diffusion characteristics of the turbulence in a 
channel indicates that the so-called diffusion coefficient is almost 
zero at the bottom and the surface, and reaches a maximum at 
about the mid-depth. Theory indicates that this is true for in- 
finitely wide channels; however, experiments in channels of finite 
width seem to verify this fact. 

Experimental data check the theory of suspended sediment dis- 
tribution based on the diffusion coefficient concept. Though this 
coefficient varies throughout any vertical section, the use of an 
average constant value for the coefficient gives a much simpler 
equation for the distribution of the suspended material. Field data 
for various large and small rivers indicate that this simpler rela- 
tion is adequate for practical engineering use. 

For equilibrium conditions there must exist some functional 
relationship between the bottom concentration of suspended ma- 
terial, velocity of fall of the material, bed composition, and the 
turbulence. Such a relationship is developed, and various field 
data seem to verify the existence of such a function. 

The items relating to the turbulence in an open channel which 
merit further laboratory investigation are: (i) Relation between 
shear, bottom roughness, and intensity of turbulence in the bottom 
region. (2) Relative values of longitudinal and transverse velocity 
fluctuations, that is, the relation between and (3) The in- 
tensity of the turbulence and the diffusion characteristics of the 
turbulence and the diffusion characteristics of the turbulence 
near a free surface. 

In regard to suspended material the most important item merit- 
ing intensive laboratory study is the mechanism of the placement 
of sediment in suspension in the bottom region. Further analyses 
relating to the prediction of the absolute amount of suspended 
material under various conditions depends on an increased knowl- 
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edge of the diffusion characteristics of the turbulence in the bed 
region. A laboratory verification of the relationship given as 
Eq. 1 6 is imperative. 

Though the analyses and studies described pertain more directly 
to the transport of suspended material by water, the fundamental 
principles should apply equally well for transport of dust and 
other settleable particles by air, or for the transport of any other 
type of solids by any turbulent fluid. 
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Flow 

By 

THOMAS K. SHERWOOD, Sc.D.* 

Chemical engineering is probably more dependent on the science 
of fluid mechanics than any branch of engineering other than 
naval and aeronautical. The various important chemical engi- 
neering operations of heat transfer, distillation, drying, humidifi- 
cation, and gas absorption involve the interphase transfer of heat 
or material, or both. Both the chemical engineer and the indus- 
tries employing these various transfer operations will benefit 
directly by the elucidation of the fundamental mechanism of inter- 
phase transfer of heat and material. Chemical engineering also 
involves the transport and flow measurement of all types of fluids, 
and chemical engineers were among the first to popularize the 
use of the Reynolds number as a basis for the correlation of pipe 
friction data. 

The important operations involving mass transfer are four: 
transfer between a solid and a gas, between a solid and a liquid, 
between a gas and a liquid, and between two immiscible liquids. 
As an example we may consider the steady-state evaporation of 
water from a wet surface into a turbulent air stream. Near the 
liquid surface the mechanism of the transfer process approaches 
that of molecular diffusion, about which much is known. In the 
main turbulent air stream the process approaches that of eddy 
diffusion, about which relatively little is known. The over-all proc- 
ess can be evaluated only on the basis of knowledge of the quanti- 
tative contribution of molecular and eddy diffusion at various 
points in the gas stream, particularly in the region near the pb^se 
boundary. 

The laws of molecular diffusion in gases have been developed 
from the kinetic theory, and have been carefully checked experi- 

* Associate Professor of Chemical Engineering, Massachusetts Institute of Tech- 
nology. 
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mentally. For steady-state diffusion of one gas (water vapor) 
through a layer of thickness j of a second stagnant gas (air) we 
have the relation 




DP _ DP 
R Ty pBM ypBM 


(I) 


where is the rate of diffusion of water vapor, g. mols/(sec.) 
(sq.cm.); D is the molecular diffusivity of water vapor in air, 
cm^/sec; P is the total pressure; aC is the difference in concentra- 
tion of water vapor in air, g. mols/cu.cm., at the two boundaries 
of the stagnant gas layer; and p^u is the logarithmic mean of the 
partial pressures of air at the two boundaries of the layer. The 
product DP is independent of total pressure. 

In the main turbulent gas stream transfer takes place predomi- 
nantly by turbulent mixing, or eddy diffusion, which is much 
more rapid than molecular diffusion. By analogy to (i) we may 
write: 




ay 


(2) 


but this relation is of little help, in that it serves only to define an 
“eddy diffusivity.” E. Dryden (6) points out that Na should be 
proportional to dC/dy for cases of diffusion over distances large 
compared with defined by Taylor as 


/2 = 



( 3 ) 


where Ry is the correlation between the values of the speed U 
(normal to the direction of diffusion) at two points distant j apart. 
If the distances are small compared with the material diffuses 
at a uniform speed, and (2) does not apply. 

Towle (17) has reported data on the spread of carbon dioxide 
and of hydrogen from a point source in a turbulent air stream. 
Equation (2) was found to apply, and the values of E obtained 
were very nearly proportional to the Reynolds number for the 
flow in the duct used. Hydrogen and carbon dioxide, differing 
22-fold in molecular weight, diffused at the same rate, suggesting 
that E may be independent of the weight of the molecules or very 
fine dusts transported by eddy diffusion. 

Woertz (14) has recently completed a study of steady-state 
transfer of water vapor across the 5.3-cm. gap between the flat 
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parallel walls of a large vertical rectangular duct 5.3 X 6i cm, in 
cross section. Water evaporated continuously from a water film 
on one wall, and was absorbed by a film of strong calcium chlo- 
ride solution on the opposite wall. The test section was preceded 
by a similar wetted-wall section to ensure steady-state conditions 
and minimize entrance effects. Hypodermic needle Pitot tube and 
sampling traverses were made to obtain the j/elocity pattern and 
the concentration gradient, respectively. The latter was found to 
be linear over the main central section of the duct, and values of 
E were obtained from Eq. 2 by dividing the rate of water vapor 
transfer by the slope of this curve representing concentration tra- 
verse. Helium, air, and carbon dioxide were used as the main gas 
stream. Velocity traverses and hot-wire turbulence measurements 
are also reported. 
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Fig. I. Water-vapor concentration traverses for steady-state diffusion of water 
vapor across a turbulent gas stream between two fiat walls 5.3 cm. apart. 


Typical concentration traverses as obtained by Woertz are 
shown in Fig. i, reproduced from the report of Woertz’s results. 
The straight central sections, and the abrupt drops near the duct 
walls are particularly evident at the higher Reynolds numbers. 
Calculated values of E for air checked Towle’s results very well, 
and the data for the three gases were correlated by plotting Ep 
vs Re^ as shown in Figure 2. 
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Various writers have pointed out that E should be proportional 
to the product of a mixing length and the deviating velocity in 
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Fig. 2. Data of Woertz on eddy diffusion compared with eddy viscosity based on 
velocity traverses in the same duct. 

the direction of diffusion. If the proportionality constant be a, 
then 

E = av'l 

where v' is the root mean square deviating velocity in the j; direc- 
tion (direction of diffusion), and I is the Prandtl mixing length. 
Since the shear stress, r, is given by 
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, , JU 
T= puv = pvl-^ 


dU 

djy 


where € is the eddy viscosity, it follows that 


(5) 


E = av'l = - (6) 

P 

Woertz’s data, as shown in Fig. 2, indicate a to be roughly con- 
stant at about 1.6, based on values of € obtained from velocity 
traverses and friction measurements in the same duct as employed 
in the diffusion tests. 


DATA ON MASS TRANSFER BETWEEN PHASES 


A large amount of data are available on solution of solids, evap- 
oration of liquids, gas absorption, etc., which may be compared 
with any theoretical equations proposed for mass transfer. As an 
example of such data we may describe briefly the results of Gilli- 
land (7) on vaporization of several liquids into a turbulent air 
stream in a wetted-wall column. The liquid tested evaporated 
from a liquid film on the inner wall of a vertical tube 2.67 cm. i.d. 
and 1 17 cm. long. Nine different liquids were used, giving a four- 
fold variation of the molecular diffusivity, D. The Reynolds 
number for the air stream was varied from 1,840 to nearly 30,000 
and the total pressure was varied from 0.14 to 3.07 atm. The re- 
sults were interpreted in terms of Equation i, the calculated values 
of B being taken as representing the thickness of a stagnant air 
layer which would have the same diffusional resistance as was 
actually encountered in the evaporation process. The data from 
a large number of tests were well correlated by the dimensionless 
equation 


^ = 0.023 (i2«) 



( 7 ) 


The correlation obtained is indicated by Figure 3, where d is the 
pipe diameter.* It follows from (i) and (7) that the rate ofvapQr- 

* Gilliland later ® obtained an equally good correlation of the same data by 
means of the equation 

^ = [903+^]^ (7a) 

The two additive terms represent the relative resistances of core and film, respect- 
ively. 
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ization is proportional to Z)®-®®, which conclusion agrees qualita- 
tively with the concept of transfer involving both molecular 
diffusion, for which is proportional to the first power of D, 



Fig. 3. Correlation obtained by Gilliland for vaporization of nine different liquids 
into a turbulent air stream in a wetted-wall column. 


and eddy diffusion, for which is proportional to the zero power 
of D. It also follows that Na is inversely proportional to as 
called for by Equation i . This suggests that the rate of eddy diffu- 
sion is also inversely proportional to Pbm- the tests described 
pBM ’was nearly proportional to p, so this conclusion agrees with 
Equation 6 (the group p/p-D is independent of pressure, and e is 
independent of pressure at a fixed value of Re ) . 

Gilliland’s data were obtained in a standard tubular shape, with 
calming section, and cover a reasonably wide range of values of 
Re and p/pi). They may be taken as typical vaporization data, 
and will be used for comparison with the theoretical relations dis- 
cussed below. 
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THE REYNOLDS ANALOGY AND THE COLBURN ADAP- 
TATION OF THE PRANDTL-TAYLOR EQUATION 

The well-known Reynolds analogy between heat transfer and 
friction is based on an assumption which may be expressed as fol- 
lows: the head lost due to friction, divided by the momentum of 
the stream, is equal to the ratio of the actual heat transferred to 
the heat which would be transferred were the fluid to come to the 
temperature of the heat-transfer surface. This assumption, as 
applied to flow in round pipes, gives 

h = \fcpUo (8) 

where h is the heat transfer coefficient, g. cal/ (sec.) (sq.cm.) (°C.), 
c is the specific heat of the fluid, p is its density, £/„ its average 
velocity (cm.sec.), and/is the usual friction factor in the Fanning 
equation. The analogous assumption for mass transfer is: and head 
lost due to friction, divided by the momentum of the stream, is 
equal to the ratio of the actual material transferred to the material 
which would be transferred were the stream to come to equilib- 
rium with the other phase. This is combined with a material bal- 
ance and the Fanning equation to give (15) 

^=K^koRT\fUo (9) 

In order to compare with Equation 7, we may take — ^ as unity, 
and assume 

/ = 0.039 (i?e)~® whence (9) becomes 

o.oig^{Rey^(j^ (10) 

This coincides with (7) exactly when ^ = 1.34, but obviously 

makes no allowance for the r 61 e of molecular diffusion {D canCels 
out when (i) and (10) are combined). The analogy (9) fits the 
data on diffusion of certain vapors, just as the heat transfer 
analogy (8) fits data on air in pipes within about to per cent, but 
agrees very poorly with data on heating water. 

A modification of the Reynolds analogy proposed by Prandtl 
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(13) and by Taylor (16) introduced the concept of a laminar film 
adjoining the phase boundary. It was assumed that heat passed 
through the laminar layer by molecular conduction, and the Rey- 
nolds analogy was applied only to the turbulent core of moving 
fluid. The resulting equation for heat transfer is 





(II) 


where k is the thermal conductivity of the fluid, and r is the ratio 
of the fluid velocity at the boundary of the laminar film to the 
average velocity of the main stream. The thickness of the laminar 
film, and hence the value of r, is determined by the critical value 
of the “friction distance parameter” jH-. If this criticial value ofj+ 
be taken as 11.6, as given by Bakhmeteflf (i), then r is given by 
the simple relation 

r=8.2Vy (12) 


Although (i i) represents a definite improvement over the Rey- 
nolds analogy (8), it is well known to be a poor approximation 
to the data on heat transfer to liquids, especially oils. 

The corresponding analogy for mass transfer has been derived 
by Colburn (4) as 


kc = 



(13) 


The relation of this equation to (9) is obviously similar to the re- 
lation of the Prandtl-Taylor equation (i i) to the simple Reynolds 
analogy (8) for heat transfer. For low partial pressures of the dif- 
fusing gas psM P are nearly equal, and the only change is the 


substitution of the group for —• 

pD k 

A further modification of the Prandtl-Taylor equation has been 
suggested by Chilton and Colburn (3), who proposed the substi- 

Yj for the denominator of 

( n \2/3 

— ] for the de- 
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nominator of (13). The two forms of the equation are the same, 

no matter what the value of r, for values of ^ or ^ of 1 .0. 

k pD 

The modified forms may be written 


h /cixY^ I . 

cUoAR) 2-^. 

(14) 

kcpBM { M I r 

~ C/> \pD) ~ J 

(15) 


Colburn (5) has presented correlations of heat transfer data in 
the form of graphs of ju vs. Re for several shapes, and Chilton and 
Colburn (3) have not only done the same for in mass transfer, 
but have shown that^i^ and j2> are approximately equal for the 
same shape at the same Re, This relation provides a valuable con- 
nection between mass transfer and heat transfer for a wide variety 
of conditions. 
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Fig. 4 . Comparison of data on skin friction and on vaporization of water from 
wet surfaces of various shapes. 


It is known that jjy approaches 1/2 /for flow past streamline 
shapes, but that may be only a small fraction of 1/2 /for flow 
past vertical cylinders or bluff objects. It is evident that the rela- 
tion between ^nd 1/2/ must be similar, but it is of interest to 
compare the two functions quantitatively. As a basis for compari- 
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son, it is convenient to use the recent data of Powell (12) who 
reports data on vaporization of water from flat surfaces, cyhnders, 
spheres, and disks placed in turbulent air streams. Other such 
data exist, but no other single investigator has reported data on 
a wide variety of shapes. The comparison is shown in Figure 4, 
which includes lines representing the vaporization data of Powell, 
together with the line FF representing Gilliland’s data on vapori- 
zation of water in a wetted-wall tube. Powell’s results are available 
only in the form of graphs, but the data appear to be well corre- 
lated by the co-ordinates of Figure 4, except in the case of spheres, 
where a slight spread of the data suggest an additional effect of 
diameter. The Reynolds number was based on the nominal di- 
ameter for the round shapes, and on the wetted length in the case 
of flat surfaces and cylinders placed parallel to the direction of 
air flow. The various curves fall fairly close together, the values 
of jj) for circular disks facing the wind being only about twice 
those for flow inside pipes. The data on cylinders placed parallel 
to the air flow represent an extension of the line for flat surfaces. 
The various shapes fall roughly in the order which might be an- 
ticipated, with the lowest curves representing the straight or 
streamline shapes. 

In comparing vaporization and friction data it is necessary to 
distinguish between the coefficient of total drag, and the skin 
friction coefficient/. The two are assumed to be equal for flow 

f 

inside pipes and over flat plates, and - for these shapes is repre- 
sented by the lines RR and QQ. The former is based on pressure- 
drop measurements in Gilliland’s apparatus, and checks approxi- 
mately with the standard curves for this shape. The line QQ, for 
flat surfaces represents the relation 

/= o.o 74 i ?^-°*20 (16) 

as given by Goldstein (9). The line PP represents the skin friction 
for a transverse cylinder, as given by Thom and quoted by Gold- 
stein (10). This last curve cannot be expected to be accurate since 
it is obtained by subtracting the normal pressure drag from the 
total drag, and f is small compared with C^), particularly at high 
values of Re. No data were located on skin friction for spheres or 
disks facing the air stream. 
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The results shown in Figure 4 are summarized briefly in Table I 
for a value of Re of 15,000. The ratio ofj^ to ^is very nearly unity, 


Table i 


Comparison of Data on Vaporization and^nction at Re = 15000 


Shape 

% 

Circular disks facing air stream (area based 
on one face only) . 

Spheres 

Transverse cylinders 

Flat surface parallel to stream 

Cylinders parallel to stream 

Inside pipes 


d or L 

• 

2Jd 

2Jd 

cm. 

JD 

/ 

Cd 

5 4-22 I 

0 0078 


0.014 

I 96-15-5 

0.0061 


0. 12 

0 16-37 2 

0 0060 

I.I7 

0 032 

1.8-24 3 

0 0055 

I .02 

1.02 

7 5-108 6 

(0 0055) 

I .02 

I .02 

2 67 

0 0043 

I.I4 

I.I4 


not only for the streamline shapes but for transverse cylinders. On 
the other hand, the ratio of j2> to half the coefficient of total drag 
is seen to be much less than unity for transverse cylinders, spheres, 
and the disks. 


THE MURPHREE ANALYSIS 

The Prandtl-Taylor concept of a laminar film at the phase 
boundary, with a wholly turbulent main stream, is now realized 
to be an over-simplification of the actual situation. Velocity tra- 
verses near the phase boundary point to a gradual change from 
laminar to wholly turbulent flow, and any successful theory of 
mass or heat transfer must be based on a reasonable guess or 
knowledge of the conditions in this boundary region. 

Although Murphree’s paper (ii) presents a theoretical relation 
between friction and heat transfer in pipes, it is not difficult to 
revise the analysis to obtain the corresponding relation for mass 
transfer. Murphree assumes the eddy viscosity to be constant in 
the main turbulent core, but in the region near the wall is assumed 
to be proportional to the cube of the distance j from the wall. If 
this region or “film” in which e varies is assumed to have a thick- 
ness a-b, where a is the pipe radius and b is the distance from the 
center to the film boundary, then in the film 
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where is the constant eddy viscosity in the main stream. The 
shear stress is given by 

(18) 

and it follows that in the turbulent core the velocity distribution 
it is given by integration as 

( 19 ) 

By combining (17) and (18) and integrating, an equation is ob- 
tained for the velocity distribution in the film. The result is com- 
bined with (ig) to obtain the velocity across the whole diameter, 
and it is found that the calculated velocity traverse agrees well 
with certain data of Stanton. Murphree then introduces Lees’ 
equation for/in round tubes and Stanton’s data for VJU„m brass 
tubes, and so obtains a general relation for and film thickness in 
terms of Re. The results are given in the form of a table of values 

<l> and -) reproduced here as Table II. The ratio - is the ratio of 
a a 

Table II 


Variables Employed by Murphree 


Re 

X 

5 

b 

a 


3000 

II . 16 

0 7203 

7.2 

5000. 

964 

0.8108 

II. 7 

10000 

7 94 

0 8835 

20 

25000 . 

6.24 

0.9318 

42 

50000 

5 28 

0-9541 

72 

100000. . 

4-54 

0.9700 

121 

300000 . 

3 66 

0 9837 

300 

500000. . 

3 36 

0 9^75 

470 

the radius out to the film to the 

pipe radius, and 4, is 

the ratio 


t'/fi. It is evident that the film thickness is determined by the 
intersection of two velocity distribution curves, that for the core 
being the parabolic velocity deficiency from (19) based on a 
constant e', and that for the film being based on an assumed 
variation of e with distance from the wall. 

In order to obtain the heat transfer coefficient, Murphree as- 
sumes the eddy conductivity to be ce, and integrates the conduction 
equation for heat transfer through the film. The coefficient h is 
defined as the heat flow through the film divided by the corre- 
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spending temperature-drop. For mass transfer we may employ 
Equation (6) and carry out the same integration to obtain 


ked 


where 


(■ - 


pu\ - 1 , 


Vs tan 




Vs 


(20) 


(-) 


and ^ represents the dimensionless group 

pij 

By employing Table II in conjunction with (20) it is possible to 
obtain a relation between the three dimensionless groups — , 

Kcd 

and Re, Murphree shows this derivation for heat transfer to check 
remarkably well with various data on heat transfer to air and 
water, and to check reasonably well with data on heating oils. 
Woertz and the author (14) showed the partial pressure gradient 
at Re = 5700 given on Figure i to agree fairly well with the curve 
calculated by the Murphree procedure, and erred in stating that 
the calculated film thicknesses became negative at high values of 
Re, Further comparison of this theory with mass transfer data 
will be presented after discussing K^rm^n’s recent contribution. 


THE KARMAN THEORY 


Von K^rm^n (18) has proposed a modification of the Prandtl- 
Taylor theory, introducing experimental data on velocity tra- 
verses as a basis for the calculation of transfer in the region of the 
phase boundary. As before, the Reynolds analogy is applied in 
the main turbulent stream. In some ways his procedure is similar 
to that of Murphree, in that the degree of turbulence is assumed 
to vary gradually through the film. 

K^rm^n’s derivation was for the case of heat transfer in pipes, 
but it will be of interest to follow the derivation for mass transfer, 
and to compare the result with vaporization data. Assuming the 
relation (6), we may write 



(22) 
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In the turbulent region, or core, both — and D may be taken as 

a. 

small compared to e. Assuming further that the ratio is 

constant, (22) and (18) may be combined and integrated to give 

U^-Uz _ a{Ci - C2) . , 

r pNa 

where the subscripts refer to any two positions in the wholly tur- 
bulent region. 

The shear stress t and the rate of mass transfer jV^ are taken 
as constant through a “film” of thickness B,, and the velocity at 
the film boundary, t/g, is obtained from (18) as 


Jo M + € 

(24) 

Similarly, from (22), 


^ ^ -^AP ny dy 

1 ^ 

(25) 


where Cg is the concentration at the wall, and Cg is the concen- 
tration at the film boundary. The variation of e through the film 
is based on the velocity gradient of Nikuradse (2), assuming the 

relation between and m+ to be 

general, and applicable to other fluids and pipes. Taking 

(26) 

then 



whence 

K^rm^n arbitrarily assumes the film boundary to correspond 
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to = 30, and approximates f (^+) between o and 30 by two 
simple functions : 

=/(/'^) =y^', fiy^) = I o <j+ < 5 

=f{y^) = 5 + i b<y^ < 3° 

5 y 

The comparison of these relations with Nikuradse’s data may be 



Fig. 5. Generalized velocity distribution, showing Nikuradse’s data and approxi- 
mation used by Kirm^n in region from y+ = 5 to y'^ = 30. (After Bakhmcteff.) 


visualized by reference to Figure 5. The combination of (29) with 
(28) gives the required variation of e, so that (24) and (25) may 
be integrated to give 


Ub=^ [5 + bln^)\/j 


(30) 


and 


= 54/^ t^+ (31) 


where, as before, = 


afi 


Now if (23) be assumed to apply in 
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the turbulent region with the average concentration replacing 

C2, and the average velocity £/<, replacing U^, then 

u.) ( 3 .) 

This is combined with (31) to give 

a{C.-Co) pUo , ./7r ..n , ^ 

— ' “^7 + sy - + 5 ' f ') - ^«6] (33) 

whence 



This result for round pipes may be seen to be of the same general 
form as the Reynolds analogy (9) and the Colburn equation (13). 

In order to compare this adaptation of the K^rm^n theory 
with experimental data on mass transfer in pipes, we may take 
the empirical relation (7) representing Gilliland’s tests. From (i) 
it is evident that if P and p^M are taken as equal, then 

d J^AdpBM kcd 
B ^ 'aCDP ^ W 


whence 


o^fUo ^ afU^ ^ fB 

2kc 2D 2 ^ 


(35) 


Substituting B from (7) and assuming/ = 0.039 R~^ 


2k 


:- o 4 T 


(36) 


In order to compare with this empirical form representing the 
vaporization data, the Kdrm^ form (34) is rewritten 


fUo 

•zkc 





(37) 


In making the comparison « is taken as 1.6, the approximate 
value found by Woertz. 
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Figure 6 shows the results graphically with plotted against 

Qikc 


pD 


for a Reynolds number of 10,000. The various theories indicate 



/OD V<=C/ 

Fig. 6. Comparison of the various analogies with data on vaporization in a wetted- 
wall column at Re - 10000. Curves: A, Colburn, Eq. 13; B, Chilton-Colburn, Eq. 
15; C, Murphrce, Eq. 38; D, Kiirmdn, Eq. 37; E, Gilliland empirical, Eq. 36; 
F, Reynolds analogy, Eq. 9. 
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relatively little variation with Re, so it is not surprising that 
Gilliland found no variation with Re (up to 30,000) of the ex- 
ponent 0.44 on the group The points represent the vaporiza- 
tion data for the nine liquids, and curve E represents Gilliland’s 
empirical equation (36). The Kirm^n equation (37), represented 
by curve D, is seen to be in remarkable agreement with the ex- 
perimental points. The maximum deviation is in the case of water, 
and is about 7 percent. The introduction of the factor «, not used 
by K^rm^n, is obviously of considerable importance in bringing 
about this agreement. The linear equation given in a previous 
footnote as Eq. 7a would be represented by a straight line on 
Figure 6 passing through the experimental points. This form 
would agree more closely with the K^rm^n curve D than does 
the alternative empirical form represented by curve R. 

Figure 6 also shows curves representing the various other theo- 
retical equations. The Reynolds analogy leads to the horizontal 


line F, independent of both Re and The Colburn equation 

pD 

(13) is represented by the straight line A, the position of which 
varies somewhat with Re due to changes in r. The Chilton- 
Colburn equation (15) is represented by line B, and is independent 
of Re, The curve representing the Murphree equation (20) is 
plotted by making the following rearrangement: 


fUo 

2ke 


= = 



6 


(38) 


The result varies appreciably with Re, and curve C would be low- 
ered about 15 per cent for Re = 50,000. 


EXTENSION OF THE KARMAN ANALOGY TO STEADY- 
STATE TRANSFER 


The remarkable agreement of the K^rmdn equation with the 
data, as evidenced by Figure 6, indicates that it is definitely su- 
perior to the other theoretical equations. In effect, Gilliland’s 
extensive data on vaporization of various liquids have been pre- 
dicted on the basis of Nikuradse’s velocity traverses in water, the 
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known relation between / and Re^ and a theory. It should be 
noted, however, that several questionable assumptions were in- 
troduced, and it may be asked if the errors due to wrong assump- 
tions could have counterbalanced each other. It is first assumed 

that IX and — are both small compared to c. This is certainly sub- 
stantiated by the results of both Woertz and Towle, and cannot 
be seriously questioned. At Re = 20,000 in air E and c are ap- 
proximately 6.3 sq.cm./sec. and 0.043 g./(sec.)(cm.) respec- 
tively, so IX is less than one half of one percent of c, and — is less 

a 

than 4 percent of e. It is next assumed that t/Na is constant in 
the turbulent core, whereas r is actually linear iny and the varia- 
tion of Na with j is different at different points along the pipe 
length. At the pipe inlet Na is large near the wall and small in the 
core; near the outlet of a long pipe Na is small at all values of^. 
In the analysis of the film r is assumed constant up to y+ = 30, 
although at low flow rates this value o{y+ may correspond to film 
thicknesses of half the radius, and the variation in r is actually 
50 percent. Similarly, Na is assumed constant in the film, although 
it must vary even more than r. In applying the Reynolds analogy 
(23) to the core, one limit was taken as the film boundary (y+ == 30) 
and the other as an indefinite point where U = and C = 
This would appear to involve the assumption that these equalities 
occur at the same point. Finally, the general applicability of the 
relation between and is assumed, and an approximate repre- 
sentation of this relation is based on Nikuradse’s data. This last 
must be accepted until better explorations have been made of the 
region near the phase boundary. 

Many of the assumptions introduced by K^rmdn may be elimi- 
nated if the same general reasoning is applied to the derivation 
of a theoretical equation for steady-state transfer across a turbu- 
lent stream flowing between parallel plates. This was the case 
studied experimentally by Woertz, and his data may be used to 
check the relation obtained. For the core assume: 

^ _l_ ^ -j_ g = /; = constant (39) 

a 

and write 
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where a represents half the distance between the plates. Integra- 
tion gives 


U^-U = 


zKa 


{a-yY 


Similarly, from (22) 
aK 


(C„ - Q = -JVCi(a -y) 


(41) 


(42) 


The film is taken as the region from jH- = o tojr*' = 30, = o to 

= 13.9, whence 

a Rey/jf 

where Re is defined as — — 


(43) 


From (42) and (43) 


(C„ Cb) - 0 

Combining (40) and (27) gives 


(44) 


K 


(f) (-3 




(45) 


so that the integration of (22) through the film takes the form 

dy 


c.-c..^r 

cc Jo 


pD 


+ .AT — /i 




dy* 


pD I 

+ rt/ 




This is integrated, as before, with the substitutions (29) to give 

I 


Cs-Cb = -In- 


a/JL 


I — 


5^ 

p « 


+ 


JVa /p^. (6o/> + 0.2 — g)(io/> + 0.2 + q) 
aq y To[ 6 op + 0.2 + ?)(iq/» + 0.2 — q) 


(47) 



MASS TRANSFER AND FRICTION 


75 


and q = 4/-I . /iZA 

y 25^i?ev7\ 'P ) 

The total concentration drop is 2 {Cg- CJ) and will be repre- 
sented by aC. It is obtained by combining (44) and (47) and re- 
arranging, to give 


UoAC 


Re I 
“ 2a ” 28.2 


ReVf 


A* / 169 
2fia\Re-\/j' 


■) 


. , (6o/> + 0.2 - q)(iop + 0.2 + g) 

^ aq\/f ( 6 oj> + 0.2 + q}(iOp + 0.2 — q) ' 

In this equation JT has been replaced by Re, since Woertz has 
shown both E and « to be nearly proportional to Re. 

It is possible to compare this relation with Woertz’s data, since 
he measured not only and the concentration gradient, but 
noted the temperatures and concentrations of the liquid films on 

the walls of the duct. Thus calculated from his results does 

not employ his concentration traverse data, and his values of /3 
and a are not based on the values of C at the walls. For water in 


air is 0.60, and taking a = 1.6, ^ is 0.96. From Figure 2 

the value of /3 is approximately 0.018 X io“®. As an approxima- 
tion, ^ may be taken as i.o for both CO2 and air, and (48) re- 
duces to 

6 --^ 

UqAC _ ^ I o Rey/f 

Na 3.2^^! — 28.2 vy” 169 
ReVf ^ ~ ReVf 

For helium ^ and fx are approximately 3.0 and 0.00020, respec- 
tively, and an equation similar in form to (49) is obtained. 

The comparison of these relations with Woertz’s data, based on 
values of f as measured by Woertz, is shown on Figure 7. The two 
solid lines are the calculated curves for carbon dioxide and air, 
and for helium, based on a = 1.6. The dotted curve is the calcu- 
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lated curve for air and carbon dioxide, with an approximate 
allowance for an increase of a at low values of Re, as indicated 
by Figure 2. The general agreement is good with the exception 
of a group of points for air at high Reynolds numbers. 
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Fio. 7. Comparison of data on steady-state vaporization across a rectangular 
channel with revised Kdrmin analogy (Eq. 48). 


In making this comparison it is important to visualize the nature 
of the experiment. Air passed vertically up through a rectangular 
passage 5.3 X 6i cm. in section. A film of water covered one wall 
and a film of strong calcium chloride solution covered the opposite 
wall. The gas was recirculated, with a wetted-wall calming section 
below the test section, so that there was very little net transfer of 
water to the gas stream in passing through the test section. aC 
was baised on values of C corresponding to saturation at the tem- 
perature and vapor pressure of the main body of each liquid layer. 
Since there is doubtless some temperature and concentration 
difference between the surface and the main body of the liquid 
an error is introduced in aC, varying both with liquid temperature 
and with or Re. These errors were not involved in Woertz’s 
determinations of E, which were based wholly on Mj, and the con- 
centration traverse. The same type of error was present in Gilli- 
land’s work, but was presumably not nearly so large because the 
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Reynolds numbers were much lower and because pure liquids 
were vaporized, giving no possible concentration gradients as in 
the case of absorption of water by strong solutions of calcium 
chloride. Considering the difficulties of the experimental technique 
the agreement with the theory is quite good. Insofar as K^rmin’s 
principal assumptions were eliminated in deriving (48), the result 
may be considered as lending further support' to the general con- 
cept on which K^rm^n’s analogy is based. 
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TABLE OF NOMENCLATURE 

a = pipe radius or half clearance between parallel plates, cm. 
b = distance from center line to film boundary, cm. 

B = equivalent film thickness, cm. 
c = specific heat of fluid, g. cal/(g.)(°C.) 

Cl = a constant 

C = concentration, g. mols/cu.cm. 
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Cj} 

CjD 

C. 

d 

D 

E 

f 

h 

jp 

Jh 

k 

K 

ka 


= concentration at film boundary, g. mols/cu.cm, 

= drag coefficient = total drag per unit surface area/~pt/^2 


== concentration at centerline of channel, g. mols/cu.cm. 

= concentration at phase boundary or pipe wall, g. mols/cu.cm. 
= pipe diameter, cm. 

= molecular diffusivity, sq.cm./sec. 

= eddy diffusivity, sq.cm. /sec. 

= friction factor = ^To/ pU^ 

= heat transfer coefficient, g. cal/(sec.) (sq.cm.) (®C.) 

= defined by first equality of Eq. 15. 

= defined by first equality of Eq. 14. 

= thermal conductivity, g. cal/(sec.) (sq.cm.) (®C. /cm.) 

== mass transfer coefficient, g. mols/ (sec.) (sq.cm.) (g. mol/cu.cm.) 


= mass transfer coefficient = 


ko 

RT 


= g. mols/(sec.) (sq.cm.) (atm.) 


K - fx+ or - — 
a 

I = Prandtl mixing length, cm. 

I2 = Taylor mixing length, as defined by Eq. 3. 

L = length of wetted surface, cm. 

Na = rate of mass transfer, g. mols/(sec.)(sq. cm.) 

Rey/f 

pji = partial pressure of diffusing gas, atm. 

pBM ^ logarithmic mean partial pressure of inert gas, atm. 

P = total pressure, atm. 



<1 

r 

R 

R, 

Re 


= ratio of fluid velocity at the boundary of the laminar film to 
average velocity 

= gas constant, (cu. cm.) (atm.) /(g. mol)(®K.) 

= correlation between two values of the speed U. 

= Reynolds number = — for pipes; = for flat plates and 


longitudinal cylinders; = 


4aUoP 


for flow between parallel 


plates. 

T = absolute temperature, ®K. 

tt' = root mean square deviating velocity in direction of main flow, 


cm. /sec. 



U = fluid velocity (time mean at a point), cm. /sec. 
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Uo 


y 


= fluid velocity at film boundary, cm./sec. 

= maximum or center-line velocity, cm./sec. 

= average fluid velocity, volumetric flow rate divided by channel 
cross section, cm./sec. 

= root mean square deviating velocity in direction of diffusion, 
cm./sec. 

= distance from phase boundary or pipewafl, cm. 


a = proportionality constant in Eq. 4, dimensionless 

^ ~ Re 

AC = driving force, or difference in concentration, g. mols/cu. cm. 

6 = eddy viscosity, g/ (sec.) (cm.) 

e = eddy viscosity in core, or main body of turbulent stream, 
g./(sec.)(cm.) 

H = fluid viscosity, g./(sec.)(cm.) 

p = fluid density, g./cu. cm. 

T = shear stress per unit area parallel to flow, dynes/sq. cm. 

Tf = shear stress at phase boundary or pipe wall, dynes/sq. cm. 

<l> = e/n 
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By 
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1. INTRODUCTION 

Modern statistical methodology is the result of the application of 
mathematical concepts and methods to the description and the 
interpretation of sets of observations of a repetitive character 
which arise in scientific experiments, in mass production, distribu- 
tion and consumption, in census work, and in many other fields. 
Statistical activity may be conveniently classified into two general 
categories. To one of the categories belong the routine collection, 
tabulation, description and presentation of masses of numerical 
data. This work has been greatly routinized within recent years 
by the use of machines; simple mathematical methods such as 
averaging, graphing, curve fitting, etc., are employed for purposes 
of condensing and presenting the data, and making simple com- 
parisons. In the other category belongs a body of statistical meth- 
odology built up over the last forty years which now constitutes 
an important part of certain scientific and technological proced- 
ures. The methodology in this category has been developed for 
making predictions and inferences from observations together 
with appropriate measurements of confidence. Here we find the 
mathematical methods somewhat more advanced with the theory 
of probability playing a very fundamental r61e. In the present 
paper we shall confine our discussion to the principles of this class 
of methodology and their scientific significance. 

* Associate Professor of Mathematics, Princeton University. 
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11 . NOTIONS AND PRINCIPLES OF PROBABILITY 

THEORY 

Approaches ranging from formal mathematical ones to empirical 
ones have been made to the theory of probability, none of them 
entirely without objections. From the point of view of applied 
mathematics, probability theory has been devised for describing 
certain empirical phenomena. Probability theory is also recog- 
nized as an important field in pure mathematics. Perhaps the most 
realistic approach is to examine informally a simple empirical 
situation and see how we should proceed to build a “model” to 
describe the aspects of the situation in which we shall be interested. 
It will be found that a model can be constructed which can be 
used to advantage under certain conditions for describing and 
making calculations in problems dealing with repetitive phe- 
nomena in a manner analogous to the way in which geometry can 
be used in describing and making calculations in problems of 
surveying, construction, etc. 

An Example. Suppose 30-watt electric light bulbs are made ac- 
cording to a certain set of specifications, and that we are interested 
in the length of life of the bulbs under continuous burning. If we 
let X denote the number of hours as determined under a given 
set of conditions that one of the bulbs will burn before expiring, 
then, of course, X will have different values for different bulbs. 
The sequence of bulbs coming from a production line is very 
similar in certain respects to a sequence of throws of a die. The 
variable X can be regarded as jumping from value to value as the 
bulbs emerge from the production line. In case of a die, we could 
let X denote the number of dots appearing on the uppermost face 
of a die. Then X would take on one of the values i, 2, 3, 4, 5, 6, 
at each throwing of the die, and of course would jump from value 
to value as the die is successively thrown. X is referred to as a 
random variable and it is used to measure whatever property or 
quality possessed by each member of a sequence of repeated 
events we may be interested in. The collection of events generated 
by such a sequence may be referred to as a population of events. 
In the case of the bulbs the value of X cannot be determined except 
through the destruction of the bulbs. In view of this destruction, 
a study of the population of bulbs emerging from the production 
line would ordinarily be made on a basis of a “random sample” 
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taken from the line — random as far as JY is concerned, although 
it may be taken by regularly choosing every 1 000th bulb. Now 
if we test a sample of n cases, a certain proportion of them, say 
Fn{x) will have values of X, i.e., life lengths, less than or equal to 
X. This statement will hold for all values of x. F^i^x) is known as 





the cumulative distribution function of X in the sample. The situ- 
ation may be represented graphically as shown in Figure i , where 
the n dots along the x axis represent the distribution of the n values 
of X in the sample, and the heavy step-like series of horizontal 
lines represent the graph oiy — F^{x), If there is only one observa- 
tion at each of the points shown on the x axis, then the jumps in 
F^{x) will all have the same value, namely i jn. It is a matter of 
experience that if n were increased the graph oiy = Fr^{x) would 
under certain conditions appear to be approaching a smooth 
curve something like the dotted curve in Figure i. Now what are 
these conditions? 

The idea of randomness , Before we could idealize and apply 
probability theory in connection with this example, we should 
have to obtain some assurance that no scheme of drawing a sub- 
sample from the actual sample of bulbs (the scheme of selecting 
a particular bulb depending only on information available before 
the value of X for the bulb is determined or used) would yield a 
cumulative distribution function, which, as n increases, appears 
to be “significantly” different from Ffx). Due to the difficulty of 
objectifying randomness and because of the practical necessity of 
having to deal with relatively small values of w, the problem of 
obtaining this assurance is, to a considerable extent, an empirical 
one relying on experience and sound judgment for its practical 
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solution. It has been found by Shewhart’ and others who have 
applied statistical methods in quality control work that the very 
process of securing randomness or ‘‘statistical control’^ in a sam- 
pling sequence of products such as bulbs yields extremely valuable 
information in discovering sources of variation in quality, and 
consequently in determining what steps should be taken in the 
manufkcturing process to achieve an increase in the level of 
quality in the units. For example, in the case of light bulbs, a 
study of the relationship between the burning life X and the in- 
formation available on the bulbs before X is made available will 
frequently indicate which factors in the manufacturing process 
are associated with relatively short-lived bulbs and which ones are 
associated with relatively long-lived bulbs. Steps can then be 
taken to eliminate those factors associated with short-lived bulbs. 
A new sequence of bulbs may then be drawn and tested, repeating 
the process until it is no longer economically feasible to remove 
causes associated with short-lived bulbs. 

The problem of establishing a state of statistical control in the 
sense used in quality control work is essentially a process of weed- 
ing out controllable causes of variation and securing greater 
homogeneity of product with respect to the variable under con- 
sideration. It should be pointed out, however, that there are situ- 
ations in which well-defined finite populations of objects exist, as 
contrasted to those generated by processes or operations, which 
are per se important subjects of study by sampling methods. For 
example, the population of eligible voters of a given state may be 
sampled in order to get an estimate of its preference toward a 
given gubernatorial candidate; the population of readers of a 
certain magazine may by sampled in order to obtain an estimate 
of its opinion of certain editorial matter; the population of wheat 
farms in the United States may be studied through samples. In 
such situations as these, the main problem is to leave the popula- 
tions as defined, however heterogeneous they may be, and devise 
methods of sampling which will yield sequences of drawings pos- 
sessing the characteristic of randomness which we have already 
informally described. 

Probability Functions. Once a state of statistical control or ran- 
domness has been established (largely by empirical methods) in 
the sampling sequence, one is in a position to idealize with prob- 

^ W. A. Shewhart, Economic Control of Quality of Manufactured Product. D. van 
Nostrand and Co. (1931). 
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ability theory, for a state of statistical control or randomness 
simply means that the cumulative distribution function for any 
scheme of subsampling remains essentially unchanged regardless 
of how the subsample is selected from the sample, provided the 
individual items are not selected on a basis of their respective X 
values. We would then assume- the existence of a cumulative prob- 
ability function F{x) which may be described, ^s far as its connection 
with the actual observations are concerned, as though it were a 
function “apparently” being approached by the cumulative dis- 
tribution function of the sample as a limit as n increases indefi- 
nitely. F{x) may be regarded for purposes of interpretation as the 
idealized proportion of objects in the population for which X x 
and is called the probability that X ^ x. Note that we are making 
no attempt to formalize in a purely mathematical fashion the con- 
nection between the observational function Ffx) and the assumed 
function F{x), The formal mathematical discipline begins only 
when the existence of F{x) is assumed. It is clear that the assumed 
cumulative probability function F{x) would be required to satisfy 
the following conditions: F( — 00) = o, 00 ) = I, F{xi) 
when Xi and X2 are any two values of X such that Xi<iX2. It will be 
observed that F{x2) — F{x^ is the probability that Xi<X^X2- In 
statistics, there are two especially important types of functions 
F{x) which may be described as follows: 

(a) That in which F{x) is a continuous function which can be 
differentiated with respect to x, whose graph is an ogive-shaped 
curve of the type shown as the dotted curve in Figure i. This case, 
known as the continuous case, arises in situations where the value of 
X is obtained from reading a continuous scale in terms of inches, 
pounds, etc. 

(b) That in which F{x) is a step function whose graph is a 
series of horizontal lines of the type shown in the jumps occurring 
at isolated values. The jump occurring at a particular value of x, 
say a;', is the probability that x takes on the value x\ This is the 
discrete case, and pertains to situations in which the value of X is 
obtained by counting, or more generally where X can take on 
only certain isolated values in the continuum of real numbers. 

In the continuous case and in view of the fact that F{x^ — F(^i) 


is the probability that Xi<iX^X2, it is evident that 


mdx 


=zf[x)dx, say, is (except for terms of order {dxY) the probability 
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that xciXc^x + dx. The function = — is known as the distribu- 
tion function or elementary probability law of X, The graph of 
dF . . 

y — IS, in most statistical situations, an inverted bell-shaped 

type of curve, which lies above the x axis. The area between the 
X axis and the curve and to the left of the ordinate x= x^h F{x^). 

In the discrete case, if we let p{x^ be the jump of F{x) at 
X = X,-, then y = p{x^ is the probability that X = a:,, and the 
graph oiy = p{x) would consist of ordinates erected at the isolated 
“possible” values of X and o otherwise. The mean value of the 
random variable X is defined by taking the product of each value 

of X and the probability corresponding to this value of X and sum- 

/ oo 

xf{x)dx or 
00 

for the continuous and discrete cases respectively. More 
generally, the mean value of any function of Xy say <p{X)y is given 

/ oo 

tp{x)f{x)dx or ^^ip{x^p{x^. In particular, the 

CO 

variance of X is defined by E\{X — £(A')) 2 ]. 0-2 ^ y^jy useful 

index and describes the amount of dispersion or spread in the dis- 
tribution of X, 

The extension of our discussion to situations in which two or 
more random variables are required to describe each element in 
the population is straightforward, but it will not be undertaken 


here because of lack of space. 

The problem of measuring on the basis of the evidence furnished 
by a sample of size n the credibility of the hypothesis that F^ix) 
is behaving as though it were converging to the assumed F(x) 
would have to be considered before we should know how much 


confidence could be placed in F{x) as the appropriate cumulative 
probability function for the population in the given situation. 
There are many sets of circumstances, such as are found in games 
of chance, genetics, etc., where one can calculate appropriate 
probability functions F{x) on an a priori basis by combinatorial 
methods with remarkable success, that is, one can by deductive 
methods arrive at functions F{x) which are approached very 
closely by the corresponding for practical values of n. How- 
ever, in problems of a more highly statistical type, such as the 
light bulb problem, the determination of the appropriate F{x) is 
not subject to such clear-cut a priori procedure. The usual pro- 
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cedure here is to assume F{x) to be of a certain functional form 
which involves parameters or constants which have to be estimated 
for the particular problem at hand. The problem of deducing the 
functional form often involves a considerable amount of mathe- 
matical ingenuity, and keen insight into the subject matter of the 
problem. ' 

The most commonly used cumulative distribution function is the 

rx 

normal or Gaussian form F{x) = / f{x)dx where f{x) is the well- 

known normal distribution function — 7^=- exp 

V27rO’ 

The parameters are a and <7, the mean and the standard devia- 
tion respectively of the distribution. By determining suitable 
values of a and <7, it has been found that the normal distribution 
function is satisfactory for many practical problems in statistics. 

III. SAMPLING THEORY 

At this point it should be made perfectly clear that probability 
theory as far as formal mathematics is concerned begins only when 
F{x) is given. Mathematical probability is concerned with the 
application of technical mathematical operations of various kinds 
to the idealized cumulative distribution function F{x ) . The theory 
of sampling and significance tests in statistics are built up in this 
manner. In the theory of sampling, we are usually interested in 
some average or index or ‘‘statistic” calculated from the values 
of X {xi, ^2? • • • 5 Xn) of a sample of n items or individuals. The value 
of fluctuates from sample to sample and in a set of jV repeated 
samples each of size n, has its own cumulative distribution 
function iv(^n) to which corresponds a cumulative probability 
function F(^„). The function F{Sn) is rigorously determined by 
technical mathematical methods from the underlying cumulative 
probability function F(:v) which is assumed to exist for the popula- 
tion of objects from which the sample is assumed to be drawn. 
As a matter of fact, F{Sn) in the continuous case is obtained by 
integrating the product f{x^f{x^ . . . f{x^ over all sets of values 
of the Jf’s for which ^ while in the discrete case F{s^ is 
obtained by summing the product • • . p{x^ over all 

sets of values of the xh for which ^ The field of theoretical 
sampling — the building of sampling distribution models of the type 
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F{Sn ) — has been developed at an extremely rapid rate by mathe- 
matical statisticians in England and in this country within the 
last quarter century. 

The problem of relating probability theory to empirical situa- 
tions is much more elusive and requires much more caution than 
that of building up the theory itself. The construction of the theory 
is a deductive process, while the task of applying the theory to 
statistical inference is a complicated mixture of trial and error, 
induction and deduction. The practical value of probability theory 
when its validity in a given statistical situation is established, lies 
in the fact that it furnishes an experimentally verifiable method of 
estimating the relative frequency with which it can be expected 
that the value of a random variable X will fall between any given 
pair of limits in repeated observations. The variable X may be a 
measurement on a single item or object; or more generally, it may 
be an average or some other index calculated from the individual 
Xs of a sample of n objects. 

IV. STATISTICAL ESTIMATION AND SIGNIFICANCE 

TESTS 

Sampling theory together with its application in statistical esti- 
mation and in the construction of significance tests constitutes a 
very significant contribution to scientific methodology. The sim- 
plest kind of situation in which sampling theory is used is that in 
which a population of objects exists and it is desired to estimate 
the value of some measurable characteristic of the population, for 
example the average of some variable or the proportion of the in- 
dividuals possessing a certain trait. Polls such as those conducted 
by Gallup are examples of simple sampling methods for estimating 
the proportion of eligible United States voters in favor of a given 
issue or candidate. The type of sampling carried out in such work 
is known as representative sampling, and the method consists in 
dividing up the population into smaller and more homogeneous 
groups and sampling at random from these sub-groups, the sizes 
of the sub-samples being proportional to the numbers of individu- 
als belonging to these sub-groups. Applications of representative 
sampling in this and more general forms are becoming more and 
more numerous in crop estimation, market research, national re- 
source surveys, etc., where it is impracticable or even impossible 
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to make an exhaustive survey of the population under considera- 
tion. 

A more general type of problem in estimation arises in situa- 
tions where there are grounds for assuming the cumulative distri- 
bution function F{x) associated with the underlying population, 
to be of a given functional form.4iepending on one or more un- 
specified parameters which have to be estimated from the sample. 
Thus if e is such a parameter, F{x) will be of a form F(x, 0) let us 
say. A fundamental assumption here, of course, is that there exists 
a value of e which will make F{x^ $) the actual cumulative proba- 
bility function for the population under consideration. We may 
refer to such a value of d as the “true” value. 

The Method of Maximum Likelihood, There are many ways of cal- 
culating estimates of B from a given sample. If we conceive of a 
repetition of samples of a given size w, then of course the estimate 
of B for the given method of estimation will fluctuate from sample 
to sample about the “true” value of 0 as a mean. This suggests 
that in constructing our mathematical model of sampling, we 
should try to find a method of estimation which, for a given sam- 
ple size, would yield estimates subject to a minimum amount of 
fluctuation from sample to sample. Such a method which is ap- 
plicable to large samples from populations having assumed under- 
lying distribution functions of a fairly wide class was discovered by 
R. A. Fisher^ in 1922. This method of estimation, now known as 
the Method of Maximum Likelihood, consists of determining the value 
of B which maximizes the probability of obtaining the set of values 
of X observed in the sample, that is, of maximizing B)f{x2^ B) 

• • • f{^n9 the continuous case and p{xi, B)p{x2, 0) . . . 

p{Xn, d) in the discrete case, /(a:, b) being used to denote — — 

in the continuous case and p{x', B) to denote the jump (proba- 
bility) at X = x' in the discrete case. The Method of Maximum Like- 
lihood provides not only the “best” estimate of in the minimum 
variance sense, but also the limiting form of the distribution law 
of this estimate in large samples from which one can determine 
the probability of obtaining a sample in which the estimate of B 
deviates from the “true” value by more than any given amount. 
Numerous improvements and refinements have been made by 

* R. A. Fisher, “On the Mathematical Foundations of Theoretical Statistics,” 
Phil, Trans, Roy, Soc, London, Series A, Vol. 282 , 1921. 
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Fisher, Neyman, and others on the general problem of estimation 
since Fisher’s original work was done. For example, the method 
of fiducial argument or interval estimation has been devised to 
provide in certain cases “observable” confidence limits of param- 
eters, even for small samples, which will include between them with 
a given preassigned probability the “true” value of the parameter 
to be estimated. 

The significance of estimation theory in the methodology of cer- 
tain types of scientific work is fairly evident as far as fitting distri- 
butions is concerned. The value of estimation theory goes much 
farther than this, however. 

The ^^Studenf^ Test of Significance in Experimental Work. A great 
deal of the work going on in experimental science is concerned 
with the variation of experimental factors and making comparisons 
of observations obtained under the different sets of experimental 
conditions. To return to the example of the 30-watt electric light 
bulbs, suppose it is desired to test the effect of a new filament de- 
sign on the life of the bulbs. A number of bulbs would be manu- 
factured in which the new filament design would be used — ^this 
being the only essential change introduced. Suppose a sample of 
new bulbs is selected and their individual burning lives deter- 
mined. Let X be the average burning life of the sample bulbs. 
Suppose the average burning life in bulbs with “regular” fila- 
ments as determined from a large number of tests is a. The ques- 
tion which now arises is the following: Is x “significantly” larger 
than a, or could the difference have “reasonably” arisen as a result 
of sampling fluctuations, under the assumption that the sample of 
new bulbs has been drawn from a normal population with the 
same average burning life as that for the bulbs with regular fila- 
ments? The answer to this question is furnished by sampling theory. 
Clearly, if the sample can be regarded as a not-too-improbable 
sample from a population which has the same average burning 
life as that of “regular” bulbs, then there are no grounds, on the 
basis of the sample of available information on the new filaments, 
for regarding the new filaments as providing longer burning life 
on the average. The usual statistical procedure for deciding whether 
or not i in a sample of size n for example is significantly greater 
than fl, assuming the sequence of new bulbs to be in a state of 
statistical control and that the burning lives of the bulbs are dis- 
tributed according to a normal distribution function, is to set up 
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a “significance level,” say .95, and calculate the two confidence 
limits 


X ± 



I ** __ 

where — xy, and t is 


such that 


r 

Values of t for various values of n and various significance levels 
have been tabulated and widely spread through statistical litera- 
ture. For example if w = 25, t = 2.02, and as n increases indefi- 
nitely, t decreases to 1.96. This test was, in its essentials, devised 
in 1908 by a British industrial statistician who wrote under the 
name of “Student.”^ It can be shown by sampling theory that 
under “ideal” conditions of repeated sampling about 95% of the 
samples will be such that the mean burning time of the popula- 
tion of new bulbs will be included between their respective pairs of 
confidence limits. If the smaller of the confidence limits exceeds 
the average burning life a, of the regular bulbs, we conclude that 
X is significantly larger than a at the .95 probability level, and that 
the filament is conducive to longer burning life on the average. 

It is true, of course, that in the event that “practically all” of 
the sample values of X are larger than a, we should be inclined to 
conclude on the basis of common sense that we were getting longer 
burning life by using the new filament. Although many experi- 
ments yield clear-cut results of this type in which common sense 
can be trusted to detect a genuine difference between experimental 
results, it is true that many classes of experiments do not fall into 
this category. Examples are those experiments in agriculture, in- 
dustry, market research, etc., where attempts are being made to 
discover ways and means of improving yield, quality, and per- 
formance, and where the improvements occur in very small intre- 
ments. Significance tests when properly applied will prevent, to within 
the significance level adopted, the error of interpreting a difference 
as genuine, when in fact it is essentially a difference due to sam- 
pling fluctuation. Of course there is a possibility that a significance 
* “Student,” “The Probable Error of a Mean,” Btometrika^ Vol. 6 (1908). 
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test will fail to reveal a genuine difference when it exists, but this 
source of error can be controlled to a large extent by obtaining 
assurance that the functional form of the distribution function 
F{x) and the associated significance test are appropriate and that 
a state of statistical control or randomness exists so that they can 
be applied. 

The Pearson Chi-Square Test of ^^Goodness^^ of Fit between Observa- 
tions and Hypothesis. We have discussed only one of many signifi- 
cance tests which have been proposed and used in practice for 
various types of problems. There are problems in which it is 
desired to test whether or not two or more samples (perhaps repre- 
senting as many variations of experimental factors) could reason- 
ably have originated from populations having identical distribu- 
tions without caring what the distributions are. Again, we may 
wish to test whether or not the sample cumulative distribution 
function Fn{x) is sufficiently close to a completely specified distri- 
bution F{x) deduced from a priori considerations to justify the 
statement that the difference between Fn{x) and F{x) can be “rea- 
sonably’’ accounted for as a sampling fluctuation. One of the most 
widely used criteria for making tests such as these is the Pearson 
Chi-Square Test. For example, we may wish to test after 600 
throws of a die which actually yield the frequencies 126, 96, 90, 
84, 96, 108 for I, 2 , 3, 4, 5, 6 respectively whether the die is be- 
having like a “true” die. That is, are the observed frequencies 
significantly different from the “expected” frequencies (each equal 
to 100) to reasonably conclude that the behavior of the die is 
consistent with that of a true die? The appropriate criterion for 
making this test is the Pearson^ Chi-Square criterion defined by 


^ -r 




100 


12. 1 


where /i = 126,^^ = 96, /a = 90, etc. It will be noticed that is 
an index for measuring how far the actual frequencies depart from 
the “expected” frequencies — the greater the collective discrepancy, 
the greater the value of and the smaller the collective discrep- 
ancy, the smaller the value of x^- Under “ideal” conditions of 
sampling, it can be shown that, in repeated samples of 600 throws, 
less than 5% of the samples are expected to yield values of x^ 

^Karl Pearson, “On the Criterion that a given Set of Deviations from the 
Probable in the Case of Correlated Variables is Such that Can Reasonably be 
Supposed to have Arisen from Random Sampling,” Phil. Mag., Vol. 50 (1900). 
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larger than that observed in the present example, i.e., 12.1. We 
should conclude on the basis of the 600 throws that the die is 
behaving in a biased manner (aces too frequent) although there is 
a small probability (less than .05) that our exceptionally large x* 
is being produced by a ‘‘true’’ die. 

Problems in which the ^ test ia one form or another is applic- 
able in testing the “goodness of fit” between observation and hy- 
pothesis arise in almost every branch of science. 

V. ANALYSIS OF VARIANCE AND ITS SIGNIFICANCE 
IN EXPERIMENTAL WORK 

There is a growing tendency in various fields of experimental 
science where statistical methods are applicable to make use of 
the principles of statistics in the actual layout and design of in- 
vestigations before any data is collected so that the results of the 
experiment can be tested and interpreted in a fairly simple man- 
ner, and in many cases so that the basis for inference will be 
broader. By using suitably arranged experiments and by making 
proper application of a technique known as the analysis of variance^ 
it has been found that one can in a single experiment obtain in- 
formation for inferring the effects of each of several experimental 
factors on the characteristic being studied in the experiment. To 
take a simple example described by Tippett, ^ suppose that in a 
cotton mill it is desired to study the extent to which weight per 
lea (yarn count) of a cotton yarn fluctuates from bobbin to bobbin 
and from doffing to doffing for a given spindle, holding other 
factors such as the raw material, operators, etc., as constant as is 
economically feasible. The variation due to these factors can be 
studied simultaneously by selecting a lea from each of the five 
bobbins in a given spindle and repeating this several times (thus 
obtaining yarn from several doffings) for the same five bobbins. 
If five doffings are used we then have 25 leas, five coming from 
each of five bobbins and five belonging to each of five doffings. 
The variable in which we are interested is weight per lea. If we 
let be the weight of the lea from the jth bobbin in the*zth 
doffing, we can arrange the data from this investigation in the 
following table: 

® L. H. C. Tippett, “Some Applications of Statistical Methods to the Study 
of Vaiiation of Quality of Cotton Yarn,” Supplement, Jour. Roy. Stat. Soc., Vol. 2 

(1934)- 
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Bobbin 



where the at the side are doffing means and those at the bottom 
are bobbin means, while x is the mean of all 25 cases. We can set 
up as indices of inter-bobbin and inter-doffing variation the sums 

of squares of residuals 55= 5 S(a:., — and 

1 1 

6 5 _ 

The total variation in the data is S = 'E — x)^, and the resid- 

Jxj tx, 

ual or “error” variation after inter-bobbin and inter-doffing 
variation have been eliminated is S -- Genuine 

variation in yarn count due to differences in bobbins will be re- 
flected by excessively large values of Sf, as compared with 5 ^; a 
similar statement holding for doffings. The question now arises: 
When can Sj, or be regarded as excessively large? This question 
can be answered by sampling theory under certain assumptions. 
The assumptions are that (a) each x^j can be regarded as being 
the sum of a general mean m, a constant due to the ith doffing 
and a constant Cj due to the jth bobbin, and an “error” residual 
(b) the are all distributed according to the same normal 
law having a = o, but unknown variance and (c) Sr, = Sc, = o. 

* j 

To test whether or not there is excessive variation from bobbin to 
bobbin we make the assumption that the c, are all o, and compare 
Sj, and i /4 and 1/16 can both be regarded as independent 
estimates of the same unknown quantity, namely the variance 
(T^ of the population to which the are regarded as belonging. 
Under the assumptions that the c, are all zero, the sampling dis- 
tribution of the ratio F = is known, from which, 

for a given probability level, it can be determined whether S'j is 
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too large as compared with to lend credibility to our assump- 
tion that the r, are all zero. If is found to be significantly large 
in relation to S^y we refute the assumption that the Cy are all zero, 
and conclude that there is a genuine inter-bobbin variation. A 
completely analogous procedure can be carried out for determin- 
ing whether or not the inter-dofl&ng variation can be considered 
genuine. It should be remarked that one can test the plausibility 
of the additive assumption (a) by duplicating the experiment. 
The significance of establishing sources of variation in this manner 
as a step in certain kinds of problems in science and technology 
is self-evident. 

It is possible to design the experiment so as to get information 
regarding the existence of further sources of variation of weight 
per lea by using more elaborate layouts. The analysis of variance 
is a comparatively recent development due largely to R.A. 
Fisher® and his followers, and it has rich possibilities in scientific 
investigation in many fields when properly applied. 

VI. SUMMARY 

There are two main categories of statistical activity: one deals 
with the routine collection and presentation of masses of statistical 
facts, and the other is concerned with a methodology founded on 
probability theory for making inferences and predictions from 
observations. The present paper deals with the latter category. 
Notions and principles of probability theory and their application 
to problems of sampling, statistical estimation, and significance 
tests are developed. The significance of sampling theory, statistical 
estimation, and significance tests in scientific and technological 
procedures is described and illustrated by examples. The value 
of laying out various types of experiments and investigations ac- 
cording to certain statistical principles is briefly discussed. 

® R. A. Fisher, Statistical Methods for Research Workers^ Oliver & Boyd, Edinburgh. 
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INTRODUCTION 

The potential contributions of statistics to the science of engineer- 
ing are an important national asset; an asset of interest to all of 
us because it makes possible the most efficient and effective use 
of natural resources and human effort to satisfy human wants; an 
asset, however, that has for years remained frozen and is only now 
beginning to be utilized; an asset that can be used to the full only 
when engineers and others learn to use it as they have learned, 
in the past, to use the product of the scientist. 

Much has appeared in the literature to indicate some of the 
contributions of statistics to date in the field of engineering and 
manufacturing. My object is not so much to review what has 
been done as to survey the potential contributions of statistics to 
the science of engineering. In doing this, I shall follow the old 
advice that the easiest way to reach the top is to go to the bottom 
of things, and I shall go to the bottom of the difference between 
engineering with, and without, statistics. 

Let us recall how the applied scientist has wrought so many 
wonders for you and me to enjoy. He has done much of this with 
a comparatively simple but extremely powerful tool, namely, 
scientific method based upon the concept of physical laws of 
nature that assume perfect or certain knowledge of a set of facts 
and then state exactly what will happen at any future time. This 
method consists of three essential steps: hypothesis, experinjent, 
and test of hypothesis. The fundamental difference between engi- 
neering with and without statistics boils down to the difference 
between the use of a scientific method based upon the concept of 
laws of nature that do not allow for chance or uncertainty and a 
scientific method based upon the concept of laws of probability 

* Research Statistician, Bell Telephone Laboratories. 
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as an attribute of nature. When viewed in this way, the potential 
contributions of statistics become quite simple indeed. Statistically 
scientific method is in fact a fundamental discipline that includes 
all of the customary scientific method based upon the concept of 
exact laws as a limiting case in which repetitive operations of any 
given kind always give identically the same results. 

All that I shall try to do here is to sketch in broad outline how 
the statistician has helped to refine scientific method by over- 
hauling and reworking each of its three fundamental steps and 
then to indicate briefly how this new tool can be used by the 
engineer. In place of hypotheses based on exact laws, the new 
scientific method introduces statistical hypotheses. The introduc- 
tion of statistical hypotheses makes it necessary to conduct the 
experiments in such a manner that the statistical hypotheses may 
be tested, and the new method provides means of testing the sta- 
tistical hypotheses with which it starts. 

Needless to say, scientists have long realized that their dis- 
covered “laws” do not always fit observed phenomena exactly. 
All that they have claimed is that the scientific method based upon 
the concept of exact laws has enabled them to make remarkable 
progress in understanding the world and to attain knowledge that 
could be used by engineers and applied scientists. Deviations from 
the assumed exact laws were simply dismissed as errors by the 
pure scientist and allowed for in factors of safety by the engineer. 
However, we are now coming to realize that many of the errors 
of the pure scientist and factors of safety of the engineer are more 
properly designated as factors of ignorance. We are beginning to 
see that we must refine our scientific methodology if we are to 
minimize these factors of ignorance in the scientific explanations 
of the world and if we are to make the most efficient and economic 
use of natural resources in the development of things to satisfy 
human wants. Statistically scientific method provides the scientist 
with an improved tool by which to extend his knowledge, and the 
engineer with a means by which to extend his useful service to 
mankind. 


BASIC ENGINEERING PROBLEM 

The engineer’s job is to devise and develop the operations that, 
if carried out, will produce things that people want.^ To do this, 

^ This is also the job of the applied scientist in the development of ways and 
means of making everything that we use. 
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he must be able to make things that have quality characteristics 
lying within previously specified tolerance ranges. Hence, a basic 
engineering problem is to devise an operation of using raw and 
fabricated materials that, if carried out, will give a thing wanted. 
The specified tolerance ranges for the quality characteristics of 
the thing wanted define a target for the engineer. He devises an 
operation and predicts that, if carried out, it will hit the target, 
but, since he does not have certain or perfect knowledge of facts 
and physical laws, he cannot be certain that a given operation 
will hit its target; in fact the best that he can hope to do is to know 
the probability of hitting the target. Here then is one fundamental 
way in which probability enters into everything that an engineer 
does. 

Furthermore, if the thing produced fails to meet tolerance re- 
quirements, the engineer is penalized in one way or another. For 
example, if the quality of any piecepart fails to meet its tolerance 
requirements, a loss is incurred through rejection or modification 
of the defective part; if the time-to-blow of a protective fuse fails 
to meet its tolerance range, loss of property and even loss of life 
may result; if the time-to-blow of a fuse in a shell fails to meet its 
tolerance range, the shell may burst prematurely and kill mem- 
bers of the gun crew and, in any case, the round of ammunition 
will fail to fulfill its function of destruction within the ranks of the 
enemy. This means that when the engineer undertakes to use 
probability theory it is essential that he thoroughly understand 
the conditions under which its use will lead to valid predictions. 

In what follows we introduce the term operation to include any 
experimental procedure designed to produce a previously specified 
result. In this sense, a production process is an operation, and a 
method of measuring is also an operation. Furthermore, an en- 
gineering operation or a production process may almost always 
be broken down into component operations. It should also be 
noted that even if only one thing of a kind is to be made, the oper- 
ation devised by the engineer for producing this one thing is pre- 
sumably capable of being repeated again and again so that the 
one thing to be produced may be thought of as but one of a class 
of an indefinitely large number of things that might be produced 
by repeating the operation again and again at will under the same 
essential conditions. In this way we may reduce the basic engineer- 
ing problem of devising an operation to hit previously specified 
tolerance ranges to one that can be treated statistically, in that 
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statistical theory treats of the properties of certain kinds of repeti- 
tive operations. 

To illustrate, let us consider a simple example in which only 
one quality characteristic of a thing is specified, let us say the 
length of a piecepart or the time-to-blow of a fuse. Let us symbolize 
this quality characteristic by X The operation for producing a 
thing of quality X within specified limits,^ if repeated again and 
again, would give rise to an indefinitely long sequence of values 
of quality X if taken in the order in which the things were pro- 
duced, and these may be represented symbolically as follows: 

Xx^ X^,^ . . . , Xiy • . . , X • • • , Xny . • • , • • • (l) 

Hence, every operation may be characterized not only by a word 
description but also by the characteristics of the potentially in- 
finite sequence (i) corresponding to an indefinitely large number 
of repetitions. In fact, we shall again and again make use of the 
characteristics of such a sequence in what follows. If an operation 
is developed to produce only one thing of a kind, then the engineer 
is interested only in the first term of this potentially infinite se- 
quence (i) but if the engineer is interested in developing an opera- 
tion or production process to turn out an indefinitely large number 
of things of the same kind, then he is interested in all terms in 
this sequence. 

In what follows, we shall try to see how statistics can help the 
engineer to solve his basic engineering problem of developing an 
operation that, if carried out, will produce an object with qualities 
that lie within previously specified tolerance limits. 

BASIC CONTRIBUTION OF CLASSICAL 
STATISTICAL THEORY 

Basic Statistical Hypothesis, As a background for viewing the con- 
tribution of statistics to the solution of the basic engineering prob- 
lem, we may state the fundamental hypothesis of applied classical 
statistical theory in the following way: 

Hypothesis /. Some repetitive operations exist that obey laws of 
probability. These are called random. The probability that such a 
random operation will give a previously specified events as for example^ 
the occurrence of a value of X within a previously specified tolerance 
range^ is a definite number associated with that event. 

* Or any operation of measurement of some objective quality characteristic. 
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If we know the law of probability or chance that controls a given 
operation, we may use the mathematical distribution theory of 
the statistician to describe how statistics of samples of size n given 
by successive repetitions of such an operation will be distributed. 
Likewise, if we know that an observed sample has been given by 
a random operation, the statistician has established valid rules of 
procedure for using the sample as a basis for estimating the param- 
eters in the law of probability underlying the random operation 
that gave the sample. In other words, the applied mathematical 
statistician has provided us with rules for making valid predictions 
if we know the law of probability and with efficient rules of dis- 
covering the functional form of such a law, including the values 
of the parameters, if we simply know that it exists. 

Of course, the work of the mathematician is purely formal; for 
example, what he calls the operation of drawing samples of size n 
at random consists essentially of acting upon some given mathe- 
matical law of chance or distribution function in accord with 
previously specified mathematical rules. Hence, if such a statistical 
hypothesis is to be of any value in engineering or applied science, 
it is necessary to know what is meant in an operationally verifiable 
manner by drawing at random. This necessitates our study of the 
second or experimental step in a statistically scientific method. 

Basic Statistical Experiment or Operation of Drawing at Random, Un- 
less an experimentalist knows what it means to draw a sample at 
random, he is not in a position to make use of statistical hypotheses 
because he cannot get the data with which to make valid tests of 
them and, having accepted a statistical hypothesis as valid, he 
does not know what kind of events he can predict with validity. 
Without such knowledge, he would be somewhat like a physicist 
who knew all about mathematical physics but did not know how 
to distinguish and measure the physical properties appearing in 
his equations. For example, an engineer may wish to select a 
random sample of 50 pieces of a new kind of product from the 
first 1,000 pieces produced. Very often he will ask under such con- 
ditions if he can take every twentieth piece as it is produced. Some- 
times the engineer will propose other schemes but, in general, it 
is found that none of those proposed can be used with much 
assurance of giving a random sample. Hence the starting point 
in the use of statistical theory is a clear understanding of what is 
to be taken as the meaning of a random operation. 

Let me begin with a description of two operations that I shall 
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choose to call random,® Let us assume that we have N physically 
similar chips on each of which is written a number and that these 
N chips are placed in a bowl. A blindfolded experimentalist thor- 
oughly mixes the chips in the bowl, draws a number, and has his 
assistant record the number. The chip is then returned to the bowl 
and the blindfolded experimentalist, after thorough mixing of the 
chips, again draws one and has the number recorded. This ran- 
dom operation of drawing a number can theoretically be repeated 
again and again under the same essential conditions so as to give 
an infinite sequence like (i). 

The other important operation for our present study is the fol- 
lowing one of randomizing a finite set of N numbers. In this case 
the blindfolded experimentalist follows the same procedure as 
described above except that he does not return a chip to the bowl 
after it has been drawn. The N numbers drawn in this manner 
may be written down by the assistant in the order drawn. The 
operation of putting the N chips in a bowl, thoroughly mixing 
them, drawing them one at a time, and writing the numbers down 
in the order observed, may be repeated again and again at will 
so that the result of the operation of drawing one sequence is but 
one of the infinite number of sequences that might be obtained 
by repeating the operation again and again.^ As an example. 
Fig. I a records one such random drawing of the 144 values of 
thickness of inlay given in Table i. 

Of course, you will note many elements in my description of 
the experiment that are not operationally definite. What, for ex- 
ample, are symmetrical chips? What is thorough mixing? How 
can the experimentalist maintain all other conditions essentially 
the same? However, this kind of indefiniteness in the operation of 
drawing at random is much the same as exists in defining any 
experimental procedure. 

The engineer without statistical training is likely to ask what 

* A more comprehensive discussion of the difference between the mathematical 
concept of random and the operationally verifiable meaning of random has recently 
been given elsewhere. 

Cf. W. A. Shewhart, Statistical Method from the Viewpoint of Quality Control, Grad- 
uate School of the Department of Agriculture, Washington, D. C. From the more 
technical viewpoint there presented, the two random operations here may be 
treated as one. For our present purpose they can best be considered separately. 

See also the discussion of what the mathematician defines as random in the paper 
by Prof. Wilks, also contributed to this symposium. 

* Of course there would only be N different possible orders. 
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TABLE I 


Thickness in Arbitrary Units of Inlay on 144 Relay Springs 


55 

13 

45 

27 

33 

59 

43 

24 

43 

52 

39 

62 

72 

31 

38 

36 

38 

45 

26 

30 

28 

42 

38 

50 

93 

37 

23 

18 

34 

41 

2 

28 

38 

53 

49 

70 

54 

38 

58 

43 

59 

47 

17 

71 

28 

53 

38 

60 

30 

43 

52 

10 

75 

38" 

28 

57 

23 

42 

35 

38 

41 

38 

35 

26 

59 

58 

51 

37 

T8 

29 

35 

85 

40 

38 

50 

8 

43 

39 

*9 

43 

39 

43 

22 

50 

45 

21 

50 

53 

41 

57 

42 

58 

52 

17 

30 

75 

31 

30 

23 

39 

57 

22 

35 

33 

43 

39 

24 

90 

25 

34 

68 

47 

59 

43 

3 

62 

31 

48 

87 

98 

30 

38 

53 

48 

48 

33 

27 

27 

38 

58 

54 

80 

23 

41 

26 

34 

42 

20 

27 

26 

42 

59 

73 

88 


there is of significance to him about the sequence in Fig. la. Well, 
the answer is that there is something about that sequence that is 
of very great importance to him. In the first place, most sequences 
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FIG 1b SEQUENCE OF VALUES OF THICKNESS OF INLAY ON RELAY SPRINGS 


obtained by a random operation possess certain characteristics 
that almost no sequence of results from repetitive engineering 
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operations is found to possess until after assignable causes have 
been eliminated through the application of the operation of sta- 
tistical control. Furthermore, if anyone not an experienced statis- 
tician were to try to write down a lot of sequences of 144 different 
numbers or if he were to try to arrange the 144 numbers in Table i 
in what he would instinctively call random order, most of these 
sequences would fail to possess the characteristics possessed by al- 
most all of the class of sequences obtained by drawing the 144 
numbers again and again at random as described above.^ Many 
of my colleagues have told me that their first real feeling for the 
meaning of a random operation came after they had tried to 
juggle with their eyes open, as it were, a set of numbers into what 
they would instinctively call random order only to find that the 
sequences thus obtained did not possess the characteristics pos- 
sessed by most of the sequences obtained by drawing from a bowl 
with their eyes shut.® Before we can go further in characterizing 
the quantitative differences between such sequences, we must 
consider first the problem of testing a statistical hypothesis and 
later that of testing the hypothesis that an operation is in a state 
of statistical control. 

Basic Test of Statistical Hypothesis. To determine whether any 
operation such as that of drawing from a bowl, as illustrated by 
the data in Fig. la, gives a sequence with the characteristics of 
one defined as random by the mathematician, we may choose 
one or more of the indefinitely large number of criteria that have 
been or may be established mathematically. If the operation gives 
a sequence that fails to meet the chosen criteria, the hypothesis 
is rejected, but if it meets the criteria, the hypothesis is accepted. 
Needless to say, no experimental test will prove or disprove a 

® Since there are different possible orders in which N numbers may be drawn, 
one might argue that any arrangement whatsoever that one chooses to make of 
the N numbers is a random arrangement in that it would be one of the possible 
orders obtained by a random operation of drawing. The important point to note 
is that I do not speak of a random number or a random arrangement except in 
the sense of a number or arrangement given by a random operation. What I am 
contrasting is the class of arrangements given by the oi>eration of arranging the 
numbers in what one intuitively may feel is a random manner with the class of 
arrangements given by the random operation of drawing from a bowl. For a more 
comprehensive discussion of this point cf. Shewhart, op. at. 

® The reader may wish to try for himself different orders of the data in Table i 
and see if his trials pass the three criteria considered later in discussion of statistical 
control theory. ' 
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statistical hypothesis, and in fact the test thereof constitutes a rule 
of behavior that must be justified upon the basis of extensive experi- 
ence showing that in the long run we shall not be too often wrong. 
In fact, the test of any statistical hypothesis is subject to the follow- 
ing two kinds of errors known as errors of the first and second 
kinds: (i) sometimes the hypothesis will be rejected even though 
true and (2) sometimes the hypothesis will be accepted even 
though false.^ 

Experience reported in the literature from many different 
sources justifies the conclusion that we may use with confidence 
the deductive distribution theory of the statistician to predict the 
distribution of any observed statistic of samples of size n or of any 
one of the many characteristics, such as lengths of runs-up and 
runs-down, of the infinite sequence of numbers that we may expect 
to get by repeating again and again without limit the operation 
of drawing a number from a bowl. Certain other operations as, 
for example, the use of tables of random sampling numbers also 
give results that have been found by experiment to possess the 
properties predicted by the mathematical statistician. 

Such studies show that a few specific kinds of operation exhibit 
properties described by the mathematical statistician as random. 
The ability to randomize a set of numbers or a set of objects by 
means of some distinguishable physical operation provides the 
scientist with a powerful technique for making valid predictions, 
and we shall now see how this can be used by the engineer. 

Four Uses of the Basic Contributions of Classical Statistical Theory, 
(i) Obviously the method of testing a statistical hypothesis can 
be used for testing the hypothesis that a sample assumed to be 
random came from an assumed law of chance, or it can be used 
to test the hypothesis that two samples, both of which are random, 
came from the same law of chance. Such is the nature of statistical 
tests of significance of observed differences between two or more 

^J. Neyman and E. S. Pearson have contributed many important papers on 
testing statistical hypotheses starting with one in Biometnkay July 1928, “On the 
Use and Interpretation of Certain Test Criteria for Purposes of Statistical Infer- 
ence.” Some of their latest contributions are given in Statistical Research Memoirs^ 
Vol. I, 1936, and Vol. 2, 1938, Cambridge University Press, London. Also see 
“On the Problem of the Most Efficient Tests of Statistical Hypotheses,” by J. 
Neyman and E. S. Pearson, Philosophical Trans. Royal Society of London, Series A, 
Vol. 231, pp. 289-337. It is to be noted that all such tests depend upon the assump- 
tion that the sample used in testing is random. Later in testing the hypothesis of 
statistical control, we start with testing the hypothesis that the sample is random. 



io6 STATISTICAL METHODS IN ENGINEERING 

samples of data, but it should be noted that these tests are valid 
only if the samples tested are produced by random operations. 

(2) If one after another of the operations of the engineer, such 
as an operation of measurement of some physico-chemical prop- 
erty or that of producing a given kind of product, could be shown 
to satisfy the hypothesis that it obeyed a law of chance,® then each 
such law and the estimates of the parameters therein could be quite 
simply obtained by means of well-established statistical proced- 
ures. Thereafter an engineer would be justified in using the dis- 
tribution theory of the statistician in predicting the outcome of 
any future repetitions of that operation with the same degree of 
assurance that he has in using statistical distribution theory in 
predicting the outcome of future drawings from a bowl universe. 
This constitutes a goal highly to be desired from the viewpoint of 
design, for then it would be possible to use mathematical distribu- 
tion theory in establishing the economical overall tolerance limits 
in terms of those for raw materials and pieceparts. 

However, as long ago as 1924, abundant evidence, since sub- 
stantiated on an even larger scale,® was obtained to show that 
few, if any, operations of the engineer obey laws of chance, even 
when carried out under presumably the same essential conditions. 
Two courses of action were open. One was to use available statis- 
tical technique as a curative measure in the sense of screening the 
product obtained by repeating a given production operation when 
this operation does not obey a known law of chance. The other 
was to develop a statistical technique to be used as a preventive 
measure in the sense of providing a means of detecting and elimi- 
nating assignable causes of variability that need not be left to 
chance. 

(3) Next, then, let us see how classical statistical theory pro- 
vides a means of screening the results already obtained by repeat- 
ing an operation, as in the production of a product of a given 
kind. The basis for the technique lies in the empirically estab- 
lished fact that the statistician has justified the use of certain opera- 

® Needless to say, in practice we can never be sure that an operation obeys a law 
of chance. 

* Cf. Economic Control of Quality of Manufactured Product^ W. A. Shewhart, New 
York, D. Van Nostrand and Company, 1931; Statistical Method from the Viewpoint 
of Quality Control^ W. A. Shewhart, Graduate School of the Department of Agricul- 
ture, Washington, 1939; and The Application of Statistical Methods to Industrial Stand’- 
ardimtion and Quality Control, E. S. Pearson, British Standards Institution, London, 
* 935 - 
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tions, like drawing from a bowl or from a set of random numbers, 
by means of which to randomize a set of objects. With this opera- 
tion of randomizing, one may break up any quantity of product 
into a number of lots and establish the most economical plan of 
sampling these lots to insure both that the producer’s risk of hav- 
ing a lot of satisfactory quality rejected will not exceed some 
previously specified value and that the consumer’s risk of accept- 
ing a defective lot will not exceed some previously specified value 
irrespective of the quality of the product in the different lots. 
These producer and consumer risks correspond to the errors of 
the two types always involved in testing a statistical hypothesis 
and were introduced into commercial use within the Bell System 
as early as 1925 in the development of sampling plans for screen- 
ing product.^® 

(4) Too much emphasis cannot be laid upon the practical im- 
portance of the fact that tests of statistical hypotheses are strictly 
valid only for random samples. This fact is taken into account in 
the design of sampling plans for screening, and it must also be 
taken into account in designing an experiment to test the signifi- 
cance of observed differences in the results obtained by submitting 
the results obtained by one operation to different subsidiary oper- 
ations. For example, a given kind of product may be subjected 
to laboratory or field tests under different conditions. Unless the 
operation of producing the product obeys a law of chance or, in 
other words, is in a state of statistical control, it is necessary to 
randomize the samples submitted to the subsidiary operations in 
order to obtain a valid test of the significance of the observed 
differences. If this is not done, observed differences between the 
results of two or more subsidiary operations may have arisen from 
assignable differences in the results of the first operation. Applica- 
tion of the operation of randomization is particularly important 
in the comparison of new designs, new materials or alloys, study 
of contact phenomena under diflFerent conditions, corrosion of 
materials under different atmospheric conditions, and field trials 
of equipment, to mention only a few.^^ 

10 Cf. “Method of Sampling Inspection,” by H. F. Dodge and H. G. Romig, 
The Bell System Technical Journal, Oct., 1929. 

The contributions of R. A. Fisher (see references to follow) are of particular 
importance in such studies. A simplified treatment of the elementary principles is 
given in Field Trials: Their Layout and Statistical Analysis, by John Wishart, School 
of Agriculture, Cambridge, England, 1940. 
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It is important for the engineer to keep in mind when reading 
all the literature on the randomization of the results of the first 
operation, that the validity of the tests for significant differences 
between the effects of different kinds of subsidiary operations rests 
upon the condition that the latter must be in a state of statistical 
control although this limitation is not explicitly stated. Hence, to 
be sure of the validity of tests of significant differences in the effects 
of subsidiary operations, we must first show that these operations 
are in a state of statistical control. This caution is necessary be- 
cause otherwise an engineer may accept the results obtained from 
small samples at their face value. 

By and large, classical theory treats of random fluctuations. It 
is a theory that tells how phenomena in nature would happen if 
they happened at random as do the results of drawing from a 
bowl and certain molecular phenomena treated in kinetic theory 
and statistical mechanics. It tells us how to discover such laws of 
chance if they exist and how to use them for purposes of valid pre- 
diction when discovered. It tells us how to make valid test of the 
hypothesis that a random sample came from a given universe or 
two random samples came from universes that have certain char- 
acteristics in common. It also provides us with a very important 
experimental operation of drawing at random that can be used 
in drawing random samples from the results of repetitive opera- 
tions already made. 

It is not, however, a theory designed primarily to tell us whether 
or not observed phenomena happen at random; or how to attain 
a state of statistical control (or randomness) of the cause systems 
underlying the operations of the engineer if these are not already 
in such a state. Whereas classical statistical theory gives us a very 
useful curative operation of randomizing results already obtained 
by a repetitive operation, control theory attempts to give us a 
very useful preventive operation of modifying the cause system 
underlying a physical operation or production process until it be- 
comes random in the sense of being in a state of statistical control. 

BASIC CONTRIBUTION OF STATISTICAL 
CONTROL THEORY 

That an ounce of prevention is worth a pound of cure holds for 
the application of statistics. For example, if a manufacturing 
process can be made to produce a quality of product distributed 
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in accord with a law of chance, one can then establish once and 
for all with the requisite degree of assurance by means of a suffi- 
ciently large sample, the probability that the quality of a piece 
of product will lie within any previously specified tolerance limits 
Lx and L^, Now, it may be shown rigorously that if one knew this 
probability sampling of lots would not tell us anything more 
about uninspected portions of the lots than" we knew before we 
sampled them.12 Hence if we could attain this idealized condition, 
the necessity of sampling would be completely eliminated and 
there would be no need of applying a screening process. Moreover, 
if it can be shown under such conditions that it is not possible to 
modify the manufacturing process further by simply removing 
assignable causes, then it follows that we have also minimized the 
percent to be rejected because of failure to meet the tolerance re- 
quirements. Let us now see how we may approach this idealized 
condition with its associated advantages. 

Basic Hypothesis for Statistical Control Theory, Perhaps, in simplest 
terms, the fundamental hypothesis is that even though the occur- 
rence of an engineering operation exhibiting a state of statistical 
control is as rare as the proverbial hen’s tooth, it is feasible and 
often desirable to establish a scientific method of modifying an 
existing operation until it obeys a law of chance. We shall consider 
the hypothesis of control in three parts, the first of which is: 

Hypothesis Ha: The maximum attainable degree of validity of pre^ 
diction that an operation will give a value X lying within any prevu 
ously specified tolerance limits is that based upon the prior knowledge 
that the probability of this event is f or more generally upon the prior 
knowledge of the law of chance underlying the operation. 

This part of the hypothesis is in line with the definite abandonment 
of the causal laws of classical physics and chemistry in favor of an 
indeterministic theory that includes the idea of probability in the 
ultimate laws. It is also in line with current theories of knowledge 
of the world in that no way has yet been devised for arriving at 
certain knowledge. So far as the statistical theory of control is con- 
cerned, no special attempt is made to justify this part of the funda- 
mental hypothesis. Instead it is simply taken over from modern 
physics and modern logic. As such it represents the limiting knowl- 
edge to which an engineer may hope to attain in giving assurance 

For an amplification of this point, see the paper contributed to this symposium 
by Captain Simon. 
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that his engineering operations will give results lying within speci- 
fied tolerance ranges. 

Hypothesis Ilb: The maximum degree of attainable control of the 
cause system underlying any repetitive operation in the physical world 
is that wherein the system of causes produce effects in accord with a 
law of probability. 

Such a state of control has been termed a statistical state of control. 
The hypothesis that such a state represents the limit to which one 
may hope to go in controlling a given operation by finding and 
removing a few assignable causes of variation was originally sug- 
gested by the second law of thermodynamics and its interpretation 
in kinetic theory and the theory of statistical mechanics. In much 
the same way that entropy measures the degree of run-downness 
of a physical system at a given energy level, so the degree of ap- 
proach to a state of statistical control measures the run-downness 
of the cause system underlying a given operation. In much the 
same way that it would take a Maxwell Demon to reverse an other- 
wise irreversible process, so it may be shown by a study of chance 
cause systems that it would take the equivalent of this Demon to 
modify a cause system already in a state of statistical control with- 
out changing, as it were, the whole operation and hence the whole 
system of causes. 

On many occasions, I have heard an engineer say on reaching 
a state of statistical control of some operation that, without chang- 
ing the kind of operation, he was going to decrease the range of 
variation still further by simply asking the men performing the 
operation to take greater pains to reduce the variability, but I 
have not yet witnessed one case where such effort succeeded. Hence 
the engineer is pretty safe in taking the state of statistical control 
as a limit to which he may go in reducing the tolerance range for 
that particular operation. 

Hypothesis lie: It is assumed that some criterion or criteria may be 
found and methods developed for their application to the numbers ob- 
tained in a sequence of repetitions of any operation such that whenever 
a failure to meet the criterion or criteria is observed^ an assignable 
cause of variability in the results given by the operation may be dis- 
covered and removed from the operation. It is further assumed that^ by 
the removal of a comparatively small number of causes^ a state of sta- 
tistical control is approached where the results of repetitions of the 
operation behave in accord with a law of chance. 
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The development of an operation of statistical control and its use 
in justifying hypothesis lie is an empirical contribution of mass 
production because only in such a process would it be economi- 
cally feasible to make the wide range of trials of experimental 
techniques for the purpose of finding one that works satisfactorily. 

Let us now consider the requirements in the way of experi- 
mental data for attaining control. 

Basic Experimental Data for Attaining Control, Let us consider the 
second or experimental step in the scientific method of using con- 
trol theory. The crucial difference between the experimental 
technique in control work and that in applying the classical theory 
of statistics to provide a screen is that in control work we pay 
attention to the condition under which an operation gives a 
value whereas in the screening operation we ignore this con- 
dition; in control work, the experimentalist must keep his eyes 
very much open, but in screening product he must, as it were, 
keep them blindfolded; and in control work, interest centers in 
controlling the product not yet made through modification of the 
underlying cause system, but in the screening process, interest 
centers in the product already made. 

To symbolize this situation let us rewrite sequence (i) and at- 
tach to each value its associated symbol for condition C^. Also 
let us divide the sequence into two parts representing that already 
observed and that observable in the future. Then we have 


JSTi, Z2, . . . , J 


f., . . . , . . . , Zn, 


Cl C2 a C, Cn 

Past 


X’fi-i-l) * • • 5 • • • 


Cn+1 Cn^k ( 2 ) 

Future 


Present 


The screening process applies to the X's of the past and eliminates 
the information contained in the C’s through the operation of 
drawing samples at random. It also completely ignores the results 
of future repetitions. In contrast, the control process focuses atten- 
tion on both the X^s and the C’s of the past in the hope of detecting 
and removing assignable causes in the C’s of the past so that these 
causes will not enter into future repetitions of the operation. Of 
coitrse, in practice the screening and control processes may be 
carried on simultaneously and the data obtained in the screening 
process may also be used for the purpose of detecting lack of con- 
trol. For example, inspection samples of lots taken in the order of 
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their production have long been used in some industries in apply- 
ing the operation of statistical control shortly to be described. 

Just as the batter must keep his eye on the ball first to make a 
hit and afterwards to determine what to do, so does the experi- 
mentalist have to keep his eye on the condition underlying the per- 
formance of any operation first to be able to apply this information 
most effectively in testing the hypothesis that he can do something 
to change the condition, and afterwards to determine what to do if 
the test is positive. By keeping his eye on the condition associated 
with each of n repetitions of an operation, the engineer may dis- 
tinguish the following three situations: 

( 1 ) In the absence of any a priori reason for distinguishing be- 
tween any two conditions, he may judge them to be essentially 
the same as may be symbolized by the equivalence 

C. C, (3) 

Operations of drawing with replacement from a bowl illustrate a 
situation where the experimentalist is practically forced to con- 
clude that the conditions are essentially the same from drawing 
to drawing. However, even under such circumstances, it is possible 
to arrange the results obtained by repeating an operation again 
and again in the same order as the operation was repeated. 

Now, if the repetitions of the operation do not obey a law of 
chance or, in other words, do not arise under a state of statistical 
control even though the experimentalist considers them to be 
essentially the same, as symbolized by (3), the sequence of values 
of Z, when taken in the order that they were observed, is not likely 
to pass the criteria of randomness established by the statistician. 
In fact, the order of observation in such instances provides the 
only quantitative basis for testing the hypothesis that the observa- 
tions came from a state of statistical control, and in practice it is 
usually found that this order indicates lack of randomness. Hence 
the statistical control engineer should insist that the record of the 
order of repetitions be preserved even though the experimentalist 
judges the conditions to be essentially the same. 

(2) The engineer may have a priori reasons for believing that 
the conditions do not remain essentially the same, as can be sym- 
bolized by 

C* C,. (4) 

For example, he may surmise that there are erratic effects or pos- 
sibly trends superimposed upon the effects of a law of chance. In 
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such a situation, the experimentalist may be in a position to 
suggest certain ways of ordering the results of n repetitions of an 
operation, independently of the magnitudes of the associated 
values of Jf, and solely upon his knowledge of the conditions under 
which the observations were taken, in order to reveal, if present, 
trends or erratic eflFects that he considers likely to exist. 

For example, in some recent work, interest 4:entered in attaining 
a state of statistical control of the variation in thickness of a rolled 
inlay of contact metal on a particular kind of relay spring. These 
springs were cut from a long strip and there were a priori reasons 
for expecting both trends and erratic effects in the thickness along 
the strip. Hence the springs were numbered in the order that they 
were cut from the original strip so that the measurements of thick- 
ness on the individual springs could be arranged in the order of 
their original position in the strip. These observations have already 
been given in Table i. The order in which the observations were 
arranged along the strip is that of the table, beginning at the left 
and reading down the columns. The set of 144 ordered measure- 
ments is given in Fig. ib and shortly we shall see how useful this 
experimental order is in giving clues to the presence of assignable 
causes .13 Hence, because of the importance of order, the statistical 
control engineer must insist that the experimentalist suggest orders 
in which to arrange the results of a series of n repetitions when he 
does not have any a priori reason for believing that any two repe- 
titions have been carried out under the same essential conditions. 

(3) The engineer may have a priori reasons for believing that 
the conditions may be divided into rational groups such that 
within each group the conditions are essentially the same but such 
that the conditions for each group are not essentially the same as 
those for any other group. In this situation, the engineer must 
insist that the experimentalist indicate the groupings and also the 
observed order within each group so that this order as well as the 
differences between the groups may be tested for indications of 
assignable causes. 

Thus far we have considered three kinds of information that the 
statistical control engineer needs to know about the conditions 
associated with the repetitions of an operation so as to make the 
quantitative results useful in testing the hypothesis that the opera- 

I am indebted to my colleague, Mr. E. B. Ferrell, for permission to use these 
data and for many helpful suggestions and criticisms about the scope of applications 
of a statistically scientific method as treated in this paper. 
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tion is in a state of statistical control. However, the contribution 
of the statistician does not necessarily end with the attainment of 
such a state, because even then the variance and hence the eco- 
nomic tolerance range may exceed the desired value. 

Accordingly, it may be desirable to change the whole operation 
or some part thereof. For example, even though the variations in 
the thickness of the inlay of contact material shown in Fig. ib were 
found to have risen from a state of statistical control, the variance 
might still be so large that many units of product would have to 
be rejected in assembly. Under such conditions, there would be 
no use trying to find and remove assignable causes of variability 
because these would have presumably been removed. If because 
of too large variance, the whole operation of producing the inlay 
is to be changed, the job is primarily one for the engineer, but if 
only some of the component parts of the operation such as the 
kind of inlay or fuse metal used, the rolling process, method of in- 
serting the inlay, or heat treatment of the material are to be 
changed in order to reduce the overall variance, the statistician 
may be of great service in showing how to obtain at minimum cost 
the data necessary for analyzing the total variance of the opera- 
tions into the component parts associated with the component 
parts of the operation of production. 

If some one of the component variances is significantly larger 
than the others, then it may be most economical to change, if 
possible, the corresponding component of the operation of pro- 
duction. In any case, the technique of analysis of variance makes 
it possible to estimate the maximum reduction in total variance 
that may be expected as a result of changing any component of 
an operation in a state of statistical control even though no sub- 
stitute for this component is known at the time. It is very impor- 
tant, however, for the engineer to note that, in accord with 
hypothesis Ilb, the problem of reducing the variance in the results 
of an operation already in a state of statistical control by changing 
one or more of the component operations is a fundamentally 
different^® problem from that of eliminating assignable causes 

^*The principles underlying the design of efficient experiments for analyzing 
variance have been set forth by R. A. Fisher and his co-workers. See, for example, 
Statistical Methods for Research Workers, 7th ed., 1938, and particularly The Design 
of Experiments, 2nd ed., 1937, Oliver and Boyd, London. Also see The Methods of 
Statistics, by L. H. G. Tippett, 2nd ed., Williams and Norgate, London, 1937, 
for applications in the cotton industry. 

This significant practical difference has, in general, been overlooked by many 
students of statistical theory. 
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without changing the component part of the operation such as 
kind of inlay used, the rolling process, and the like in the example 
considered above. 

When assignable causes are present, the engineer may reason- 
ably expect to find and eliminate these without changing the pro- 
duction operation or any of its component parts. For example, if 
assignable causes of variation are present in the process of produc- 
ing the thickness of inlay on relay springs as shown in Fig. ib, the 
engineer may find that they could be removed by further adjust- 
ment of the controls of component operations of rolling, heat 
treating, and the like. When assignable causes are shown not to 
be present, and it is found that one of the component operations, 
say the rolling operation, contributes a large share of the variance, 
it will usually be most economical to secure a reduction in variance 
by changing this component operation, perhaps by using another 
type of rolling mill. 

Enough has been said, I hope, to indicate some of the things 
that an experimentalist must do in taking and recording data if 
he is to obtain data that can be used efficiently in testing the hy- 
pothesis of statistical control and in analyzing the total variance 
under controlled conditions into its component parts or, in other 
words, if he is to make progress in the removal of assignable causes 
of variability that need not be left to chance in almost every field 
of science and engineering, thereby extending the potential use- 
fulness of raw and fabricated materials. 

Even in the face of this situation, however, it is still customary 
engineering practice to neglect the importance of order and in- 
stead group together into a frequency distribution all data whether 
C^^Cj or Ci4=-Cj and irrespective of whether or not subgroups of 
conditions might be suggested. Then most likely engineers will 
use the average of this distribution and in addition they may pos- 
sibly plot the distribution as a frequency curve or ogive. 

For example, in making field studies on new materials or de- 
signs, studies of the electrical properties of contact materials, or 
the corrosion of materials, the engineer’s faith is likely to be pinned 
on the averages of large numbers of observations whereas such a 
procedure is almost sure to mask the very differences that he is 
looking for. In fact, we shall shortly see that where there are 
assignable causes present, as soon as we group data together we 
are almost certain to destroy thereby all clues to the presence of 
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these causes, the effects of which must be eliminated before we can 
make valid comparisons of materials, designs, and the like. 

Basic Test of Statistical Control Hypothesis. The test of the hy- 
pothesis that a technique can be established for approaching a 
state of statistical control must be empirical. One sets up criteria of 
control and then searches for and eliminates, if possible, the causes 
producing the deviations that fall outside the limits fixed by the 
criteria. If such search reveals the presence of assignable causes 
and if, as these are removed, one approaches a state where the 
criteria of control are satisfied, we accept this evidence as an em- 
pirical justification of the hypothesis. 

A satisfactory technique for attaining a state of statistical con- 
trol has been described in considerable detail elsewhere^® and con- 
sists of the following five essential steps that are referred to in 
control theory as the operation of statistical control: 

(1) Specify in a general way how an observed sequence of n 
data is to be examined for clues to the existence of assignable 
causes of variability. For example, it is essential that the order in 
an observed sequence always be tested for randomness whether 

or 

(2) Specify how the original data are to be taken and how they 
are to be broken up into subsamples upon the basis of human 
judgments about whether the conditions under which the data 
were taken were essentially the same or not. 

(3) Specify the criterion of control that is to be used and indi- 
cate what statistics are to be computed for each subsample and 
how these are to be used in computing action or control limits for 
each statistic for which the control criterion is to be constructed. 
Three of the conditions that such criteria should satisfy are as 
follows: the limits in the criteria should be as nearly independent 
as possible of the functional form of the law of chance when the 
state of statistical control is attained; the criteria should in so far 
as possible minimize the error of accepting the hypothesis when 
false and should keep the error of rejecting the hypothesis when 
true less than some prescribed value fixed by economic considera- 
tions; and the criteria should indicate as closely as possible the 
condition under which the assignable causes enter the operation 
and as much as possible about the nature of these causes. 

(4) Specify the action that is to be taken when an observed 
statistic falls outside its control limits. The general action required 

Statistical Method from the Viewpoint of Quality Control^ loc. cit. 
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is to look for assignable causes whenever the criteria are not 
satisfied. 

(5) Specify the quantity of data that must be available and 
found to satisfy the criterion of control before the engineer is to 
act as though he had attained a state of statistical control. 

To illustrate the role played by statistical criteria in the process 
of detecting and removing assignable causes of variation and to 
show how important it is not to group data together so as to de- 
stroy all possibility of ordering them in terms of the conditions 
under which they were taken, let us see what happens when cer- 
tain of these criteria are applied to the two sequences of Fig. i. 
Will the criteria tell us anything about these two sequences that 
we cannot see by just looking at the sequences themselves? To 
comprehend the significance of order, it should be kept in mind 
that the data in Fig. la are the same as in Fig. ib except that in 
Fig. I a they appear in the order drawn from a bowl. Anything 
that we can find out about the presence of assignable causes by 
studying the order in Fig. ib is thus lost just as soon as this original 
order is destroyed by grouping the data together. Will the criteria 
of control indicate the presence of assignable causes of variation 
in both Fig. la and ib? Of course, we know that if they give an in- 
dication of such causes for the drawings from a bowl (Fig. la), 
this will likely constitute a false lead because experience shows 
that we should not expect to find any assignable causes in the 
operation of drawing from a bowl, particularly when carried out 
by an experienced observer. 

As a starting point, let us apply the control chart criterion based 
upon averages of successive fours.^^ Fig. 2 shows what happens: 
No false lead is given for the drawings from the bowl, Fig. 2a, but 
the presence of assignable causes of variation in thickness is indi- 
cated, Fig. 2b. 

For a detailed statement of the method of constructing a control chart, see 
Economic Control of Quality of Manufactured Product^ loc. cit.^ pp. 309-313. It should 
be kept in mind that, for reasons that I have given in the literature, it is desirable 
to use small subgroups as is here done. See, for example, Chapter i of Statistical 
Method from the Viewpoint of Quality Control, loc, cit. Of course, the practical man 
wants to know what is a small sample, and to him it may be of interest to know 
that in much of my own work, particularly in laboratory research, I have found 
it desirable to use where possible subgroups of size four. One reason for not using 
a larger sample is that, as we shall soon see, runs of seven or raoi e are very unlikely 
for a statistically controlled process. Hence, if we are to be sure that at least one 
subgroup will be completely within a run of seven, we should not use a subgroup 
larger than four. 
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Next let us see if other criteria may be found that will tell us 
something about the nature of assignable causes to be expected. 
For example, do you detect any evidence of causes producing 



FIG 2a CONTROL CHART FOR AVERAGES OF 



SAMPLE NUMBER 

FIG 2b CONTROL CHART FOR AVERAGES OF FOUR 
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trends in the data of Fig, ib? Is there a downward trend at the 
left and an upward one at the right? Or would you expect to find 
that the assignable causes are of the type that produce discontinu- 
ous and erratic shifts in the expected values? Is there any evidence 
for believing that both kinds of causes are present? It should be 
kept in mind that we seek answers to these questions so that we 
may be better able to detect and remove the assignable causes. 
Recent studies covering a broad field of research problems have 
shown that two very simple criteria may be used successfully in 
helping to distinguish between causes producing trends and those 
producing discontinuous and erratic effects. A complete report on 
the application of such criteria cannot be given here and it must 
suffice to indicate in broad outline the simple nature of these 
criteria and how they have been found to work. Basically, what is 
done is to note two kinds of runs in any sequence: runs up and 
down and runs of numbers above and below the average of the 
sequence where runs up and down are defined as follows. In any 
sequence, certain numbers are greater than either of their im- 



SCIENCE OF ENGINEERING 


1*9 


mediate neighbors and form, as it were, maxima in the sequence. 
In a similar way, we have minima. A run-down is the interval 
between a maximal number and the next succeeding minimal 
one, and the length of run-down has been defined in the literature 
as the number of numbers involved, including the maximal and 
minimal numbers. Correspondii^g definitions hold for runs-up. 
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FIG 3b SEQUENCE OF VALUES OF THICKNESS OF INLAY ON RELAY SPRINGS 

In Fig. 3 where the data of Fig. i are replotted, runs-up are shown 
by solid lines connecting observed values, and runs-down by 
dotted lines. Also in Fig. 3 the points above the average of the 
sequence are filled in and those below are not. 

W. O. Kermack and A. G. McKendrick^s have recently given 

“Tests for Randomness in a Series of Numerical Observations,” Froc. of Roy. 
Soc. of Edinburgh, Vol. LVII, pp. 228-240, 1937. If the total number of observed 
runs up and down be N, these authors show that the expected number of lengths 
2, 3, 4 . . are 0.6250000A; 0.2750000^^; o.o79i667jV; o.oi726i9jV; o.oo305o6jV; 
0.0004547^^ and 0.0000587JV. See also R. A. Fisher’s note, “On the Random 
Sequence,” in the Quarterly Journal of the Royal Meteorological Society, July, 1926, 
p. 250. Of course, the probabilities here given apply to a frequency distribution 
of runs in an infinite sequence. By making use of results given by Kermack and 
McKendrick in a paper, “Tests for Randomness in a Series of Numerical Observa- 
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the probability of a run of any length I occurring in an observed 
set of JV* runs up and down from any infinite random sequence in- 
dependent of the law of distribution of the variable. Likewise, 

W. G. Cochran^^ has recently given a formula for deciding whether 
sequences of similar meteorological events, such as runs of con- 
secutive wet months, may be expected under a state of statistical 
control. If the result of an operation must be either Ei or and 
if we know the probability p that an operation will give the event 
Ely and the probability {i — p) == q that the operation will give 
the event E^y then Cochran’s work gives the expected number of 
runs of length r (counting runs of both Ei and E^) in any sequence 
of m repetitions of the operation. If now we call X the average of 
a sequence of m observed numbers and if we call event Ei the 
occurrence of a number in the sequence greater than X and event 
£2 the occurrence of a number less than X, we may use Cochran’s 
theory to compare the observed number of runs of any length r 
with the corresponding theoretical number to be expected upon 
the assumption that the observed fraction of the numbers greater 
than the average X is equal to the probability p that the under- 
lying cause system will give a value of X greater than the number 

X, considered simply as a number and not necessarily as the 
average.^® 

Now, if the assignable causes simply produce erratic shifts in 
the expected values, it is not likely that the distribution of lengths 
of runs-up and of runs-down will be much disturbed, but such 
shifts will tend to give extra long runs of numbers above and below 
average. However, if the assignable causes produce trends with 


tions,” Proceedings of the Royal Society of Edinburgh, Vol. LVII, Part III, No. 17, 
pages 352-356, it may be shown that there will not be appreciable differences in 
this frequency distribution so far as we are here concerned with its use even though 
we are dealing with an observed distribution of runs from a finite sequence of 
144 numbers. 

Cf. “An Extension of Gold’s Method of Examining the Apparent Persistence 
of One Type of Weather,” in Quarterly Journal of the Royal Meteorological Society, 
Vol. LXIV, pp. 631-634, 1938. The number of runs of length r of the events E\ 
and El out of m trials is 

/m,r = ■\’pq'‘) + (m-r-i) when i S r ^ {m-i) 

where ^ = i . 

Of course, one might arbitrarily choose any value Xi and any value p and use 
Cochran’s theory to compare the observed number of runs of length r both above 
and below Xi with the corresponding theoretical number to be expected upon 
the assumption that the probability of occurrence of a number greater than Xi is p. 
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slopes large in comparison with the fluctuations produced by the 
superimposed causes acting at random, the distribution of lengths 
of runs-up and runs-down will be disturbed. Likewise, the distri- 
bution of runs of numbers above and below the average of the 
sequence may be somewhat modified. 

Of course, we should expect close agreement between theory 
and experiment for the lengths of runs of both kinds in Fig. 3a, 
but since the control chart Fig. 2b gave evidence of the presence 
of assignable causes, we may expect discrepancies between theory 
and experiment for the lengths of runs in Fig. 3b. The data for 
such comparisons are given in Table 2. As expected, there is close 

TABLE 2 

Drawings from bowl Measurements of inlay thickness 


Runs above and Runs up and Runs above and Runs up and 
below average down below average down 


Length 

of 

runs 
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served 

number 

Ex- 

pected 

number 
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Ex- 
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I 

42 

36 



29 

36 



2 

18 

18 

65 

61 

13 

18 

53 

57 

3 

10 

9 

22 

27 

3 

9 

28 

25 

4 

4 

4 

10 

8 

5 

4 

7 

7 

5 

4 

2 

I 

2 

0 

2 

2 

2 

6 

0 

I 

0 

0 

4 

I 

I 

0 

7 

0 

I 

0 

0 

4 
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0 

0 

8 

0 

0 

0 

0 

I 

0 

0 

0 


agreement for the drawings from the bowl. For the sequence of 
thickness measurements, there is excellent agreement between 
theory and experiment for the distribution of lengths of runs-up 
and runs-down but not for lengths of runs above and below aver- 
age. These results constitute good evidence for believing that 
the assignable causes underlying the measurements of thickness 
in Fig. 3b do not produce trends such as might be produced by 
lack of symmetry in the rolls but simply produce discontinuous 
and erratic shifts in expected values of thickness such as might be 
produced by slippage at the cleavage planes of the crystals in the 
inlay. 

Not only does the application of these criteria to the observed 

Of course, the total number of expected runs above and below average will 
not usually be the same as the coi responding number of observed runs. 
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order indicate the presence of assignable causes and give a clue 
to their nature, but it also indicates when the assignable causes 
enter. For example, we see from the table that runs above and 
below the average of length greater than 7 are very unlikely, 22 
and all we need to do is to note where these runs occur, keeping 
in mind that such indications are subject to errors of the first and 
second kind as described above. 

Since all of this information about assignable causes is lost just 
as soon as we group the data together in the form of a frequency 
distribution, as is so often done by the engineer, we see how very 
important it is to pay attention to order in taking and recording 
data. For example, years of experience have shown that Nature 
almost never gives us a sequential order that passes even the con- 
trol chart criterion and is almost certain not to give us one that 
passes all three criteria and that the failure to meet such criteria 
almost always is traceable to an assignable cause. Yet, if we were 
to take these same data without reference to observed order as 
might be done by shaking them up, drawing them at random, 
and then applying the three criteria, almost all of them would 
slip through the net, thereby failing to indicate the presence of 
assignable causes that must be detected and removed in order 
to attain the advantages of a state of statistical control. 

Since, as noted above, experience has shown that one can al- 
most always find and remove assignable causes when indicated 
by such criteria and thereby approach a state of statistical control 
beyond which it is not possible to go except by some fundamental 
change in the operation itself, we have good grounds for accepting 
the basic hypothesis of statistical control theory as stated at the 
beginning of this section. 

SUMMARY OF POTENTIAL CONTRIBUTIONS OF STA- 
TISTICS TO THE SCIENCE OF ENGINEERING 

The basic contribution of statistics to the science of engineering 
is an improved scientific method to fit the world of probability in which we 
live. Classical theory contributes the hypothesis of a repetitive 
operation obeying a law of chance, the knowledge of which com- 
bined with the knowledge of formal mathematical distribution 
theory enables an engineer to make valid predictions of the out- 

**The expected number is actually 0.5682 instead of the rounded value of unity 
given in Table 2. 
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come of future repetitions of the operation. Statistical control 
theory contributes the hypothesis that it is humanly possible to 
remove assignable causes of variability in the repetitive operations 
of the engineer until such operations approach a state of maxi- 
mum control and obey laws of chance, a knowledge of which 
provides maximum assurance that the results of repeating an 
operation will fall within previously specified tolerance limits. To 
test these two hypotheses, statistical theory provides the necessary 
experimental techniques outstanding among which are (a) the 
operation of randomization and (b) the operation of statistical 
control. 

Broadly speaking, statistical theory treats of repetitive opera- 
tions and provides the engineer with a method of regulating such 
operations to his best interest. Fundamentally, the engineer’s job 
is to devise operations that, if carried out, will give results within 
previously specified tolerance limits. Sometimes the operation, like 
that of building a bridge, is to be carried out only once or at most 
a few times, and sometimes the operation, like that of mass pro- 
duction, is to be carried out an indefinitely large number of times. 
Inherently, all such operations are potentially repetitive, and in 
this sense, differ only in the number of repetitions carried out. It 
has long been recognized that one of the most revolutionary prin- 
ciples ever introduced into manufacturing was that of interchange- 
ability dating back at least to Eli Whitney in 1 798. The introduc- 
tion of that principle prompted the engineer to consider the 
advantages of introducing repetitive operations into production 
processes, and the contribution of statistics to engineering may be 
thought of as a means of maximizing the advantages to be attained 
by interchangeability. 

The basic contributions of statistics to scientific method make 
possible the attainment of the following objectives that are not 
otherwise attainable and that are of interest to all of us: 

(1) Even before a repetitive operation has reached a state of 
statistical control, the application of statistical theory makes pos- 
sible the establishment of sampling plans that will screen at mini- 
mum cost the output of such an operation so as to meet previously 
specified tolerance requirements and previously specified producer 
and consumer risks. 

(2) The use of statistical theory provides efficient experimental 
techniques based upon the operation of randomizing the results 
obtained from one operation or production process that is not in 
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a state of statistical control before submitting them to two or more 
subsidiary operations or treatments for the purpose of comparing 
the effects of these subsidiary operations, as in the case of field 
tests and the like. Such a procedure minimizes the chance of con- 
cluding that observed differences in the effects of the subsidiary 
operations are significant when in fact they came about because 
of assignable differences in the results of the first operation. 

(3) The operation of statistical control provides an experimental 
technique for minimizing tolerance ranges and maximizing the 
assurance that the product turned out by a given process will meet 
its tolerance requirements. Such an operation makes possible the 
most efficient use of limited quantities of raw material and pro- 
vides the maximum degree of refinement attainable by any pro- 
duction process. Preliminary studies indicate that the operation 
of statistical control also provides a useful technique for eliminating 
assignable causes of variability in certain kinds of human effort 
as, for example, typing and other forms of transcription. Both 
strategically and commercially, industrial groups and even na- 
tions often need every increment of efficiency in the use of limited 
quantities of raw materials and human effort that can be provided 
through the application of the operation of statistical control. 
Likewise they often need maximum refinement in quality through 
elimination of assignable causes, not only in pursuit of the arts of 
peace but also in time of war. As one example, the attainment of 
maximum homogeneity and hence minimum tolerance ranges in 
the properties of raw and fabricated materials may extend the 
potential carrying capacities of ships in the air and on the sea. 
Needless to say both the engineer and the consumer of the en- 
gineer’s or manufacturer’s products stand to gain through the 
increased assurance that the products will be found to meet their 
tolerance requirements. 

(4) The operation of statistical control provides a technique for 
modifying and coordinating the three fundamental steps in the 
process of mass production, namely, specification, manufkcturing, 
and inspection, so that the maximum number of pieces of product 
having a quality within specified tolerance limits can be turned 
out at given cost. It does this by showing how to minimize the 
cost of inspection and the cost of rejection. 

In conclusion it may be said that statistical theory plus mass 
production provides a means of maximizing our physical comforts 
in time of peace and our strategic factors in time of war. 
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LESLIE E. SIMON* 

INTRODUCTION 

All of you have heard of the blind man in a pitch-dark room try- 
ing to find a black cat that was not there. I believe that a man 
would be in much the same position if he tried to write a specifica- 
tion without statistics. I say this, knowing full well that quite 
workable specifications are written by men who have no knowledge 
whatever of formal statistics. However, these men are none the less 
users of statistics. Without attempting to credit all successful 
achievement to statistics, it is evident that by keen native judg- 
ment and observation, these men have indirectly attained satis- 
factory specifications. But the problem of dealing with a variable 
product is fundamentally statistical, and these men could have 
achieved their desired end much more economically, swiftly, and 
unerringly with the proper statistical tools. The distinction be- 
tween the use and non-use of statistics in this situation is perhaps 
similar to that between the two hunters, one of whom is armed 
with a spear, the other with a rifle. Both may succeed in eating, 
but one must strive for a precarious existence, while the other may 
live in relative comfort. So in writing good specifications: they 
may be eventually reached by the brute-force method of trial and 
error and hard common sense, but can be attained with greater 
ease and certainty through the use of the scientific aids offered by 
statistical theory. It is of the contribution made by this powerful 
auxiliary that I wish to speak. 

Several years ago I was given what I believe to be an unusual 
opportunity to observe the functioning of purchasing specifica- 

* Captain, Ordnance Department, United States Army; Executive Oflficer, Bal- 
listic Research Laboratory, Aberdeen Proving Ground, Maryland. 
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tions, their relationship to standards of quality, and the power of 
the statistical technique in ascertaining and defining quality with 
swiftness, clarity, and validity. 

Since this experience is of a military character I shall, for obvi- 
ous reasons, confine myself to broad scientific principles which are 
of general application both to the Ordnance problems and to 
many problems in industry. I may, however, define the experience 
by stating the general circumstances. 

It is necessary for the Army to maintain stocks of ammunition. 
Unfortunately ammunition is subject to deterioration in storage. 
It therefore follows that, after a period of years, investigations 
must be conducted to determine the nature and extent of the de- 
terioration, and lots of ammunition must be assigned to various 
grades, with a view to their ultimate replacement when their 
quality becomes too poor. 

In order to detect deterioration, one must know what the qual- 
ity was originally. It is but natural, therefore, to look to the pro- 
curement specification as at least an indication of the desired 
initial quality. Since there is frequently some doubt about the 
meaning or efficiency of the procurement specification, one supple- 
ments this knowledge by a critical examination of all the accept- 
ance tests, subsequent tests, and recorded service use which pertain 
to the ammunition. This procedure gives a keen insight into both 
the level of quality which was current at the time, and the effi- 
ciency of the procurement specification as an instrument for 
obtaining a product of the apparently intended level of quality. 
Finally, tests must be performed to determine (at least within 
limits) the existing level of quality after deterioration, if any, has 
taken place. Since the function testing of ammunition is both de- 
structive and costly, it is highly desirable that all possible informa- 
tion be extracted from such expensive test data. Statistical methods 
are therefore used in the analysis of test data. They provide reliable 
protection against the frailty of human judgment in estimating the 
quality of lots of ammunition from the results of tests of samples 
taken from those lots, and are a cogent factor in making judgments 
comparable when made from time to time and by person to person. 
This procedure throws much light on the operation of the statis- 
tical method and the results which might have been obtained, had 
these methods been in use at the time the ammunition was origi- 
nally accepted. At least some of these principles have now been 
applied to current procurement specifications. I should like to take 
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this opportunity to express my gratitude to Dr. Walter A. Shew- 
hart for his generous help and patient advice while I was devising 
the system of surveillance of ammunition which has now been in 
use for over two years. 

THE OLDER METHOD OF PROCUREMENT 

Some of the contributions of statistics to the development of 
specifications and standards of quality can best be shown by first 
analyzing the conditions which exist in the absence of the statis- 
tical method. 

The attainment of one’s professional end under these circum- 
stances is sometimes highly imperfect, slow, and costly. When it 
becomes necessary to procure a product in quantity there is gen- 
erally a considerable degree of incertitude both as to the level of 
quality which is necessary and sufficient for the purpose intended 
and as to the level which is practically attainable in manufacture. 
These levels, if clearly defined, would surely form the basis of a 
sound and economical standard of quality. It appears that the 
management frequently avoids the tedious problem of determining 
an economic standard of quality, of evaluating the statistical con- 
stants associated therewith, and of devising an inspection pro- 
cedure to insure to an appropriate degree the meeting of the stand- 
ard by what appears to be a clever and practical stratagem. They 
merely decide not to bother about either the standard or the sta- 
tistics (perhaps on the grounds that it is impractical or at least 
troublesome) , and instead write a good, stiff acceptance specifica- 
tion with a view to getting the very best obtainable. I think that 
sly old Aesop could have written a good fable about this attitude. 
It carries grave penalties. 

Of course, an acceptance specification calls for inspection. In- 
spection is always costly, and therefore some sampling scheme is 
generally adopted in the interest of economy. In the case of ammu- 
nition, sampling must suffice, since a functioning test of the entire 
product results in its total destruction. The specification therefore 
generally reduces to one or more requirements of the form, “a test 
sample of n shall be selected at random from each lot, batch, or 
consignment of the product, and not more than c shall fail.” The 
n used in most of these specifications is very small. Many of you 
have or will come in contact with this type of specification. I wish 
to show that this popular and well-known type of specification 
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offers almost no guarantee of the quality of product which will be 
accepted thereunder, and that the quality of product accepted 
under such conditions is practically that which happens to be 
offered for inspection. 

Suppose, for example, that a manufacturer submits a lot for 
acceptance of which the fraction p are good and the fraction q 
are failures, or more specifically q is the fraction defective. It is 
easy to show by the binomial distribution that such a lot has a 
definite probability, of passing the above specification where 


P6 = 


^ n! 
M{n-k)!k! 


pn-kqk 


It is thus evident that if manufacturer A makes a very good, but 
not quite perfect product, some of his lots will be rejected; and if 
manufacturer B presents only very bad lots, some of his lots will 
be accepted. In general, if a large number of lots, all having the 
fraction defective, q^ are presented for acceptance, then lOO 
of the lots will be accepted and loo (i — ^5)% of the lots will be 
rejected. The most discouraging factor, however, inheres not in 
the rejection of some good lots or in the acceptance of some bad 
lots but in the fact that the accepted lots in this case are no better 
than the rejected lots. Suppose that k manufacturers produce lots 
each of fraction effective pi, p2> • • • respectively, where pi to pj, 
represents a limited range which in practice may be of the order 
of 0.99 to 0.85. Each manufacturer will have his lots rejected with 
more or less frequency, but the accepted lots of a manufacturer 
will be no better than the rejected lots. Whereas no precise state- 
ment can be made as to the proportion of accepted product at 
any given quality level without a knowledge of the number of lots 
presented at each quality level, it is nevertheless obvious that the 
average of the accepted quality will differ little from the average 
of the presented quality, and that the limits are the same, even 
though one require that there shall be no defectives in a small 
sample of n. For the case of zero defectives in a sample of ten, and 
equal frequency of all kinds of lots from 0.99 to 0.85, the average 
quality levels are 0.92 (presented) and 0.94 (accepted). For a 
specification of not more than one defective in ten, they become 
0.92 and 0.93. It is therefore clearly evident that the quality of 
product accepted under the old type of specification is practically 
that which happens to be offered for inspection. It also appears 
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that the only way to prevent the acceptance of poor quality con- 
sists in not allowing it to be presented for inspection. More will be 
said of this subsequently. 

This unhappy situation cannot be corrected by any simple 
remedy, for the same law applies whatever the value of c. Suppose 
that one reduces c to zero, in an ^ffort to improve quality. What 
is the result? A plot of the binomials, which in surveillance is called 
the operating characteristic of the specification, readily shows that 
the frequency of rejections of lots of poor quality is somewhat in- 
creased; so is the frequency of rejections of lots just short of per- 
fection. The penalties fall on both the just and the unjust as one is 
merely subjected more mercilessly to the rigors of an unjust chance, 
and the other merely permitted to profit less from its vagaries. 
The unsatisfactory condition remains unaltered, and additional 
grounds for producer dissatisfaction may have been introduced. 
Quality becomes a function of margin of profit: the larger the 
margin of profit, the greater the number of rejections the manu- 
facturer can afford; and, consequently, if it costs more to produce 
good quality than poor quality, the more he can consciously lower 
his level of quality in the interest of making maximum profit. On 
the other hand, one does not get what he does not pay for, and 
too small a margin of profit may well be a deterrent to good qual- 
ity. It is conditions like this that must have inspired the doggerel: 

The rain falls alike 

On the just and unjust fellow, 

But mostly on the just because 
The unjust has the just’s umbrella. 

These observations become clear in the hindsight of inspecting 
the record of what actually occurred. They might well have been 
clear beforehand, had one applied no more complicated statistics 
than that of the binomial distribution. 

TWO BASIC FAULTS IN THE SPECIFICATION 

The analysis thus far, however, is wholly destructive. It neither 
points out the basic faults nor suggests their remedies. A little 
further use of engineering sense and statistical methods will ac- 
complish these constructive ends. 

In the critical examination of long records such as those en- 
countered in large contracts one can make some important obser- 
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vations. One can clearly see, for example, that manufacturer A 
made a very good product, yet a moderate number of his good 
lots were rejected because of choice of acceptance criteria. On the 
other hand, the whole record shows that manufacturer B made a 
perfectly abominable product. He suffered a greater frequency of 
rejections than Mr. A, but he certainly managed to get many of 
his poor lots accepted before he went bankrupt or was barred 
from business, or perhaps he survived despite his lack of merit. 
All of A’s lots should have been accepted: none of B’s. One basic 
fault then lies in judging each lot on its own merits on a small 
sample basis which even in the case of homogeneous lots is obvi- 
ously not an efficient basis for discrimination between good and 
bad lots. 

If then, a large sample is economically prohibitive, one must 
turn to some other source of information. This information is, in 
general, abundantly available. It exists in the past performance 
of the manufacturer. It is a part and parcel of a very valuable 
asset which is generally called experience and expert engineering 
judgment. When data are costly and scant, no scrap of pertinent 
information should be allowed to escape; and statistics supplies a 
competent means (a) of determining the validity of the prior 
knowledge and (b) of mathematically combining the accumulated 
knowledge with that supplied by the sample itself. 

The second basic fault is failure to make use of available knowl- 
edge other than that supplied by the sample, but the cogency of 
this criticism may not be apparent without a brief description of 
the mode by which the statistical method makes efficient use of 
this knowledge. 

The technique of checking the validity of accumulated prior 
knowledge and of combining it with that supplied by the sample 
is a very common-sense rather than a technical procedure. It re- 
quires for its operation only three ideas: (i) the concept of statis- 
tical uniformity, (2) the operation of an elementary statistical 
technique, and (3) recognition of a simple empirical law. 

Irrespective of homogeneity, one has almost certain (not just 
probable) knowledge of the product actually sampled. However, 
any inference drawn about the remainder of the lot inheres in 
some relationship, known or unknown, between sample and lot. 
Unless a product be homogeneous in the sense of being statistically 
uniform, inferences regarding the lot, based on the observation 
of samples taken from that lot, must be severely circumscribed. 
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since the reliability of sample quality as a witness of lot quality is 
unknown. 

On the other hand if (a) statistical uniformity exists, and (b) 
one knows the general level of a quality characteristic, and (c) the 
measure of its variation, then statistical theory can predict the 
limits within which practically all samples of any size, n, should 
lie. One can obtain the knowledge of the general level of a quality 
characteristic and the measure of its variation from the accumula- 
tion of a large number of small samples. One can then predict the 
statistical limits of sampling fluctuation for practically all samples 
of size n. If a sample fall outside these limits, it appears that a 
state of statistical uniformity does not exist. If, however, no sam- 
ples fall outside the limits, there is at least no reason for not be- 
lieving the product to be statistically uniform. This inference is 
dependent upon an empirical law. 

Extensive experience has shown that once a state of statistical 
uniformity is attained in a manufacturing process, the state of 
quality is inherent in the process itself, and cannot be changed 
without changing the process. 

With these three ideas: viz., statistical uniformity, a statistical 
technique, and an empirical law, and with these three pieces of 
knowledge: viz., a state of statistical uniformity, the quality level, 
and the measure of sampling fluctuations, the small sample be- 
comes a significant index of lot quality. The steps to significance 
are two. First, the state of statistical uniformity, level of quality, 
and measure of variation are inferred from the accumulation of 
samples. (This is the abundantly available information.) If statis- 
tical uniformity cannot be inferred from the data, then efficient 
discrimination between lots of good and bad quality by means of 
small samples is impossible. Second, if the quality be satisfactory, 
and if a subsequent sample fall within limits, there is no reason 
for believing that the satisfactory quality level has changed, but 
if it be not within limits, then the quality level appears to have 
changed from the known satisfactory level. (This inference comes 
from the information of the sample.) 

These are essentially the principles used in Dr. Shejvhart’s 
Criterion I for quality control in manufacture. Concisely then, 
one may say, that whereas one cannot judge, except most approxi- 
mately, the quality of the lot from the small sample, one can judge 
by statistical methods whether or not there is any reason for be- 
lieving that the lot is of a different quality level from its predeces- 
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sors; and knowing the quality level of the predecessors from ac- 
cumulated data, one can then infer the quality of the lot with 
considerable assurance. Hence, the second and corollary fault in 
the specification consists of not using hindsight as an intelligent 
guide to foresight. 

EVOLUTION OF A STANDARD OF QUALITY 

In the light of these observations, let us consider the r61e of sta- 
tistics in the development of a purchasing specification. One is in 
an ill position to write a specification unless one knows what one 
wants in the first place. Even if one knows what one wants, it does 
little good to specify it, unless there is at least some assurance that 
it can be met. Hence, that which is ideally desirable must be 
tempered by that which is practically obtainable. The only way 
to ascertain what is practically obtainable is through examination 
of records.^ 

From a long run of records one can readily determine the frac- 
tion defective, the average, the standard deviation, or other meas- 
ures of the principal quality characteristics of the product which 
should be expected under good manufacturing practice. In my 
own case I generally consider every recorded test of the product 
under consideration and often the tests of allied products. Without 
statistical methods a fair consideration of such mass data would 
be almost impossible. Furthermore, the statistical method is essen- 
tial as a means of inquiry into the homogeneity of the product, for 
without a state of statistical uniformity (as has just been pointed 
out), one is very limited in drawing conclusions regarding the lot 
based on the observation of samples. ^ If a homogeneous product 
of satisfactory quality level can be economically produced under 
good manufacturing practice, then one is in position to describe 
an economic standard of quality. Whereas the standard may or 
may not be the specification quality, it is nevertheless a desirable 
fundamental for writing the purchasing specification. 

^ When the product is new, resort must be made to the process described under 
“The Judgment of Statistical Control,” Chap. I, pp. 43-46, & Chap. II, pp. 50-79, 
Statistical Method^ W. A. Shewhart, Graduate School, Dept, of Agriculture, Wash- 
ington, D. C., 1939. 

* Chapter X, pp. 1 21— 144, et. seq., Bxoiwmic Control of Quality of Manufactured 
Producty W. A. Shewhart, D. Van Nostrand Co., New York, 1931. 
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CONTRIBUTION OF STATISTICS TO THE DESIGN 
SPECIFICATION 

It is impossible to specify completely a definite series of opera- 
tions which will certainly detect that a very simple product is or 
is not of standard quality. Some of the quajity characteristics are 
not capable of definite measurement, and after all there is a limit 
to detail. Hence, the specification of the intended quality must be 
limited to some chosen or principal quality characteristics of the 
product. This is termed the design specification, as it describes 
the design which it is intended that the specification shall obtain. 
Statistics makes an important contribution to the design specifica- 
tion as the intended or aimed-at values of the quality character- 
istics of the design, and their tolerance limits can be most concisely 
described by such statistical measures as the average and standard 
deviation. In surveillance work a quality description or design 
specification is prepared for each grade of ammunition of a given 
kind. Furthermore, it should be noted that a sound foundation 
for the statistics of the quality characteristics described in the de- 
sign is provided by the investigation of the record which was made 
in the interest of the standard of quality. 

CONTRIBUTION OF STATISTICS TO THE 
ACCEPTANCE SPECIFICATION 

The contribution of statistics to the acceptance specification is 
perhaps more marked. In the acceptance inspection, which is 
generally on a percentage basis, it is obviously impossible to de- 
termine with certainty that all of the product will meet the ideal 
of the design. Here again it is necessary to compromise with practi- 
cality and frankly admit that the logical mission of this specifica- 
tion is a mere statement of the quantity and kind of evidence 
which will be accepted as a satisfactory indication that the product 
will meet the standard of quality. It describes no indefinite ideals 
but confines itself to the delineation of a sequence of operations to 
be performed upon the product, the interpretation of the results 
of such operations, and a definition of the hmits within which such 
interpreted results must lie in order for the product to be accept- 
able. In surveillance this procedure is, of course, the criterion by 
which lots are admitted to the respective grades described by the 
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design specifications for those grades. The objectivity of the statis- 
tical method makes it ideally adapted to the role of drastically 
limiting the need for interpretations of requirements. The necessity 
of interpretations and vague estimates usually opens the field for 
differences of opinion, whereas the statistical method swifdy pro- 
duces the estimate of the measure under consideration in a manner 
superior to unaided judgment and without the bias always associ- 
ated with personal consideration. This is a great step toward 
placing both purchaser and vendor on solid ground of mutual 
trust and understanding. Furthermore, in the light of the statistical 
technique the probability of accepting poorer than standard qual- 
ity can be consciously minimized to any degree which appears to 
be economically desirable. One can do much to reduce the chance 
rejection of good lots and the chance acceptance of bad lots. With- 
out an exhaustive investigation into the reasons for each step, 
which would be somewhat tedious, let us see one way (a number 
of others could be suggested) that a few modifications in principle 
can be made in the popular type of specification in order for it to 
have these manifold advantages. 

A STATISTICALLY SOUND ACCEPTANCE 
SPECIFICATION 

a. From the first lot of a series k test samples of n shall be taken 
as nearly as practicable in order of manufacture, and the fraction 
defective, of each sample computed. The average of all ^’s shall 
not exceed the specification value, q\ and no q shall depart from 
the average by more than three times the standard deviation of q. 

b. From subsequent lots of a series, a test sample of n shall be 
taken at random and the fraction defective, q, computed. The 
observed q shall not depart from the average of all qh by more 
than three times the standard deviation of q (based on all accumu- 
lated samples) nor shall it cause the average of all to exceed 
the specification value, q\ 

Cn In event of a rejection, the next lot presented shall be con- 
sidered a first lot of a new series. 

Detailed steps in the specifications and provisions for retests 
have been omitted in the interest of focusing attention on statistical 
measures. 

Part a, which involves a relatively large sample has two objec- 
tives. First, to subdivide the lot into sub-groups with respect to 
time in order that undue variation in the lot from part to part can 
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be detected. Second, to provide a large sample from which quality 
can be estimated with a high degree of precision, thereby providing 
effective discrimination between the manufacturer who has an 
established process in which the inherent quality level is good; i.e., 
of specification quality or better, and the manufacturer who is 
poor and who would ride as a parasite on the sampling fluctuations 
provided by a kindly chance. Part a makes relatively efficient dis- 
crimination on a small sample basis subsequently possible. It may 
be noted that the difference produced by the relatively large 
sample of kn is only one of degree, but it is nevertheless an im- 
portant difference. 

Part J is a step designed merely to detect a change in quality. 
It should be noted that it is necesasry to take adverse action if 
there appears to be significant change in quality, even though the 
sample which is indicative of the change might be acceptable if 
the average quality level inherent in that manufacturer’s process 
had happened to be lower but still within specifications. This situ- 
ation must obtain since such an instance indicates either a change 
in average quality or a loss of statistical uniformity and hence re- 
sults in an unknown situation, and one can no longer validly use 
the accumulated inspection data in predicting from sample to lot. 

Part c is, of course, an avenue (though a somewhat painful one) 
for the manufacturer to reestablish himself. 

These measures provide for an adjustment in acceptance criteria 
whereby the probability of rejection of a product of specification 
quality can be made very remote without prohibitive sampling. 
However, in attaining this end, the door has not been opened to 
the acceptance of poor quality, for: (a) the manufacturer has to 
qualify first on a large sample basis, and (b) if he should relax his 
quality after qualification, he might succeed in passing one or 
two poor lots, but the probability of several poor lots becomes in- 
creasingly remote and, once caught, he has to start all over with 
a new first lot of a series. The operation of these principles on ac- 
cepted stocks not only shows their advantages in bold relief, but 
also shows the great gain that would be made if the influx of poor 
quality were checked at the time of its occurrence. 

SUMMARY 

Statistical methods have been invaluable in evolving standards 
of quality for ammunition; for specifying grades based on these 
standards; and for writing the specifications for admission of lots 
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to these grades. In the field of procurement specifications their 
use has not been as extensive, because the need has not been felt 
so keenly. Nevertheless, some progress has been made. Statistical 
methods have proved to be a powerful tool in the critical exami- 
nation of some ammunition specifications prior to final approval. 
Their use, either directly or indirectly, is almost essential in de- 
termining a reasonable and economic standard of quality through 
the method of comparing the quality desired with that which can 
be reasonably expected under good manufacturing practice. In 
like manner the statistical technique renders a valuable service in 
framing the acceptance specification. Through its use the quantity 
and kind of evidence which will be accepted as proof that the 
product will meet the standard of quality can be clearly expressed 
in a fair, unequivocal, and operationally verifiable way. 
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Law is a term used in three senses. The oldest meaning of the term 
among jurists refers to the body of authoritative precepts recog- 
nized or established as rules of conduct and of adjustment of rela- 
tions among men, to be applied by the tribunals and administra- 
tive agencies of a politically organized society in determining 
controversies. Another meaning refers to what is perhaps better 
called the legal order — the regime of regulating conduct and ad- 
justing relations by the orderly and systematic application of the 
force of a politically organized society. A third meaning refers to 
what is perhaps better called the judicial process, to which we 
must now add the administrative process — ^the process of deter- 
mining the controversies involved in the regulation of conduct and 
adjustment of relations by applying to those controversies the 
authoritative precepts according to an authoritative technique. 
Two of those three meanings describe processes. But it was not till 
the present century that men thought of the legal order in terms 
of a process. Kant thought of it as a condition — a condition in 
which the free will of each was limited only by the like free will of 
all others according to a universal law. As to the third meaning, 
jurists of the last century simply ignored the processes by which 
the legal order is maintained, and concentrated on the precepts 
by applying which the condition (as they took the legal order to 
be) is realized. If, however, the nineteenth-century jurists over- 
looked the processes and saw only a condition brought about by 
precepts, the most insistent type of jurists of today ignore the pre- 
cepts and apparatus and see only the processes, which, indeed, 
are two-thirds of the matter. For the legal order is a large-scale 
process, a specialized form of a still larger-scale process of social 
control, and the particular operations of the judicial process and 

* University Professor, formerly Dean of the Faculty of Law, Harvard University. 
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the administrative process are directed toward the ends of the 
large-scale processes and designed to achieve them. 

Legislation is a process of framing the precepts by which the 
operations of the judicial and administrative processes are to be 
guided. Thus whether we are concerned with law or with law- 
making we are for the most part dealing with processes and opera- 
tions, not, as men held in the last century, with conditions and 
rules. Here is the point of contact between the science of law, as 
a practical science, having to do with the attainment of ends by 
processes carried out with a minimum of friction and waste to 
plans which we call ideals, and the science of engineering, as an- 
other practical science having to do with attainment of ends by 
processes carried out with a minimum of friction and waste to 
plans and specifications usually in the form of blueprints. We 
must remember that “ideal” comes from a Greek word meaning 
picture, and is as much a picture as is the blueprint. Industrial 
engineering as a practical activity has a very significant analogy 
for the jurist. I have not hesitated to call the task of the jurist and 
lawmaker one of social engineering. Where the industrial engineer 
seeks to minimize, if he cannot eliminate, the element of chance 
and narrow the limits of expected deviation, so the juristic social 
engineer seeks to minimize, if he cannot eliminate, the personal 
element in the processes of the legal order and to keep expected 
deviation within closer bounds. 

With reference to the judicial and administrative processes, we 
speak of this side of social engineering (for it is not the whole) as 
a quest of objectivity. This quest has been pursued by jurists since 
the Greek philosophers began to speculate about social control 
twenty-five centuries ago. To the Greeks, the phenomena of 
physical nature, the regularity and predictability of the return of 
the seasons, the phases of the moon, the succession of day and 
night, were helpful as suggesting an ideal of human conduct and 
exercise of the powers of social control equally regular and pre- 
dictable. The phenomena of physical nature furnished a sort of 
blueprint. But they did not provide specifications for carrying 
out the plan. Indeed, for a long time jurists and practising lawyers 
did not get beyond the blueprints to the mode of controlling the 
operations so as to realize the blueprint plans as completely as 
possible. 

As Dr. Shewhart puts it, it is not merely how we aim but how 
we hit that counts. General Schofield tells a story of a volunteer 
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and a regular soldier who had fought side by side at the battle of 
Wilson’s Creek. The volunteer’s cartridge box was empty, and 
he asked the regular how many rounds he had fired. The regular 
looked at his cartridge box and answered “nineteen.” The volun- 
teer explained proudly that he had fired forty. “Yes,” said the 
regular, “how many rebels did you hit?” The answer came, “I 
don’t know; how many did jotf'"hit?” Said the regular, “I reckon 
just about nineteen.” Control of operations so as to make results 
conform as closely as possible to the ends sought is much the same 
problem in a factory, in an army, or in a legal order. 

For a time, objectivity was sought by attributing divine au- 
thority to the body of legal precepts, assuming that the author- 
ity behind them would sufficiently insure that they would be 
applied and their inherent objectivity would take care of the 
process of application. Later, a basis of authority was found in 
custom on the theory originally, I suppose, that long-continued 
usage showed that the gods approved of the customary precepts, 
and so it was safe to apply them and unsafe to risk others which 
the gods had not approved. Still later, the authority of custom 
was referred to the superior wisdom of an ideal past, and then 
authority for a time was held self-sufficient or was bolstered up 
by an authoritative philosophy. Next, for a time, reason was put 
as the basis of the precepts and of their application, and the an- 
alogy of geometry was relied on to show how applications could 
be deduced infallibly from authoritatively prescribed or discovered 
premises. In the last century, metaphysics and a metaphysical 
conception of history as the unfolding of an idea, and so of legal 
history as the unfolding of an idea of right or idea of freedom, 
were taken to show the path to objectivity in the legal order; and 
in the second half of the century such ideas began to give way to 
a positivist conception of laws of social development analogous 
to the laws of physical nature as they were then understood. In 
all these successive theories of law, objectivity of application is 
assumed to be involved in the nature of the legal precept, although 
the classical natural law, taken to rest upon reason, did think of 
application guided by reason as well as precepts expressing reason. 

In the last quarter of the nineteenth century, Jhering with his 
theory of end or purpose in law and social-utilitarian philosophy 
of law, not as something self-developing but as designedly framed 
and applied for social ends, started a new movement in jurispru- 
dence which is still felt, especially in America. In this theory of 
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law the process, and direction of the process consciously toward 
conscious purposes, was stressed, and legislation and administra- 
tion, which the historical jurists of the last century considered 
futile or negligible, found important places in the science of law. 
It is true this turning of attention to the processes, without more 
theory of guidance of the processes or measuring the conformity 
of their results to the assumed ends than acceptance, express or 
tacit, of the older idea that the precept somehow involved in itself 
the guarantee of objective application, led presently to a reaction 
from the whole idea of objectivity and assumption that in the 
nature of things it was not attainable. Different varieties of self- 
styled realism — for their use of the name “realism” is often more 
a boast than a description — and phenomenalism have grown up 
which urge on one ground or another — economic determinism or 
Freudian psychology, or Neo-Kantian relativism or all embracing 
skepticism — ^that all we have been trying to do and have believed 
we were doing in law since the classical Roman jurists, whose 
books have been quarries ever since, is mere illusion or self- 
deception or superstition or pious fraud. A recent reaction from 
the extreme views of realists and phenomenalists has led to a re- 
vival of scholasticism, of the official philosophy of the Catholic 
church, the philosophy of St. Thomas Aquinas, and neo-scholastics 
in many lands are turning back to a basis of philosophical-theological 
authority. 

Lawmaking and application of law have been guided by re- 
ceived ideals which are as much a part of the law as the precepts 
developed, interpreted and applied to their pattern. Indeed, ob- 
jectivity in the phenomena of the legal order is itself an ideal. It 
is always hard to attain ideals, and one may concede that no more 
than an approximation sufficiently close for practical purposes 
may be expected. Moreover, twentieth-century psychology and 
economics make the ideals of social control and so of the legal 
order appear much more difficult of attainment than we took 
them to be. Hence it is now the fashion to give up the ideals, to 
give up the age-long quest of an objective process, and yield to or 
even get behind the revival of political absolutism. But to give 
up in this way is to give up half of civilization. 

It will have been seen that at the moment the science of law is 
very much at large. There is no agreement about fundamentals 
or even whether there are fundamentals. But a type of realist has 
been insisting that somehow we may build better on a basis of 
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statistics. Little has come of this, however, and it will be instructive 
to see why, before I turn to the light thrown upon the possibility 
of control of operations of lawmaking, and application of law 
through the judicial process, by Dr. Shewhart's discussion of 
quality control. 

For one thing, the realists who have been urging that in some 
way or other properly kept and properly organized statistics of 
the operations of the judicial process would go far to help toward 
solution of problems of jurisprudence have had a faith, curious in 
those who are in general so skeptical of everything but their own 
skepticism, that a mass of figures have an intrinsic objective va- 
lidity, and that a bit of juristic writing is scientific if only it has a 
sufficient number of tables of figures. Matters which have been 
well known from long and general experience are regarded as 
requiring to be proved by such tables although the tables may 
very likely be constructed on the assumption of what they are 
called upon to prove. 

In the next place, the unreliability of our official statistics, na- 
tional and state, so far at least as they have to do with the opera- 
tions of the legal order, has stood in the way of any effective use 
of them. Indeed, it requires very little experience of trying to use 
them to convince one of the truth of Bill Nye’s saying that figures 
won’t lie but liars will figure. 

Bureaus seeking to make a showing in order to obtain increased 
appropriations, clerks set to the task of preparing reports and pro- 
ceeding to pad them with statistics for statistics’ sake, and well- 
meaning blunderers preparing tables with no clear idea of what 
is or may be significant in the materials before them, have con- 
tributed to discourage the conscientious judge or advocate or 
teacher who would like to know the facts underlying current 
assumptions as to American justice. There have come to be abun- 
dant official tables and masses of figures. But they are seldom 
useful except for the purpose for which they were prepared, 
namely, to make a showing. On inquiry they are likely to prove 
inaccurate, uncritical, contradictory, and not comparable with 
each other. For example, some statistics as to the prevalence of 
crime in different parts of the country published a few years ago 
showed a much greater amount of larceny in a rural agricultural 
community than in one of the old settled states, chiefly urban and 
industrial, and with a shifting population. The compiler had not 
noticed that in the former state by statute larceny bad been ex^ 
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tended to embrace embezzlement, obtaining money or property 
by false pretenses, and receiving stolen goods, while in the latter 
the narrow definition of the English common law still obtained 
and only a taking from the possession of another without his con- 
sent and with intent to convert to one’s own use was called larceny. 
Again, not many years ago different federal bureaus dealing with 
different aspects of federal enforcement of a certain statute pub- 
lished tables in which the same subject was treated on different 
bases so that the figures could not be made to agree. In some juris- 
dictions, tables have been drawn up so carelessly or so perfunc- 
torily that the numbers reported by one as sent to penal institu- 
tions and the other reporting the number received by the same 
institutions did not agree. Certainly it was not safe to infer from 
the serious discrepancy that a large number of those sent to prisons 
had escaped in transit. In one case within my personal knowledge 
where tables published by different bureaus in the same jurisdic- 
tion, which should have agreed in the totals, were quite out of 
accord, when attention was called to the matter and revision was 
promised and supposed to have been made, it turned out that the 
revision had brought the totals into accord by a simple process of 
wrong addition. A mean was struck between the discrepant totals 
and the arbitrary result was put at the foot of the respective col- 
umns. Nothing could be more discouraging than experience with 
statistics published by state bureaus and departments in this 
country. 

Of late, judicial councils in a number of states have been putting 
out some reliable statistics of the work of the courts, carefully 
compiled to an intelligent plan and directed to certain ends and 
so useful for purposes clearly manifest. But even when statistics of 
that sort become general, the difficulty will remain that the re- 
alists who have been urging a science of law based on statistics 
have expected too much of them. We can’t expect to use them to 
solve the fundamental problems of jurisprudence. They can’t give 
us a measure of values of competing claims or a criterion of justice 
or a theory of what we are seeking to bring about by means of the 
legal order. But it does not follow that we have no use for statis- 
tics. On the contrary, the stress which we now put upon the legal 
order as a process and upon the judicial process as a significant 
meaning of the term law, and along with the administrative 
process as entitled to a place in the front rank in the science of 
law, things which the realists have been prominent in urging, in- 
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dicates where statistical method is to be made use of. We must learn 
how to use statistics to control the quality of the output of the oper- 
ations by which the legal order is maintained and carried on. 

As long ago as 1823, Jeremy Bentham, the founder of the science 
of legislation in the English-speaking world, urged a ministry of 
justice as the means of improying legislative lawmaking. Since 
Bentham it has been urged by nftny of the leaders in movements 
for the improvement of the administration of justice, notably in 
recent years in the United States by the late Mr. Justice Cardozo. 
It has been urged also that such an institution could do much to 
improve the quality of the output of the judicial process. But those 
who have urged this necessary agency of improvement have al- 
ways been asked how we may be sure that a ministry of justice 
will be able to bring about the quality of lawmaking and of appli- 
cation of law which we desire. Answer has not been easy. It is now 
commonly asserted that judgments of quality can never be more 
than subjective opinions of no scientific value. When it is urged 
that a ministry of justice will give us properly kept and intelli- 
gently organized statistics of the operations of the legal order, an- 
swer is made that we have still to develop adequate technique and 
well worked out principles before we may be assured that our 
ministry of justice, if such things can be measured or controlled at 
all, will be able to pass objectively upon the quality of the legisla- 
tive and judicial product and enable us continuously to improve it. 

When Dr. Shewhart was good enough to send me his illuminat- 
ing paper “Some Aspects of Quality Control,” I saw at once its 
bearing on this question of a ministry of justice. Quality control 
in the legislative and judicial and administrative processes was 
what jurists had been looking for. “Verifiable operational quality” 
and “objective judgment of quality” were phrases of familiar 
sound putting from the standpoint of the industrial engineer prob- 
lems which had long troubled the social engineer in the fields of 
law and legislation. Very likely my enthusiastic response to Dr. 
Shewhart was the causa causans of my appearance here in the ca- 
pacity of an engineer. 

Let me endeavor to point out some of the possibilities of*an en- 
gineering theory of quality control to the work of a ministry of 
justice. 

It is significant at the outset that quality is not defined as some- 
thing absolute but as related to a process— * a succession of perceiv- 
able qualia, i.e., repetitions of attributes, associated with a previ- 
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ously conceived or specified set of operations. In the science of law 
we have to think of quality not as absolutely determined by ab- 
stract reasoning or by logical deduction from some unchallengeable 
metaphysically given datum, as in the legal philosophy of the past, 
but as relative to processes of making or finding precepts or of de- 
termination of controversies or to continuous demands which are 
analogous to processes. As to the lawmaking and determining proc- 
esses, the succession of qualia we seek are in the same way associ- 
ated with a previously conceived and specified set of operations. 
Even in the case of demands or claims, which make one of the most 
difficult problems of our science, there is a suggestive analogy in 
the types of quality laid down for the industrial engineer. If we are 
considering the quality of a demand or claim, type one has to do 
with its character independent of other demands or claims, e.g., 
the intensity of the demand. Type two has to do with its character 
as one of a group or body of claims or demands, e.g., individual, 
public, or social. Type three has to do with its character of being 
asserted generally, widely, by many persons or otherwise. Can 
these things be determined objectively? One should ask this in 
comparison with the practical method of the past, namely, deter- 
mination by experience of what will satisfy the most of the whole 
scheme of demands or claims with the least friction and waste. 
But more can be made of the method of quality control with refer- 
ence to the processes of the legal order. 

What can we derive from Dr. Shewhart’s theory and technique 
of quality control? In the judgment of quality, we are told, the 
first step is legislative in character; specification, to be compared 
with legislative or traditional fixing of the ends to be attained by 
a precept. The next step is one of judging whether a determination 
when made by applying the precept is likely to meet the specifica- 
tions and whether there is any evidence that the specifications 
should be changed so as better to adapt them to realizing the ends 
prescribed. For the judgment of quality we need not only specifi- 
cations of form but also some definite indication of the degree of 
assurance that we are to have, in rendering the judgment, that the 
quality of a particular application of a precept or judicial or ad- 
ministrative determination will meet its associated specification. 
There must be an element of uncertainty in fixing the specifica- 
tions. The function of fixing them is legislative and is tempered in 
practice by a judgment based on experience as to the attainable- 
ness of any specified goal. It will be seen here that I have put one 
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of the chief tasks of a ministry of justice in almost the very language 
of engineering quality control. Again, the task of judging of quality 
involves seeing to it that any failures to hit the mark prescribed 
are analyzed in such a way as to indicate whether or not it is de- 
sirable to modify the specifications. This is as much the job of a 
ministry of justice as it is of an industrial engineer. In the one case, 
too, as in the other it is important to determine how much infor- 
mation is to be collected and how it is to be collected. Incomplete 
specifications have to be interpreted and amended so as to make 
them more complete; the evidence as to defects of specification has 
to be interpreted; allowance has to be made for counterbalancing 
effects of other specifications and a running record kept of the 
evidence. Thus evidence may be found of assignable causes of 
deviation and the deviations may be traced back to lack of control 
of specified qualities or to the insufficiency of specified character- 
istics. This could have been written of a ministry of justice or judi- 
cial council as well as of an industrial engineer. 

So, too, of the method of making a quality specification objec- 
tively verifiable. First, the qualia to be produced are to be specified 
(and herein lies much of the quarrel between lawyers and the 
ambitious young advocates of administrative absolutism), and, 
second, the detailed technique of verification is to be prescribed. 
There is always an uncertainty that any set of operationally veri- 
fiable specifications constitutes the ideal set, and so the jurist has 
no occasion to feel that he cannot achieve a reasonable practical 
approximation any more than the engineer. 

Preparation for legislation under the influence of natural-law 
ideas consisted in putting different proposed precepts and those 
already promulgated in other jurisdictions side by side and choos- 
ing the one that, subjected to the scrutiny of reason, seemed most 
abstractly just. The program of sociological jurists calls for study 
of the actual social effects of legal institutions, legal precepts, and 
legal doctrines. Rationalist natural law insisted exclusively upon 
reason. Positivism insisted exclusively upon experience. What we 
seek is experience developed by reason and reason tempered by 
experience. But experience is not something we can grasp and 
organize by offhand individual recollection and observation. The 
merit of the theory of quality-control for the jurist’s purposes is 
that it directs statistics toward definite preconceived ends instead 
of expecting statistics gathered to no clear purpose to reveal the 
something for which they are to be used. 
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It is not that a ministry of justice or a judicial council can take 
Dr. Shewhart’s pamphlets and use them as manuals. But if we 
cannot use his precise methods for our purposes, for which they 
were not designed, we can use them to indicate to us how to go 
about setting up an objective method of our own. We can perceive 
what to require of a method and what we may expect it to do. 
We can be helped to do intelligently what we have been doing by 
a halting and sometimes blind empiricism. 
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